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ABSTRACT.  This  is  a  translation  from  the  French 
of  "Recherches  sur  la  Spectroscopie  par  Transfor¬ 
mation  de  Fourier”  by  Mme.  Janine  Connes,  which 
appeared  in  1961  in  Revue  d'Optioue,  Volume  40, 
pages  45-78,  116-140 171-190 ,  and  231-255.  The 
article  considers  the  use  of  the  Michelson  inter¬ 
ferometer  for  spectroscopy  in  the  /infrared  and 
treats  the  following:  Fourier  analysis  the 
interferogram  using  high-speed  digital  computers 
and  analog  methods-,  instrumental  errors'',  analysis 
of  noise  and  procedures  for  attaining  maximum  s/n, 
ratios,  specific  applications  to  night  sky  spec¬ 
trometry  in  the  lM  to  1.6m  region  as  well  as  to 
the'  r, e-emission  spectrum  of  c/ermanium,  and  com¬ 
parison  of  this  method  with  those  using  classical 
spectroscopic  instruments. 
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SPECTROSCOPIC  STUDIES  USING  FOURIER 
TRANSFORMATIONS 
by 

Mme .  Janine  Connes 
Labor at or ie  Aime  Cotton,  Bellevue 

*  * 

INTRODUCTION 

There  are  two  principal  ways  to  look  for  improvements  in 
spectroscopic  techniques  [l].  One  method  is  constructing  in-— - 
struments  to  receive  the  broadest  spectrum  possible  for  a  given 
resolution.  Towards  this  end  several  methods  have  been  pro¬ 
posed:  double  or  multiple  passage  of  a  beam  of  light  through 

a  single  dispersive  element  [2],  use  of  the  Fabry-Perot  inte¬ 
gral  spectrometer  [3j,  or  of  the  sisam  [4]*.  But  these  methods 
make  poor  use  of  the  energy  of  the  source  in  time  because  the 
spectral  elements  are  analyzed  successively.  A  fundamental 
improvement  can  be  obtained  by  receiving  energy  from  all  the 
spectral  elements  simultaneously;  this  is  what  is  done  using 
spectrographic  techniques  and  the  Fourier  transform  method. 

The  principle  of  spectrographic  methods  which  use  a  two- 
beam  interferometer  and  a  Fourier  transform  has  been  known  a 
long  time.  Michelson  conceived  the  first  in  1891.  He  showed 
that  the  intensity  I  of  the  light  that  can  be  measured  at  the 
exit  of  a  two-beam  interferometer  with  a  variable  phase  differ¬ 
ence  is  the  Fourier  cosine  transform  of  the  spectrum  of  the 
incident  light  [5,6 J.  Since  the  Fourier  transform  (F.T. )  is 
an  invertible  operation,  it  therefore  suffices  in  theory,  to 
reconstruct'  the  spectrum  studied,  to  make  the  cosine  transform 


■  ■  in  >  i  .  ■  ■  i  — •  1  1  1  — 

Translator's  Note:  As  defined  in  [4]  s.i.s.a.m.  is  an 
interference  spectrometer  for  selecting  the  modulation  ampli¬ 
tude. 
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of  the  interf erogram  function  representing  the  variations  of 
I  as  a  function  of  the  phase  difference  6  between  the  two  inter¬ 
fering  beams.  But  taking  the  F.T.  of  a  very  complex  function  is 
an  extremely  long  procedure  if  there  is  no  high-speed  modern 
computer  at  one’s  disposal,  and  Michelson  did  not  atvack  the 
study  of  an  arbitrary  and  extended  spectrum.  He  perfected  a 
method  adapted  to  a  very  particular  problem:  the  study  at 
very  high  resolution  of  an  extremely  narrow  spectral  domain, 
made  by  studying  the  curve  of  visibility  V,  the  envelope  of 
the  interf erogram  function.  Lord  Rayleigh  explained  how  the 
knowledge  of  V  alone  permits  constructing  the  spectral  profile 
of  the  line  studied  without  ambiguity  only  in  the  case  of  the 
symmetric  line  [ 7 ) .  This  method,  however,  permitted  Michelson 
to  exceed  the  highest  powers  of  resolution  attained  at  that 
tine  and  to  make  evident  the  hyperfine  structure. 

Nowadays  to  study  the  structure  of  a  line  at  very  high 
resolution,  one  uses  large  gratings  (8),  and  plane  [ 9 ]  or 
spherical  l 10]  interferometers.  However,  Michelson’ s  method 
retains  a  certain  interest,  and  recently  TeriMen  showed  how 
the  study  of  V  determined  by  photoelectric  observations  made 
under  the  best  conditions  resulted  in  information  about  the 
profile  oi  a  line  (known  a  priori  to  be  symmetric)  which  was 
■ore  precise  than  that  obtained  by  other  methods  [11,12]. 

The  possibility  of  using  photoelectric  photoconduct ive , 
or  thermal  dctc.tors  to  record  the  "interf erogram”  function  and 
the  means  now  at  one’s  disposal  for  making  the  F.T.  have  opened 
a  field  of  application  of  the  F.T.  method  that  Michelson  never 
suspected.  Rubens  ana. Wood  in  1910  obtained  the  first  spectra 
in  the  far  infrared  b>  taking  the  F.T.  of  an  interf erogram. 

Much  more  recently  P.  Fellgett  [13  to  15]  and  P*  Jacquinot 
[1,50,51]  independently  showed  that  the  reconstruction  of  a 
spectrum  by  the  Fourier  transformation  of  an  interf erogram  has 
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a  great  advantage  over  the  direct  recording  of  the  same  spec¬ 
trum:  one  can,  while  using  a  single  channel  detector,  measure 

Simultaneously  all  the  spectral  elements  as  with  spectrographlc 
methods--hence  the  name  given  by  P.  Fellgett  to  the  method, 
"multiplex  spectrometry".  A  simple  argument  allows  one  to 
predict  the  order  of  i  -^nitude  of  gain  that  can  be  expected. 

If  with  an  exploring  spectrometer  it  is  desired  to  measure  M 
spectral  elements  during  time  T,  each  element  is  observed  dur¬ 
ing  time  T/M;  with  a  spectrograph  or  a  multiplex  method,  each 
element  is  observed  during  the  whole  time  T.  In  the  ease  where 
the  fluctuations  have  a  spectrum  of  uniform  frequency  and  where 
the  background  noise  is  not  increased  by  the  simultaneous  ar¬ 
rival  of  all  spectral  elements,  the  precision  varies  as  the 
square  root  of  the  exploration  time.  The  gain  in  sensitivity 
is  therefore  of  the  order  of  y/M.  Hence,  in  the  infrared  where 
photon  noise  is  negligible  compared  with  receiver  noise,  the 
spectroscopic  method  using  Fourier  transforr  <>  'ons  is  the  most 
powerful,  no  matter  what  may  be  the  desired  resolutions  and  it- 
is  all  the  more  interesting  because  the  spectrum  studied  con¬ 
tains  more  spectral  elements.  This  gain  in  sensitivity  can  be 
expressed  in  several  ways: 

--the  signal/noise  ratio  can  be  multipled  by-v/H  if  the 
duration  of  the  measurement  is  the  same  as  in  the  classical 
method; 

--for  equal  signal/noise  ratio,  the  duration  of  the 
measurement  is  divided  by  M; 

--rinaJi.,,  tor  a  given  s/n  ratio  and  measurement  time, 
one  can  obtain  an  inert  ase  in  resolution  and  thus  treat  some 
problems  inacceiii.Jt  by  any  other  method.. 

Furthermore,  P.  Jacquinot  llj  showed  independently  that 
in  addition  to  the  preceding  advantage,  the  use  of  a  Michelson 
type  interferometer  allows  him  to  employ  &  luminous  flux  which 
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was  far  greater  than  he  could  use  with  slit  devices ,  indeed 
as  much  as  that  given  by  a  Fabry-Perot  interferometer.  Finally, 
P.  Connes  [39]  showed  how  this  flux  could  be  increased  even 
more  if  one  inserted  afocal  systems  upon  each  of  the  inter¬ 
ferometer  arms. 

At  a  colloquium  on  interferometric  spectroscopy  in  Bel¬ 
levue  in  1957,  the  first  results  obtained  by  the  method  were 
collected.  It  should  be  noted  that  the  first  problems  treated 
were  astrophysical  problems  because  in  this  domain  there  exist 
some  feeble  sources',  the  study  of  which  is  interesting  even  at 
low  resolution.  Fellgett  obtained  the  first  star  spectra  in 
the  near  infrared  at  a  resolution  varying  from  60  to  100,  by 
using  an  interferometer  with  reflecting  trihedrals  [16], 

Strong  and  his  collaborators  Vanasse  and  Gebbie  Investigated 
the  far  infrared  region  [17  to  19].  For  an  interferometer, 
they  used  a  grating  with  lines  of  variable  depth  which  did  not 
’have  symmetry  of  revolution,  but  had  the  advantage  of  not  re¬ 
quiring  separating  and  compensating  plates  (‘difficult  to  effect 
in  this  spectral  region).  Gebbie,  with  his  instrument,  measured 
for  the  first  time  the  solar  spectrum  in  the  region  of  300  to 
1,000  ft  with  a  limit  of  resolution  of  0.2  cm*1  and  localised 
the  transmission  window  with  precision  [20]. 

In  the  visible  range  Mertz  constructed  polarization  inter¬ 
ferometer  consisting  essentially  of  a  variable  retardation 
plate  between  polarisers  [21];  the  maximum  resolutions  attain¬ 
able  by  this  system  are  low,  of  the  order  of  100.  He  hoped  to 
make  an  empirical  classification  of  star  spectral  types  based 
on  interfercgrams  without  making  F.T.;  Bakhshien  seems  to  have 
thought  of  the  same  instrument  but  he  does  not  give  a  single 
detail  about  the  theory  nor  about  the  construction  of  his  ap¬ 
paratus  [22]. 

In  each  case  the  work  took  place  in  two  stages;  the  inter- 
ferogram  function  was  recorded  and  later  the  F.T.  was  calculated 
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by  digital  computers,  using  discrete  values  of  this  function. 
Strong  made  {23]  an  analogue  computer  giving  the  intensities 
of  10  spectral  elements  for  each  passage  of  the  moving  part, 
but  this  apparatus,  requiring  several  explorations  to  cover 
all  the  spectral  domain  studied,  does  not  allow  one  fully  to 
exploit  the  advantage  of  the  method.  Different  authors  have 
mentioned  the  fact  that  the  resolution  was  proportional  to 
the  maximum  phase  difference  attained  [1,15,11];  certain  ones 
have  studied  the  conditions  under  which  the  numerical  Fourier 
transformation  should  be  made;  but  they  limited  themselves  to 
the  study  of  the  case  where  the  spectrum  extends  from  frequency 
0  to  a  maxjuaum  frequency  [15,21,24]  and  have  not  studied  the 
conditions  under  which  the  numerical  transform  ought  to  be  made 
to  obtain  the  maximum  signal/noise  ratio  in  the  spectrum. 

In  Chapter  II  of  this  paper  a  more  complete  study  is  con¬ 
sidered;  one  determines  the  form  of  the  theoretical  equipment- 
function  in  the  case  where  the  Fourier  transformation  is  made 
by  a  numerical ^calculator  and  one  treats  the  number  of  points 

Ns. 

to  pick  from  the  interf erogram  as  a  function  of  the  spread  and 
position  of  the  spectral  domain  to  be  studied. 

In  Chapter  III  we  study  the  form  of  a  real  apparatus- 
function;  it  is  different  from  the  theoretic  apparatus-function 
for  several  reasons: 

--necessity  of  x'eceiving  on  the  detector  a  beam  of  finite 
width; 

--imperfections  of  the  optical  components  and  in  the  ad¬ 
justment  of  the  instrument; 

--irregularities  in  the  motion  of  the  moving  part; 

--error  in  the  adjusting  of  the  -e ro  phase  difference  of 
the  interf erogram. 

Chapter  IV  is  devoted  to  the  calculation  of  the  mean 
square  error  of  the  fluctuations  in  the  reproduced  spectrum. 
Chapter  V  to  the  study  of  the  different  possible  ways  to 
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obtain  In  a  given  experimental  time,  the  maximum"  s/n  ratio 
In  the  ease  where  the  F.T.  is  computed  numerically. 

Chapter  VI  sets  forth  the  practical  setting  up  and  re¬ 
sults  obtained  in  two  experiments:  the  spectrum  of  night 
sky  emission  and  the  spectrum  of  germanium  recombination  at 

1.6  it. 


Finally,  Chapter  VII  is  giyen  over  to  an  attempt  to 
synthesize  the  spectroscopic  methods  called  nclassicn  and 
the  Fourier  transform  method;  there,  it  will  be  shown  how 
every  case  the  spectrum  is  obtained  by  making  the  cosine 
transform- of  the  autocorrelation  function  of  the  vibration 
of  incident  light. 


in 


II.  THEORETICAL  APPARATUS- FUNCTION 
NUMERICAL  FOURIER  TRANSFORM 


We  shall  recall  the  principle  of  the  method  before 
determining  the  relations  which  exist  between  the  reproduced 
spectrum  and  the  spectrum  studied.  Spectroscopists  call  the 
"apparatus-function”  the  distribution  of  the  spectral  density 
that  is  obtained  if  the  instrument  receives  strictly  monochro-. 
matic  radiation.  In  the  classical  methods,  the  spectrum  ob¬ 
tained  in  the  general  case  is  the  convolution  of  the  spectrum 
under  consideration  with  the  apparatus  function;  we  shall  show 
that  it  is  the  same  thing  in  the  r curler  transform  method. 

This  chapter  is  devoted  to  the  study  of  the  theoretical  appara¬ 
tus-function,  i.o. ,  the  instrument  i.  assumed  perfect  and  used 
with  a  beam  oj  zero  extent.  Likewi-.c  we  shall  treat  the  im¬ 
portant  problem.,  ©i  tne  number  of  .points,  to  choose  from  the 
interf erogram  i.i  oi'dci  to  make  a  numerical  Fourier  traiwfoSN*. 


1.  PRINCIPLE  OF  THE  METHOD 

Let  us  consider  a  two-beam  interferometer  in  which  one 


can  vary  the  phase  difference  be¬ 
tween  two  beams,  for  example  a 
Michel  son  inter! crometer.  For  the 
.at.  01  simplicity  wt  suppose  that 
the  beam  splitter  is  attached  to  a 
compensating  glass  plate  of  equal 
thickness  (Fig.  1>.  One  of  the 
mirrors  can  be-  displaced  parallel 


Fig.  1.  Diagram  a;  u.. 
Principal  01  a  Minn  1- 
son  Interferometer. 


co  o.t sc- II  with  rens'-ant  velocity  V;2. 
suppose  tb.e  interf eroiiK  ter  to  bo  il¬ 
luminated  by  a  strictly  monochromatic 


light,  of  luminance  B,  with  wave  number  n#.  -  The  beam  splitter 
divides  the  light  into  two  equal  parts  (if  it  is  not  absorbing 
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and  if -the  reflectivity  61  is  equal  to  Q.5)  which  it  recombines 
to  make  them  interfere  after  they  suffer  a  relative  delay 
6  *  2* ,  2*  being  the  difference  between  the  optical  paths  taken 
by  the  two  beams,  the  light  on  leaving  .follows  the  well-known 
law  •  . 

L  *  B  cos*  ■  (B/2)  (I  '+  cos  2re#d). 


It  is  composed  of  a  constant  part  and  a  variable  part  which,  . 
up  to  a  constant  multiplier,  is  given  by 


.with 


1(B)  «  B  cos  2*e#B  ■  B  cos  2»Vte#  »  B  cos 


%  “  v  V“  I*.  V/c. 


The  detector  gives  a  sinusoidal  electric  signal  with  amplitude 
proportional  to  the  luminance  received  and  frequency  equal  to 
the  product  of  the  frequency  of  the  light  N#  by  the  fraction* 
V/o„  For  example,  in  the  case  of  the  green  Mercu  line,  if  V 
equals  l(f*  duty's,  *  1.83  cps. 

If  the  interferometer  is  illuminated  by  several  monochro¬ 
matic  lines,  each  radiation  is  treated  in  the  same  way,  and 
the  resulting  signal  is  the  sum  of  the  signals  corresponding 
to  e  >ch  line.  In  the  case  of  an  arbitrary  spectrum  taken  Be* 
tween  and  c#,  the  interferogram  has  the  form 


or 


K«) 


B(0  cos  2rwB  <*r 


5(0 


B(r)  cos  2*vt  dm 
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according  to  the  variable  considered.  1(d)  and  3( t)  are  lumi- 
nances,  B(<r)  and  B(r)  densities  of  luminance  ( luminance/cm” 1 
i~„.-  the  first,  luminance/cps  for  the  second).  This  is  the 
Fourier  transform  of  the  spectrum  B(a)  [or  B(v)1 i -  In  the  case 
of  strictly  monochromatic  radiation,  the  interferogram  will  be 
a  sine  wave  of  constant  amplitude.  In  practice  the  lines  have 
a  finite  width  and  the  sine  wave  has  an  amplitude  which  de¬ 
creases  when  the  phase  difference  increases  (Fig.  2).  When 

the  spectrum  is  very  rich,  the 
different  sine  waves  correspond¬ 
ing  to  the  different  frequencies 
are  in  phase  only  for  zero  phase 
difference  and  the  interferogram 
appears  to  be  reduced  to  some 


fringes  in  the  neighborhood  of 


this  position  (Fig.  3).  Between 
these  two  extreme  cases,  all  inter* 
mediate  cases  are  possible,  the 
modulation  dropping  in  direct  pro¬ 
portion  to  the  increase  in  size  of 
the  spectral  interval  (Fig.  4). 

In  order  to  reconstruct  the 
spectrum  from  the  interfergram, 
it  is  necessary  to  make  tha  invars# 
F.T. 


44 m •Hum 


Fig.  2.  Interferogram  of 
the  Line  5,401  A  Produced 
by  a  Mercury  Lamp  Under 
Low  Pressure  [*d.d.m.w 
Phase  Difference!. 


n  -H® 

B(a)  «  j  1(B)  cos  2 rad  dd. 


pass  6i*  corresponding  to  the 
where  dr  ■  V  da  and  centered 


One  can  consider  making  an  harmonic 
analysis  of  the  signal  with  the  aid 
of  n  filters,  each  having  a  band 
dcLir^f  limit  of  resolution  da 
about  thi  interesting  frequencies 
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Fig.  3.  Interferogram  Produced  by  an  Incadescent 
Lamp,  the  Detector  Being  a  Lead  Sulphide  Cell. 


ffi  ’  *•'»  an’  This  analysis  can  be  made,  either  during  the 
measurement  or  afterwards  from  a  magnetic  tape  upon  which  the 
interferogram  is  recorded.  0m  can  also  record  the  interfero* 
gram  with  a  recording  pen  ami  numerically  calculate  the  F.T. 
of  the  function  thus  obtained. 

2.  THEORETICAL  'PPARATUS- FUNCTION 

Let  us  consider  a  Michclson  interferometer  illuminated  by 
a  point  source  placed  in  the  focus  of  the  entrance  collimator. 
One  of  the  mirrors  is  strictly  parallel  to  the  image  of  the 
other  formed  by  the  beam  splitter. 
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Fig.  4.  Above  ,  Inti  r» c rograra  of  a  Radiation  Band 
About  1.0  fi.  Below,  th.  jk  i  trum  computed  from 
the  interf  croaran:  cl  --  i;  ua,  t  tequcncy  in 

cm’*;  the  •  rninav.  .  inv-  nsity. 


If  B(a)  is  the  spectrum  of -tin  omitted  light,  the  intcr- 
ferogram  has  the  form. 

p  *  p  +" 

(II. 1)  1(4)  »  \  B(a)cos2r<roJflr  «.  \  B  (cr)oos2vtr4d0, 

»'.«•  P 
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Bp(e)  being  the  even  part  of  the  function  B(o)  defined  by 
Bp(<r)  *  j  [8(0  +  M-Ol  (Pig.  S).  .  “ 


AA 

K  &V  •«,  -ft 


_ZA  AV. 

-•*  -*>  if  * %  »' 

Fig.  5.  Even  Part  of  an  Optical 
Spectrum. 


If  one  could  know  the  interf erogram  completely  for  4  - 
varying  between  0  and  *°  (hence  in  fact  between  -•  and  -H® 
since  it  is  symmetric),  the  spectrum  could  be  reproduced  ex¬ 
actly  by  the  operation 

.  "  -V 

n  -H® 

(II.  2)  B  (or)  «  \  1(d)  cos  2*o6  dd. 

P 


But  in  fact,  during  the  measurement,  6  varies  only  between  0 
and  a  maximum  value  L;  the  value  B^(ai)  computed  for  a  parti¬ 
cular  wave  number  a j  will  be  only  an  approximate  value  of 
Bp(a4),  which  one  expresses,  by  saying  the  resolution  is  finite 
The  expression  (II. 2)  can  be  put  in  the  form 


(II. 3) 


p  -H® 

B’(c  )  =  \  1(d)  D(S)  cos  2TOd  dd, 

P  1  J-to 


D(d)  being  a  square  pul-e  function  such  that 


D(6)  «lif-L<d<+L, 

D(d)  *  0  for  d  <  -  L  and  d  >  +  L* 
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/ 


(1)  Form  of  the  Apparatus-Functicn 


The  cosine  transform  of  the  product  of  two  even  functions 
is  the  convolution  of  the  cosine  transforms  of  each  of  the  two 
functions,  whence 


UI.4)  B •(„,)  -  TcosI«B»  *  Tcos[D(«)1. 

The  cosine  transform  of  1(d)  is  precisely  the  observed  spec¬ 
trum  Bp(a).  Calling  f(<r)  the  cosine  transform  of  0(d),  (II. 4) 
can  be  written 

(II. 5)  BJOTj)  *  J  Bp(a)  [t(af(r)  +  f(*a«*)l  d*  4 

*  C  B  (a)  f(*.-<0  da. 

J  «  P  1 


The  spectrum  obtained  is  the  convolution  of  the  even  part  of 
the  initial  spectrum  with  the  function  f(a)  which  we  shall  call 
the  apparatus-function.  The  apparatus-function  thus  defined 
has  a  dimension,  that  of  length  in  the  case  where  the  spectrum 
B(a)  is  calculated,  that  of  time  when  B(iO  is  calculated. 

When  the  spedtrum  is  reduced  to  a  single  line  of  wave 
number  aQ,  of  luminance  unity,  the  reproduced  spectrum  hat  the 
form 

(11.5b)  B^(o)  -  l  [f(cr-a0)  +  f(<r4*r#)J. 

If  the  F.T.  i&  tiken  by  a  method  permitting  the  calculat¬ 
ing  of  luminance  densities  corresponding  to  negative  wave  num¬ 
bers  (which  have  no  physical  reality),  the  apparatus-function 
contains  two  maxima  centered  about  ff0  and  -  and  of  height 
f/2  (Fig.  6). 

When  the  F.T.  of  a  constant,  say  one, is  taken  the  only 
maximum  is  found  to  be  centered  about  the  zero  frequency  emu 
of  height  double  f # ,  the  two  preceding  maxima  baing  superposed. 
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One  can  define  a  normalized  apparatus-f unction  which  will  be 
the  ratio 

(11. 6)  r(o-a0)  =  Bp(ff)/Bp(cr0) 

b?tween  the  density  of  luminance  calculated  for  the  wave 

number  a  near  a  to  that  calculated  for  the  number  c  cor- 

o  t 

responding  to  the  line  under  consideration.  The  apparatus- • 
function  can  be  put  in  the  form  • 

(11. 7)  f(a)  «  f  #f  •  (»)♦ 

with  .}  ' 

(II. 7b)  f0  -  2qL  , 

q  being  a  numerical  factor  he.e  equal  to  the  mean  value  of 
the  function  0(6),  hence  one,  and  which,  more  generally,  will 
be  equal  to  the  mean  value  of  the  function  A(6)  defined 
later,  whence 

(11. 8)  f(o)  »  2Lf(ff)»  2L  S-fScir*  • 

This  function  f’(ff)  has  zeros  for  “  V^2L  (Fig.  7,  curve  a) 
it  lias  very  important  secondary  maxima  and  one  may  find  it 
desirable  to  apodir.c  it.  To  this  end  it  is  sufficient  to  re¬ 
place  the  square  pulse  function  D(6)  by  another  function  A(6) 
which  has  a  profile  better  adapted  to  the  problems  of  spectro 
scopy.  This  is  especially  simple  in  the  case  where  the  F.T. 
is  made  numerically.  There,  it  is  sufficient  to  replace  the 


*rf**r» 
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values  of  1(0)  by  tlioso  of 
the  product  1(6)  A(6);  this 
is  numerical  apodizing.  The 
apodir.ing  problem  has  been 
tin  subject  of  numerous 
papers ,  in  particular  those 
ol  Mile.  b.  Dossier  [26] 
and  J.  Arcac  [27]. 

We  ^arc^  going  t  o  con  - 
.aider  two  particular  forms 
of  the  weighting  function 
A<6). 

1^.  A  (6)  is  a  triangle 
function  such  that 


Fig.  7.  Theoretical  Apparatus - 
Function. 


‘e:  > 


A^  (6  )  =  1  -  Jd/l !  for  -  L  <  6  <  +  L 


and  A^fd)  -  0 

Under  these  condition; 


f  or  J  ^  •  L  and  A  ^  +  L. 


(II.9)  q  -  ^  §  +L  (l  -  \i\)  a*  -  0.5, 

This  is  the  same  apparatus-function  as  that  of  the  grating 
spectrometer  (Fig.  7,  curve  b). 

If  one  chooses  [27]  (Fig.  >0 

A^fd)  *  [1  -  (6a/ L*  )]*  a  or  -  L  <  6  <  +  L, 

q  -  -^r  \  *  0.533, 

2L  \ 
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Fig.  8.  The  A  (6)  Func¬ 
tion.  " 


•(II. 10)  »(e)  •  1.066  L  f*(w), 

with 

j: ...  •  f’(<0  *  (2tla)~%/%  Jf/>(2»l4r). 

This  function  (Fig.  7,  curve  c)  has  its  first  zero  for  0  very 
close  to  l/h\  its  breadth  at  mid-height  is  practically  the 
same  as  that  of  theTf unction  ((sin  »L«r)/»Lir)* ,  but  the  ampli¬ 
tude  of  the  first  minimum  is  only  a  twenty- fifth  of  that  of 
the  first  maximum,  and  that  of  the  next  maximum  Is  only  a 
hundredth. 

(2)  Limit  of  Resolution 

The  Rayleigh  criterion  can  be  applied  only  in  the  case 
where  the  apparatus -function  has  zero  minima;  hence  in  general 
one  selects  a  criterion  which  gives  nearly  the  same  results  as 
t  in  the  ordinary  cases.  One  is  led  to  take  as  a  measure  of  the 
limit  of  resolution  6a  the  breadth  at  0.404  of  the  maxlnum  of 
the  apparatus -function.  Since  the  F.T.  of  a  function  gets 
narrower  as  the  function  gets  wider,  one  sees  immediately  that 
the  limit  of  resolution  varies  inversely  with  L,  the  exact  law 
being  determined  for  each  particular  case. 

The  study  of  limits  of  resolution  rests  on  the  study  of 
the  breadths  at  mid-height  of  the  functions  f '(r), 
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If  there  is  no  apodization,  6a  -  3/2L,  whence  a  theoreti¬ 
cal  resolution  Rfl  =  2Lff.  With  the  weighting  function  A l(6)  or 
A#(6)  one  has  very  nearly  6a  =  l/L ,  with  a  theoretical  resolu¬ 


tion  R 


L a 


( i ) 


0  —  .  Hence,  the  resolution  is  proportional  to  the 

maximum  phase  difference  attained.  This  result  is  completely 
analogous  to  that  of  the  classical  spectrographic  method. 

Figure  84b  shows  two  spectra  of  the  night  skytaken  from 
the  same  interf erogram.  For  the  first  at  resolution  35(1,  a 
F.T.  has  been  made  of  that  portion  of  the  interf erogram  for 
which  6  varies  between  0  and  0.4  mm.  For  the  second  at  re*^ 
solution  1,000,  the  length  of  interf erogram  used  is  nearly 
triple.  '  '  >  *V.  i  ■'  V  ' 


(3)  Numerical  Fourier  Transforms 

Any  spectrum  studied  is  limited- -whether  by  the  emission 
itself,  by  the  transparency  of  the  materials  making  up  the 
interferometer,  or  by  the  sensibility  zone  of  the  detector. 

It  is  this  important  fact  that  allows  us  to  make  F.T.  from 
discrete  values  of  1(6}  and  to  use  modern  high-speed  calcula¬ 
tors. 

A.  FOURIER  TRANSFORM  AND  THE  THEORY  OF  SMFLINC.  Accord¬ 
ing  to  the  sampling  theorem,  all  information  concerning  a 
function  of  a  bounded  spectrum  is  contained  in  a  countable  in¬ 
finity  of  discrete  values  of  the  function,  and  interpolation 
formulas  can  be  found  which,  based  on  these  points  permit  the 
reconstruction  of  the  function.  These  questions  are  treated  in 
detail  by  Shannon  f28j.  Woodward  [29],  and  Kohlenberg  (30} . 

• 

And  so  all  the  bits  of  information. contained  in  the  intcr- 
ferogram  which  is  Known  co  have  a  bounded  spectrum  arc  given  by 


(,)In  this  case  one  can  say  further  that-  the  theoretical 
resolution  is  equal  to  the  number  of  recorded  fringes. 
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discrete  values  of  1(d).  As  a  consequence  it  is  sufficient 
in  calculating  the  spectrums  to  deal  with  a  particular  number 
of  equidistant  values  taken  from  the  int'erferogram.  The 
calculation  of  the  Integral 


(H.ll) 


B;<v  -  2 


£ 


+L 


1(d)  A(ft)  cos  2»o.«  dft 


is  replaced  by  that  of  the  sum 

(11.12)  "Bp(^1)s=Ml0A0+2IfcA/tcos2rfflfc  +  ...+2InfcAnfccos2»an»J. 


-v  V 


Inft  are  the  values  on  the  interferogram  for  the 
values  0,  ft,  . ..,  nft  of  the  phase  difference;  A#,  Afc>  ..., 
are  the  corresponding  values  of  the  apodizing  function  ^ 

The  sum  (11.12)  is  equivalent  to  the  integral  '7 

j1  1(d)  A(d)  R^ft)  cos  2w t6  dft , 


(11.13)  BJ^)  -  A 


Rft(d)  being  a  Dirac  distribution  of  periodic  sfepport  with 
interval  ft  (Fig.  9);  we  deduce  immediately 


f 


,  ft 


Lu^t 


Fig.  9.  Grating  function  K^(ft). 


/»  T*> 

(11.14)  BJO^)  -  i  J  Bp  (a)  lF(oJ-a)  +  P(.el.a)J  do. 

* 

F(a)  is  the  F.T.  of  the  product  of  the  two  functions  A(d)  and 
1:  Rft(d);  since  the  F.T.  of  a  periodic  Dirac  function  with 
interval  '■  a,  w.-.icpt  for  a  factor  l/h ,  another  periodic  Dirac 
distribution  with  interval  l/h  (Fig.  10),  a  new  apparatus- 
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I 


M 


i 


Fig.  10.  Grating  Function  .R^C*) . 

i  • 

function, can  be  found: 

F(<r)  «2t  F'(a), 

with 


(11.15)  F 1 (a) 


■  £  f 

M=.« 


s  in  2th  '  r-tf 

k 


2th 


(?-“)] 


in  case  A(d)  is  a  square  pulse  function  and  F(e)  «  1.066L  F’(e) 


where 


Nll.16)  F *  (ff )  *  S  £  2th  (f-a)]’8/*  Jt/t[2#4 


with  the  apodizing  function  Aa(d). 

The  new  apparatus-function  is  made  up  of  a  series  of 
identical  maxima  which  occur  at  every  multiple  of  l/k  (Fig.  11). 
It  can  be  said  that  the  even  part  of  the  spectrum  is  examined 
simultaneously  by  two  series  of  apparatus-functions  centered 
one  on  <r  ,  the  other  on  -  a ^  (Fig.  12).  The  examination  of 
the  spectrum  is  made  here  by  using  a  calculation  like  that  of 
(11.12)  for  a  new  valut  <*a  of  o’.  During  the  spectral  examina¬ 
tion  the  two  series  oi  apparatus -f unctions  have  translations 
in  opposite  directions  (Fig.  15). 
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B.  CHOICE  OF  INTERVAL  ft  .  This  is  a  very  important 
problem  because  the  computed  spectra  are  going  to  repeat  at 
multiples  of  l/fc  and  because,  for  the  result  to  be  useable, 
they  ought  not  be  entangled.  This  amounts  to  saying  that 
a  single  maximum  of  the  apparatus-function  ought  to  explore 
the  spectrum.  We  shall  call  <r M  and  <r  the  limiting  wave 
numbers  of  the  spectrum  and  A<r  *  m  °m  the  spectral  interval 
used  (fig.  14). 


Ac 


Fig.  14. 


a.  Relation  Between  ft  and  Ao.  Three  cases  can  dis¬ 

tinguished  according  to  the  position  of  A<r  with  respect  to 
the  zero  frequency. 

a)  <r  *  0.  a  varies  between  0  and  in  order 
tn  n 

that  there  be  no  overlapping  o.f  spectra,  it  is  sufficient  that 

l/ft  =  2  da  or  ft  *  l/(2Ae)  *  l/(2o-M)  (Fig.  15a).  This  is  the  1 

result  usually  stated  by  saying  that  one  ought  to  record  two 

points  per  period  of  the  signal  coi'responding  to  the  highest 

* 

freqv  ■■ncy  contained  in  the  spectrum. 

/J)  *  X^ff  (x  being  an  integer).  It  is  easily 

seen  in  this  case  also  that  the  orders  do  not  overlap  if 


l/h  «  20c  or  ft  «=  2/{2&tr)  (Fig.  15b).  The  distance  between 
two  points  taken  fro*  the  interferogra*  in  centimeters  is 
one  half  the  inverse  of  the  spectral  Interval  used,  measured 
in  cm"1. 

y)  Any  <rM  whatever.  “  X  +  X  being  an 

integer  and  f  a  fraction.  Then  U r  is  replaced  by  another 
quantity  Aa’  defined  by  &r'  *  °yt/X •  "*  is  larger  than  Ao  and 

has  the  same  upper  bound.  This  reduces  to  the  preceding  case 
with  the  choice  ft  »  (Fig.  15e). 
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Fig.  15.  Repetition  of  Numerically  Calculated 
Spectra. 


Numerical  Examples.  The  first  example  illustrates 
the  case  where  the  spectrum  extends  from  0  to  aM;  we  have 
computed  the  F.T.  of  the  function  sin*  *0^6/ (*0^6  j*  using  36 
discrete  values,  the  distance  between  any  cwo  successive  values 
being  l/(2o ^) .  The  calculated  spectrum  repeats  every  2ff^;  the 
resolution  is  *  18  because  we  have  used  only  a  short  length 
of  1(6);  also  the  apexes  of  the  triangle  are  slightly  rounded 


off  (Fig.  16). 


Fig.  16. 
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In  the  second  example  the  calculated  spectrum  is  that  of 
an  incadescent  lamp  whose  light  passes  through  filters  which 
isolate  a  band  of  radiation  in  the  neighborhood  of  1  ft. 

Knowing  the  transmission  curve  of  the  filters  we  can  de-  . 
duce  from  them  that 

dff  *  3,410  cm'1, 

<rM  =  12,850  cm'1  and 

•'  S'- 

ff^/acr  *  3  +  0.8.  j 


The  domain  do*  is  then  defined  such  that 


Pig.  17.  Spectral  Domain  Occupied,  3,410  cm'2; 

h  -  0.11695  x  io' 3  cm;  Number  of  Points  Taken 
From  the  Interferogram,  25. 


Figure  17  shows,  for  this  value  of  fc,  the  agreement  be¬ 
tween  the  position  of  the  observed  spectra  and  thoap  of  the 
computed  spectra. 


]  - 

NAVWEPS  REPORT  8099 


C.  RELATION  BETWEEN  THE  NUMBER  OF  POINTS  CHOSEN 
AND  THE  NUMBER  OF  SPECTRAL  ELEMENTS  STUDIED" 

Wc  have-  just  seen  that  the  interval  ft  depends  only  on 

the  breadth  and  the  position  occupied  by  the  spectral  interval** * 

In  addition,  the  limit  of  resolution  depends  essentially  on  the 

maximum -phase  difference  L  attained  and,  secondarily  on  the 

form  of  the'  apodizing  function  chosen.  Hence,  the  exploration 

of  a  given  spectral  domain  with  a  given  resolution  poses  two 

distinct  problems:  determination  of  L  and  choice  of  interval  ft. 

• 

The  total  number  of  points  to  choose  is  given,  within  a  unit,  by 
n  =  L/fc,  and  there  exists  a  relation  between  n  and  the  number  M 
of  spectral  elements  contained  in  the  spectrum  under  considera¬ 
tion.  The  spectral  element  is  that  portion  of  the  spectrum  in¬ 
cluded  between  two  points  that  can  be  considered  as  barely 
separated  by  the  instrument;  it  has  the  value  6<r  and  M  =  /iff/do. 

A  general  argument,  due  to  Shannon,  allows  one  to  predict  that 
starting  with  n  bits  of  information  contained  in  a  function, 

V^one  can  find  n  independent  values  of  its  F.T.  [31j.  We  are 
going  to  show  that  the  relation  between  n  and  M  varies  with  the 
spectrum  considered,  and  in  certain  particular  cases,  we  find 
n  *  M. 

Consider  at  first  the  case  where  vm  *  0  or  better  yet 
that  where  *  mCt/r t  m  being  an  integer.  Then  ft  »  l/(2d<r>  and 

n  -**  1 /ft  *  2L Mr. 

If  the  interf erogram  is  not  apodized,  *  l/2L  and 
n  *  aa/fia  *  M,  that  is  to  say  that  the  number  of  points  to  be 
chosen  is  equal  to  the  number  of  spectral  elements. 


(,,This  is  true  in  the  absence  of  noise.  In  practice, 
to  obtain  the  best  signal/noise  ratio,  one  decreases  ft  con¬ 
siderably,  the  number  of  points  chosen  being  multiplied  by 
a  number  of  the  order  of  5. 


< 
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If  apodizing  is  done  as  indicated  above,  then  obviously 
60  «  l/L  and  n  *=  2M’,  M1  being  the  new  number  of  spectral 
elements.  Its  value  is  half  that  of  M. 

In  the  most  unfavorable  case,  is  very  little  under 
(Fig.  18);  then  1  =  is  nearly  equal  to  2 ta  and 
q  *  Aff’/fia  is  very  close  to  2M1  or  4M1 
according  to  whether  or  not  apodizing 
has  been  used.  Hence,  in  conclusion. 


with  apodizing,  the  number  of  points 
to  choose  varies  between  2  and  4  'times' 


Fig.  18. 


the  number  of  spectral  elements. 


Remark.-- We  have  just  seen  that  for  a  given  resolution 
the  number  of  points  to  choose  is  greater  the  more  extended 
the  spectrum.  Although  the  advantage  of  this  method  over 
classical  methods  is  greater  the  higher  M  is,  one  can  consider 
u^ing  it  for  the  high  resolution  study  of  a  weak  spectral  do¬ 
main,  for  example  the  hyperfine  structure  of  a  line  in  the 
infrared.  If  M  is  of  the  order  of  100,  the  number  of  points 
to  take  will  be  rew,  in  the  order  of  hundreds. 


l>.  COMPUTING  METHOD 

We  constructed  a  program  for  computing  Fourier  transforms 
on  the  IBM  704 1 3  * .  In  order  to  achieve  a  reasonable  computing 
time,  the  problem  became  one  of  finding  a  recurrence  formula 
for  calculating  cos  2a,ajn«  based  on  cos  2*0^.  Indeed,  com¬ 
puting  cosines  by  using  either  series  or  tables  resulted  in 
prohibitive  calculation  times. 

We  used  the  Chebyshev  formula  (32]  • 

cos  (p  +  l)x  -  2  cos  x  cos  px  *  cos  (p  -  l)x. 


<>,The  program  is  now  in  the  IBM  general  library. 


I , 
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If  the  cosines  are  computed  to  six  significant  figures,  the 
error  between  the  value  computed  from  a  series  expansion  and 
that  obtained  after  12,000  iterations  is  It  thar  3/50,000. 

Program  Checking. --We  have  calculated  the  F.l.  of  a 
function  containing  24.5  sine- wave  arcs  of  unit  amplitude, 
and  two  points  per  period  have  been  taken,  the  interval  h 
being  taken  equal  to  unity.  The  maximum  height  of  the  calcu¬ 
lated  function  is  f#  in  this  case  and  not  \ Jl  [see  (II. 5. b) 
and  (II. 10)]  because  of  the  overlapping  of  the  orders  due  to 
the  par ticoTar^ choice  of  k .  -  ■  * 

2qL  «  1.066  x  49  *  52.2634.  The  number  calculated  by 
the  machine  is  52.25-82,  that  is  to  say  an  error  of  the  order 
of  1/10,000. 

The  calculation  is  made  in  two  parts: 

1" During  the  very  rapid  first  part,  in  each  address  the 
value  1^1  of  the  interferogram  at  the  point  6  -  pfc  is  replaced 
by  the  corresponding  value  1^  A^  of  the  weighted  interferon 
gram.  I 

2* .  Then  the  calculation  proper  of  tne  F.T.  begins.  For  each 
new  value  of  u,  the  first  cosine,  cos  2»irk ,  is  calculated  using  I 
series;  the  cosines  which  fellow  are  computed  by  a  re  cur  re  it  oe 
formula  and  the  sum 

B’(<0  *  k  [ IA  +21  A  cos2ir<rft+. .  .+21  A  cos2»enk  ]  , 

<  ii  n  n 

is  computed.  This  calculation  is  performed  in  I  •  0.6  x  I0~a  sec 
per  input  point,  per  output  point. 

If  n  points  are  taken  from  the  interferogram  and  k  spectral 
densities  calculated,  the  computing  time  will  be  t  *  0.8  nk  X 
10“*sec. 

Assuming  that  it  is  necessary  to  compute  the  intensity 
for  three  points  per  spectral  element  in  order  to  facilitate 
the  drawing  of  the  spectrum,  the  total  time  t  in  the  most 
favorable  case  is  t  ■  0.8  x  io~*  sec  x  2M  x  JM  -  0.48  Ma  10~a*ac 
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and  in  the  most  unfavorable  case*  t  =  M*  10”*  sec. 

It  is  interesting  to  remark  that  these  calculation  times 
correspond  to  the  case  where  the  T.T.  is  made  from  an  inter- 
ferogram  containing  no  noise.  In  practice,  to  obtain  the  maxi¬ 
mum  s/n  ratio,  it  is  necessary  to  pick  more  points  from  the 
interf erogram^  and  the  calculation  time  increase.  This  ques-_ 
tion  will  be  studied  in  detail  in  Chapter  V.  < 

Note . --It  suffices  to  calculate  one  point  per  spectral 
element  in  order  to  collect  every  bit  of  information  contained 
in  the  interf erogram,  or  M  points  in  all.  But  we  have  granted 
that  it  was  convenient  to  calculate  Q  points  per  element  (Q 
being  at  least  equal  to  3)  in  order  to  facilitate  drawing  the 
spectrum.  But,  the  spectrum  is  the  convolution  of  a  sequence 
of  Dirac  functions  with  the  apparatus-function.  The  Dirac 
functions  are  spaced  6<t  apart  and  have  heights  the  M  values 
which  were  calculated  previously  (theorem  of  Dirichlet). 

Since  it  is  possible  to  perform  convolutions  with  a 
numerical  calculator  (see  Chapter  IV  and  V)  and  since  the  cal¬ 
culation  is  relatively  rapid*  one  can  consider  calculating 
intermediary  points  in  making  the  convolution  of  the  M  Dirac 
functions  with  the  apparatus- function  (limited  to  the  domain 
where  the  latter  has  a  noh-fi^ylihible  value).  One  would  thus 
obtain  a  reduction  of  the  calculation  time  by  a  factor  of  Q. 

We  propose  to  write  this  as  a  computing  program  in  the  near 
future.  J 
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.  III.  REAL  APPARATUS- FUNCTION 

Under  conditions  attainable  in  practice,  the  spectrum- 
obtained  is  always  different  from  the  one  that  would  result 
from  the  convolution  of  the  spectrum  under  consideration  with 
the  theoretical  apparatus-function.  The  preceding  results, 
in  fact,  have  been  established  for  a  beam  of  zero  extent. 

They  are  modified  when  the  detector  receives  a  line  of  finite 
extent.  We  shall  study  these  modifications  in  the  case  where 
a  diaphragm  of  finite  dimensions  is  placed  in  the  focal  plane 
of  the  exit  lens  of  a  Michelson  interferometer. 

We  shall  show  that  there  results,  as  in  the  classical 
method,,  an  increase  in  the  apparatus-function,  with  a  corre¬ 
sponding  decrease  in  resolution.  The  study  of  the  product 
luminosity  x  resolution  allows  one  to  determine  the  best  condi¬ 
tions  for  using  the  instrument. 

In  general,  the  spectrum  obtained  when  the  source  is  re¬ 
duced  to  a  single  linp  of  negligible  width  is  different  from 
the  real  function  defined  above.  In  fact,  an  instrument  is 
never  perfect  and  the  precision  of  a  measurement  is  always 
limited.  We  shall  study  some  of  the  causes  of  perturbation 
of  the  apparatus-function,  such  as  irregularity  in  the  displaco 
ment  of  the  moving  element  and  the  maladjustment  of  the  inter¬ 
ferometer.  Finally  when  the  Fourier  transform  is  made  by  a 
method  sensitive  to  phase,  the  necessarily  limited  precision 
of  recording  the  zero  phase  difference  causes  ti  distortion  of 
the  apparatus-function. 

,  !  : 

1.  APPARATUS-FUNCTION  FOR  A  FINITE  BEAM.  WIDTH 

Several  authors  have  already  been  concerned  with  the  form 
of  the  diaphragm  that  has  to  be  placed  at  the  exit  of  a  two- 
beam  interferometer  {33 ]  and  with  the  dimensions  it  must  be 


I 
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given  in  order  to  have  a  suitable  depth  of  modulation  [34  ,35) 
or  to  measure  fractional  excesses  under  the  best  conditions 
| ll ] .  We  are  going  to  consider  this  question  with  the  view 
-of  using  the  Michelson  interferometer  as  a  spectrometer  by 
means  of  Fourier  transforms.  . 

The  effect  of  the  use  of  slits  of  finite  width  upon  the 
luminosity  and  resolution  of  a  grating  spectrometer  or  a  Fabry- 
Pcrot  spectrometer  has  been  treated  by  P.  Jacquinot  and 
C.  Dufour  [36]  and  P.  Jacquinot  [37];  the  same  study  for  the 
sisam  has  been  made  by  P.  Connes  [4)«  These  are  the  chief 
results  of  the  studies:  a  spectrometer  with  slits,  grating 
or  prism  is  an  apparatus  which  can  operate  a t^  very  different 
resolutions,  possibly  very  much  below  the  theoretical  resolu¬ 
tion  R  .  When  R  «  RA  the  product  luminosity  x  resolution  is 
o  o 

practically  constant,  and  simply  by  enlarging  the  slits,  one 
can  increase  the  luminosity  but  lose  in  resolution.  A  state¬ 
ment  of  the  results  is  made  easier  by  the  introduction  of  the 
reduced  variables: 

R/R  ratio  of  the  real  resolution  to  the  theoretical 

'  o 

resolution; 

Z  reduced  diaphragm  breadth  or  ratio  of  the  diaphragm 
"width  to  the  width  at  a  height  0.404  times  the  maximum  of  the 
theoretical  apparatus-function; 

•£  ratio  of  the  real  luminosity  to  the  maximum  luminosity 
r 

that  can  be  obtained  by  increasing  the  opening  of  the  diaphragm. 
In  the  case  of  slit  spectrometers,  has  a  simple  interpreta¬ 
tion.  It  is  the  ratio  of  the  real  luminosity  to  tnat  which  can 
be  obtained  in  the  absence  of  diffraction  with  slits  correspond¬ 
ing  to  l  -  1. 

Figures  10  and  20  show  for  the  Fabry-Perot  and  the  sisam 
the  variations  of  ^Ro»  **r  R/Ro  function.,  of  Z.  This 

last  product  passes  Through  a  maximum  for  a  normalised  diaphragm 
width  near  1  and,  for  best  results,  these  instruments  should  be 
used  near  their  theoretical  resolution. 
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Fig.  19.  Fabry- Perot  Spectrometer. 


Y75 

// 
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Fig.  20.  Sisan. 


The  use  of  the  Michel- 
son  interferometer  in 
spectrometry,  like  that 
of  the  Fabry-Perot  etalon, 
requires  a  circular  dia¬ 
phragm  (entrance  or  exit), 
centered  about  the  normal 
to  the  mirrors  and  li'mit- 
_ ing  a  bundle  of  rays  with 
solid  angle  Q  which  is 
slightly  inclined  to  the 
normal.  Under  these  con¬ 
ditions,  the  recorded  in- 
terferogram  *(d)  i.s  dif¬ 
ferent  from  the  theoretical 
interferogram  1(A)  and  its 
F.T.  is  different  from  the 
theoretical  apparatus-func¬ 
tion. 


Computation  of  the 
fee.oraed  Flux  and  of 
the  Real  Apparatus- 
function. 

The  calculation  of 


the  transmissed  flux  4>(6)  when  the  spectrum  studied  is  reduced 
to  a  single  line  of  negligible  width,  wave  nuuber  and  lu¬ 
minance  B  is  classical  [22,24,28]: 

sin(ifi«x  0/2i) 

(iii.i)  ♦($)  -  bu  —rirwzr—  cos  2*VU  -  <o/4*> h 


90  ■  U  is  the  beam  width.  It  has  already  been  shown  [38]  that 
using  a  finivo  width  has  two  consequences}  widening  the  appara¬ 
tus  function,  with  consequent  lowering  of  resolution,  and  , 
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displacement  of  the  spectrum  towards  the  lower  frequencies. 

The  correction  Aa  to  apply  to  the  calculated  spectrum  is 
&r  *  ofl/4  w. 

Numerical  Example .--As  we  shall  see  in  Chapter  VI,  the 
solid  angle  useable  for  studying  at  resolution  1,000,  a  spec¬ 
trum  about  1.5  n  is  Q  *  ra*  =  i  x  16  x  10"4  steradians.  Under 

these  conditions  /\s  =  6 . 500  x  4  x  10”4  *  2.6  cm”1. 

( 2 )  Determination  of  the  Bust  Conditions  for  Using  the 
Instrument 

Rigorously,  the  width  of  the  diaphragm  to  be  used  depends 
upon  the  weighting  function  A(6)  chosen.  We  consider  two  types 
of  apodization. 

a.  Weighting  by  a  Triangle  Function  A(6)  =  1  -  (6/L). 

The  variations  of  £r,  R/R0  and  R/Rfl  are  exactly  the  same 
as  for  a  sisam  used  with  a  diamond- shaped  diaphragm.  In  both 
cases  the  apparatus-function  is  like  that  of  a  classical  grat¬ 
ing  spectrometer  with  the  sisam,  it  is  necessary  to  choose  a 

diaphragm  such  that  l  «  1.  Then  R  =  0.80  R#  and  ^  *  0.S2 

(Fig.  20).  I 

b.  Weighting  the  Interf program  by  A(6)  -  fl  -  (^Vl*)!*. 
The  variation  of  R/R0  is  then  nearly  the  same  as  in  the  pre¬ 
ceding  case  because  the  two  apparatus-functions  have  the  same 
width  at  mid-height;  the  variation  of  is  slightly  different. 
The  product  R/Rq  still  passes  through  a  maximum  for  l  very 
near  to  1.  Under  these  conditions,  the  effective  resolution 
again  is  R#  =0.82  R#  and  *  0.87  (Fig.  21). 

The  angular  field  of  the  diaphragm  is  given  by  the  rela¬ 
tion 

dir.2)  a  =  /27r^  a  0.9  /57R. 

Let  us  suppose  that  for  the  maximal  phase  difference  L  attained, 
the  rings  in  the  exit  plane  of  the  objective  have  bright  centers 
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Fig.  21.  Variations  of  R/R^, 
and  R/R0  in  the  Michilson  Case. 


the  diameter  of  the  diaphragm  to  be  used  is  very  close  to 
9/10  of  that  of  the  first  bright  ring. 

(3)  Correction  to  the  f .ilculatcd  Frequencies 

The  wave  number  a c  or  the  radiation  for  which  the  spec¬ 
tral  density,  B,  is  calculated,  and  the  interval  k  of  the 
grating  function  occur  as  a  product  in  the.  calculation  of 
cosine.  The  wave  number  c  of  th*  radiation  having  the  same 
density  B  in  the  spectrum  studied  is  related  to  »c  by 

(111. 3)  ffc  =  a(l  -  (U/4*)|. 

In  taking  as  the  value  of  h  ,  the  interval  of  the  grating 
function,  a  number  such  that  <rc  J:  -  a  h  ,  it  is  possible, 
when  making  the  F.T.  numerically,  to  correct  the  frequencies 
automatically. 

One  is  then  led  to  choose 

(111.4)  ht  -  h  0J0  -  Ml  -  (U/4») ) . 
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We  shall  see  a  little  later  (§  111.2.(4))  that  to  eliminate 
the  lead  errors  of  a  moving  mirror,  it  i-s  necessary  to  adjust 
"the  phase  differences  with  a  reference  line  recorded  at  the 
same  time  as  the  interf erogram.  In  this  case  the  interval  h 
of  the  grating  function  R;?  (6)  [see  (11.13)]  is  measured  by 
taking  as  unity  the  wave  length  of  this  reference  line. 

Now  the  reference  beam  has  a  finite  width  generally  dif¬ 
ferent  from  that  of  the  beam  carrying  the  signal  to  bo  studied. 
Let  us  consider  the  case  where  the  reference  beam  follows  the 
same  path  in  the  interferometer  as  the  light  to  be  analyzed 
and  let  us  call  Q’  the  solid  angle  subtended  by  the  diaphragm 
in  front  of  the  reference  source  seen  from  the  center  of  the 
exit  objective.  The  reference  interf erogram  is  a  sine  wave 
corresponding  to  a  signal  of  wave  number  a .  -  a(l-(n’/4ir)]. 

From  this  it  can  be  deduced  that  if  h  is  the  value  of  the 
period  measured  in  the  wave  lengths  of  the  reference  line, 
its  exact  value  h  ’  is  given  by 

fc’  =  Mi  -  (Q’/4*)J. 

,  To  correct  automatically  and  simultaneously  the  fact  that  the 
beam  to  be  studied  and  the  reftrence  beam  have  finite  widths, 
it  is  necessary  to  take  as  measure  ^  of  the  period  which  occurs 
in  the  numerical  calculating, 

(III.  5)  -  Ml  -  (0/40  ]/[l  -  (C’/4»)] 

where,  as  Q  and  Q’  arc  always  small  compared  with  one 

(III. 6)  ht  =  i/l  +  —71). 

Remark .  --P .  Cortnes  [5‘t]  showed' 'that  the  phase  difference 
variation  of  interiering  rays  with  angle  of  incidence,  which 
necessitates  limiting  the  width  of  the  beams,  is  useful  when 
the  exploration  of  th«_  spectrum  is  made  in  space  (prism  and 
grating  spectrograph,  and  spectrometers,  Fabry-Perot  used  in 
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connection  with  photographing  the  rings),  but  is  detrimental 
when  it  is  made  in  time  (F.P.  or  Michelson  photoelectric  inter¬ 
ferometer).  If  one  inserts  $n  afocal  system  in  each  of  the  two 
beams  of  a  Michelson  interf erometer ,  it  is  possible  to  combine 
the  relative  displacements  of  the  mirrors  and  of  the  afocal 
systems  in  such  a  way  that  the  apparatus  no  longer  produces 
fringes.  The  interference  state  is  constant  in  the  entire 
space.  Luminosity  increase  is  practically  limited  only  by  the 
aberrations,  and  the  exit  aperture  if  it  cannot  be  completely 
omitted  will  be  considerably  enlarged;  all  the  previous  con¬ 
siderations  about  the  variation  of  the  product  •4  */*„  are 
obviously  no  longer  valid  in  the  case  where  afocal  systems  are 
used. 

2.  CHANGES  IN  THE  APPARATUS-FUNCTION  CAUSED  BY  DIVERSE  ! 

ERRORS  “  ' 

(1)  Effects  of  Imperfections  in  the  Optical  Components 

As  with  the  sisam,  the  aberrations  of  the  entrance  and 
exit  objectives ,  if  they  remain  small  with  respect  to  the 
angular  diameter  of  the  diaphragm,  have  practically.no  effect 
on  the  theoretical  resolving  power.  On  the  other  hand,  laminar 
-defects  in  the  interferometer  (mirrors  or  trihedrals,  separating 
and  compensating  plates  of  a  Michelson  interferometer,  lamellar 
gratings  of  Strong’s  apparatus  l 13])  bring  about  a  lowering  of 
the  modulation  which  varies  with  wave  length,  whence  a  loss  of 
luminosity  in  the  reconstructed  spectrum  which  increases  with 
decreasing  wave  length. 

With  grating  or  Fabry-Pcrot  [3]  spectrometers,  the  apparatus 
function  recorded  by  imperfect  optical  parts  is  the  convolution 
of  the  theoretical  apparatus -function  with  a  function  represent¬ 
ing  the  surface  faults.  For  example,  with  surface  having  a 
spherical  or  cylindrical  curvature,  the  fault  function  is  a  step 
function  which  increases  as  the  ratio  of  the  sagitta  to  the 
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wavelength  under  consideration  increases.  This  is  a  Gaussian 
function  in  the  case  of  poorly  inicropolished  surfaces.  These 
faults  are  accompanied  by  a  lowering  of  the  resolution  and  a 
loss  of  luminosity. 

In  the  F.T.  method,  it  is  the  recorded  interf erogram 
which  is  the  convolution  of  the  theoretical  interf erogram  with 
the  laminar  fault  function  D(6).  The  calculated  spectrum 
finally  is  the  product  of  the  theoretical  spectrum  by  the  F.T. 
d(5)  of  the  fault  function— Dfd ) ,  i.e.,  in  general, the  spectrum 
will  be  attenuated  with  respect  to  high  frequencies.  There  is 
no  widening  of  the  apparatus-function  nor  lowering  in  resolu¬ 
tion. 

Example . - - Consider ,  for  example,  the  case  of  a  circular 

mirror  having  a  spherical  curvature,  the  sagitta  being  Xg/2; 

that  of  the  wave  surface  will  be  X  .  which  means  that  instead 

.  o ' 

of  having  at  a  given  instant  a  well-determined  phase  difference 
6,  one  has  all  values  between  6  and  6  +  A  .  The  fault  function 
is  a  square  pulse  function  o(  width  Xfl  (Fig.  22b).  The  recorded 
interf erogram  I 1 (6 )  (Fig.  22c)  is  the  convolution  of  1(6)  (Fig. 
22a)  with  D(6)  and  the  spectrum  obtained  will  be  the  product 
(Fig.  22f)  of  the  theoretical  spectrum  B(a)  (Fig.  22d)  with  the 
function  (sin  (xo/o#)  (Fig.  22c),  the  F.T.  of  this  D(6) 

function  which  has  zeros  for  a  =  k  a#,  calling  o#  »  l/A#. 


Fig.  22.  Effect  of  a  Spherical  Curvature 
of  the  Surfaces. 
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In  practice,  it  will  obviously  be  necessary  that  the 
magnitude  of  the  flaws  be  small  with  respect  to  the  smallest 
wave  length  contained  in  trhc  spectrum;  this  amounts  to  say¬ 
ing  that  the  maximal  frequency  of  the  spectrum  will  be 
small  with  respect  to  aQ  (Fig.  22e);  under  these  conditions, 
in  the  whole "spectrum,  the  attenuation  will  be  weak  and  vary 
little . 

( 2 )  Effect  of  a  Maladjustment  of  the  Interferometer 

Wc  distinguish  the  two  cases  in  which  the  maladjustment 
due  simply  to  an  initial  poor  adjustment  of  the  optical  compo¬ 
nents  remains  constant  or  varies  with  the , displacement  of  the 
mirror * 

A.  CONSTANT  MALADJUSTMENT .  The  effect  of  a  constant 
maladjustment  is  equivalent  to  a  surface  fault.  If  the  inter¬ 
ferometer  mirrors  arc  stopped  down  by  rectangles,  and  if  one 
of  the  sides  of  length  L  is  perpendicular  to  the  line  of  inter- 
^  section  of  the  mirror  planes  which  make  an  angle  *•'  (Fp.gr .  23), 
then  the  two  reflected  wave  front-  make  an  angle  2  c,iand,  at  a 
given  instant,  the  phase-  difference  varies  between  6  Jand 


Fig.  2S.  Maladjustment  of  the 
Interferometer. 


{4,Thc  effect  of  a  maladjustment  upon  the  photoelectric 
measure  of  a  displacement  has  been  studied  by  G.  W.  Stroke  [40 J. 
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6  +  2c  L  (This  is  the  simplest  case).  Here,  by  putting 
A#  -  2cL,  one  sees  that  the  effects  on  the  interferogram 
and  the  spectrum  are  identical  to  those  of  the  sagitta  A0; 
the  calculated  spectrum  is  further  attenuated  by  a  factor 
(sin 

B.  VARIABLE  MALADJUSTMENT 

If  the  movement  of  the”  mirror  is  not  a  translation  (as 
the  result  of  an  imperfection  in  the  ways,  for  example),  there 
is  a  variable  maladjustment  in  the  interferometer.  Now  one 
must  write  that  e,  hence” A  and  ,  are  functions  of  6;  the 
signal-,  corresponding  to  each  frequency  a  in  the  interferogram, 
is  attenuated  by  a  factor  depending  on  6  and  ar,  P(<r,d)  = 

(sin  to/9^6) ] lircr/o  (6) J ;  the  new  apparatus-function  will  be 
the  F.T.  of  P(d).  In  contrast  to  the  preceding  case,  not  only 
its  height  but  also  its  shape  can  be  modified,  the  change  al¬ 
ways  being  the  greatest  for  the  highest  frequencies.  A  more 
detailed  inquiry  into  this  case  appears  to  be  of  little  use 
since  there  exists  a  method,  due  to  E.  R.  Peck  [41],  which 
removes  nearly  completely  this'  cause  of  error:  the  use  of  re¬ 
flecting  trihedrals,  unaffected  by  a  rotation  about  thair 
vertex . 

(3)  Effect  of  an  Error  Upon  the  Determination  cf  the 

Point  on  the  Interferogram  Corresponding  to  t 

Phase  Difference  6  ■  0. 

We  are  going  to  show  that  this  point,  called  in  what 
follows  the  "zero  point"  of  the  interferogram,  ought  to  be 
located  with  great  precision  (small  error  with  respect  to  the 
wave  length  studied)'8'.  Even  a  very  small  error  strongly 


(,)This  necessity  disappears  in  the  case  where  the  F.T. 
is  made  by  a  method  insensible  to  the  Phase  of  the  signals  con 
tained  in  the  interferogram:  the  use  of  an  harmonic  analyzer 

(Cont.  on  following  page) 
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perturbs  the  form  of  the  apparatus-function,  which  becomes 
asymmetric  and  has  important  negative  minima.  Moreover  its 
shape  varies  rapidly  with  a  so  that  in  the  case  of  a  broad 
spectrum  the  results  soon  become  hard  to  interpret. 

Suppose  that  we  have  chosen  for  the  zero  point  a  point 
which  corresponds  to  a  phase  difference  e.  According  to 
whether  e  is  positive  or  negative,  one  makes  the  F.T. ,  not  of 
1(d)  (Fig.  24a),  but  of  one  of  the  two  even  functions  repre¬ 
sented  by  the  Figures  24b  and  24c. 


Fig.  24.  Effect  of  an  Error  in  the  Choice 
of  Zero- Point  on  the  Interf erogram: 
a.  e  *  0;  b.  c  >  0;  c.  <  <  0. 


(»,cont)jj^e  conceived  by  P.  Jacquinot  (lj  and  made  by 

L.  W.  Hertz  (48 ]  or  of  a  numerical  transformation  of  the  type 
studied  in  §  11.2(3). 
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A.  CALCULATION  OP  THE  APPARATUS-FUNCTION 

a.  e  >  0.  In  order  to  calculate  the  spectral  density 
corresponding  to  the  wave  number  a  ,  in  lieu  of  using  the 
equality 

(111. 7)  B(uj  *  2  l  1(6)  A(6)  cos  2*<r  6  d6, 

*  .  1 

the  following  operation  is  performed 

(111. 8)  B »(<r)  »  2  (  I(6+e)  A(6 )  cos  2»a  6  d6 . 

1  Jo  1 

To  find  the  form  of  the  new  apparatus-function,  it  is  suf¬ 
ficient  to  consider  the  case  where  the  spectrum  is  reduced 

to  a  radiation  of  wave  number  a  of  luminance  B  “  1.  Then 

o 

(111. 9)  I(6+c)  =  cos  2ra#(6+0, 
so  that  (III. 8)  can  be  put  in  the  form 


B '  (a  )  «  2[ cos  2*a  c  \  Md)  cos  2 6  cos  2W  4  dd  - 

1  L  #  •'o  0  1 

(III.10) 


-  s 


in  2j §  A(6)  sin  2i tr(6  cos  2*,aj6  ddj. 


We  find  the  integral 

f\  CO 

2  ^  A(6)  cos  2rra06  cos  2«,aJ6  d6 


(III. 11) 


*  \  wvv + >i, 


which  gives  the  two  usual  maxima  of  the  apparatus-function, 
the  cosine  transform  of  A(6)  centered  on  »#  and  -  (Fig. 
25a). 

On  the  other  hand  (Fig.  25b) , 
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Fin.  2\.  Appji  atu  -1‘unct ionti  When 
T) n  i  c  is  an  Ertor  in  the  Determina¬ 
tion  of  6  -  0. 


(ina2) 


•  to 

2  ^  At«i)  sin  2»atf6  cur,  2»aj6  d6 


i  [/-  ^  * 


fc(tr)  lujini  an  odd  fund  it  n,  tin  sin*  transform  of  A(6),  that 
wo  shall  call  the  apparatu.,- 1  unct  ion  in  sines: 

I*  ~ 

(III.  13)  *(<r)  ~  2  ^  A(«i )  in  2 na&  dd  . 


Fig.  2t> .  inc  Transform  of  a  Square- 
Pulse  Function. 
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\ 


It  is  the  convolution  of  k(a)t  the  sine  transform  of  a  square- 
pulse  function  (Fig.  26)  with  a(a),  the  cosine  transform  of 
A  (6). 

The  spectral  density  calculated  for  wave  number  a x  is 

-then 

B’Co'j)  *  j  [cos  2r aQe  *((rima0)  +  sin  2* aflc  h(ai~ao)  + 


(III. 14) 


+  cos  2* <jq€  f(-ffi-o,0)  +  sin  2vaQe  n(-ai-ff())]. 


The  calculated  spectrum  is  thus  formed  by  the  two  peaks  of 
the  apparatus-function 

(III. 15)  F(<r)  =  cos  2mtqc  1(a)  +  sin  2raoe  h  (a) , 

centered  one  on  ,  the  other  On  — (Fig.  25c). 

The  expression  (III. 15)  shows  that  the  apparatus -function 
depends  on  the  frequency  analyzed  in  the  spectrum.  The  term 
ac  or  e/x  measures  the  error  made  in  the  choice  of  the  zero 
point  divided  by  the  wave  length  studied.  .  When  e/x  is  equal 
to  0  or  1,  F(cr)  *  f(a).  When  e/x  is  eqijal  to  l/4,  i.e,,  when 
an  error  of  a  quarter  of  a  fringe  has  bqen  made,  F(ff)  *  *(<r)# 
Between  these  two  extreme  cases,  all  intermediary  ones  are 
possible.  Figure  27  gives  the  function  F(a)  for  values  of 
e/x  taken  between  0  and  0.25.  The  asymmetry  of  the  function 
increases  with  the  displacement  by  translation  of  the  princi¬ 
pal  maximum.  In  Fig.  28,  the  variation  of  the  ratio  C  of  the 
height  of  the  negative  minimum  to  that  of  the  central  maxi¬ 
mum  for  e/x  varying  from  0  to  0.25  has  been  studied.  It  is 
interesting  to  note  the  rapidity  with  which  C  increases  for  . 
low  values  of  c/X  (Fig.  29). 


b,  c  <  0.  The  case  where  c  is  negative  is  easily  de¬ 
duced  from  the  preceding  case.  FCsr  a  given  wave  number  <r 
and  a  given  negative  value  t ’  of  the  error,  the .apparatus- 
function 
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Fig.  27.  Different  Forms  of  the  Apparatus-Func¬ 
tion  for  Different  Values  of  the  Phase  Displace¬ 
ment:  Curve  1,  c/A  *  0;  Curve  2;,  e/X  *  0.042; 

Curve  3,  c/A  =  0.084;  Curve  4*  «/X  *  0.126; 

Curve  5,  e/A  -  0.168;  Curve  6,  f/X  *  0.210; 

Curve  7,  c/X  *  0.250* 


F  ’  (a)  =  cos  2irae’  f(a)  -  sin  2rc r  je'J  h(a), 

is  identical  to  the  function  F(<r)  determined  for  a  positive 
value  of  c, 

F  (v)  *  cos  2worc  f(ff)  +  sin  2i rot  *(»), 
under  the  condition  that  loc’l  *  1  -  ec. 
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Fig.  28.  Variation  in  the  Fig.  2,9. 

Asymmetry  of  the  Apparatus- 
Function  With  Phase  Dis¬ 
placement. 


When  investigations  of  an  extended  spectral  domain  are 
made,  e/x  and,  consequently,  the  apparatus-function  vary  con¬ 
siderably  from  one  end  of  the  spectrum  to  the  other.  Figure  30 


Fig.  30.  Deformation  of  the  Apparatus-Function 
for  Phase  Displacements  Taken  Between  -  2»  and 
+  2». 
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shows  the  deformation  of  F(a)  for  e/X  varying  between  -  1 
and  +  1. 

E .  PRACTICAL  PROBLEMS  RELATED  TO  THE  DETERMINATION 
OF  THE  ZERO  PHAGE  DIFFERENCE.  When  the  spectrum  studied  is 
broad,  V  e  interf erogram  in  the  neighborhood  of  the  zero  phase 
difference  has  only  a  few  fringes,  the  intensity  of  which  de¬ 
creases  rapidly  and  (under  the  condition  that  the  interfero¬ 
meter  used  is  very  symmetric,  hence,  that  the  interferogram 
itself  is  symmetric)  the  zero  can  be  located  with  great  pre¬ 
cision.  In  general,  the  precision  is  limited  for  two  reasons: 

1* .  in  practice,  on  the  ideal  interferogram,  /there  are  always 
superposed  some  fluctuations  due  to  noise; 

2* .  sometimes  the  very  nature  of  the  interferogram  renders 
the  recording  difficult.  If  the  spectrum  considered  consists 
only  of  some  rather  close  lines  in  the  neighborhood  of  the 
zero  phase  difference,  the  modulation  is  practically  constant, 
and  one  is  uncertain  by  several  fringes. 

One  solution  to  these  two  problems  consists  in  recording 
simultaneously  the  interferogram  to  be  studied  (interferogram  1) 
and  that  given  by  a  source  with  an  intense  and  broad  spectrum  - 
(interferogram  2).  Practically,  one  uses  an  incandescent  lamp 
and  a  filter  isolating  a  broad  radiation  band  in  the  same 
domain  as  the  spectrum  to  be  studied. 

When  one  studies  an  emission  spectrum,  there  is  no  in¬ 
convenience  in  replacing  a  portion  of  a  fringe  or  an  entire 
fringe,  in  the  center  of  interferogram  1,  by  the  equivalent 
part  of  interferogram  2.  At  the  time  of  the  recording,  the 
passage  of  the  zero  phase  difference  is  easily  recognized.  It 
is  enough  then  to  replace  the  light  coming  from  the  incadescent 
lamp  by  the  spectrum  to  be  studied.  This  question  will  be  con¬ 
sidered  in  detail  in  Chapter  VI.  The  important  problem  of  the 
determination  of  the  zero  phase  difference  in  the  case  where 
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the  interf erogram  is  recorded  with  a  time  constant  which  makes 
it  asymmetric  will  be  considered  in  Chapter  V. 


C.  CALCULATION  OF  THE  SPECTRUM  EY  A  METHOD  NOT  REQUIRING 
A  KNOWLEDGE  OF  THE  ZERO.  Wc  are  going  to  show  that  it  is 
possible  correctly  to  reconstruct  the  spectrum  studied  without 
being  bothered  about  the  position  of  the  zero  point,  when 
the  interf erogram  is  recorded  for  6  varying  from  -  L  to  L, 
to  make  the  cosine  transform  B]  of  the  interf erogram;  to  make 
the  sine^tr^Hfrsftjrm  ;  and  to  take  as  the  measure  of  the 

spectrum(6)  B'  * x/b ’ 2  +B^ a .  Suppose  that  the  spectrum  contains 
only  a  single  line  of  luminance  B  =  1,  of  wave  number  o#;  the 
spectrum  obtained  in  this  case,  will  be  the  apparatus-function-* 
The  interf erogram  reduces  to  a  sine  wave  having  the  form 


7 16)  =  cos  2*oo(6+c)  =  cos  2»c rdg  cos  2*oe  -  sin  2*tr#6  sin  2*ec. 

It  is  composed  of  an  even  part  Ip(6)  *  cos  2rec  cos  2wr#4  and 
of  an  odd  part  =  sin  2'*ac  sin  2*o#6. 


B«(cr)  =  cos  2*ee  ^ 


«.• 0 


A(d)  cos  2*o  6  cos  2*od  dl 


■  i  cos  2*oc  If (o-o#)  +  f (o+o0 ) ) . 
n  +w 

B^(o)  «  -  sin  2*ce  j  A(6)  sin  2**#5  sin  2*el  dl 
-  -  j  sin  2*oe  lf(o-o#)  -  f(o+o#)J. 


The  calculated  spectrum  will  have  the  value 


‘#,This  process  is  used  by  J.  G.  Braithwaite  and 
J.  Brooks  (53). 


45 


NAVWEPS  REPORT  8099 


B’(a) 
(III. 16) 


*/[  b; <») j*+| »; <»>  j* 

-  j  Clf«r-<r,j1  +  ltOr+*c)|]. 


The  new  apparatus-function  |f(ff)|  has  the  same  width  as  the 
.usual  function  f (a).  Hence  the  limit  of  resolution  remains 
unchanged.  However,  this  method  has  several  inconveniences: 


1* .  it  requires  a  double  displacement  of  the  movable  mirror, 
which  can  become  difficult  for  high  resolutions; 

2* .  for  equal  time  durations,  the  signal/noise  ratio  is 
somewhat  less  advantageous  (this  question  will  be  studied  in 
greater  detail  in  the  following  chapter); 


3* .  the  calculation  time  in  the  case  of  numerical  analysis 
is  a  little  more  than  quadrupled  (there  are  twice  as  many  in¬ 
put  points  and  there  are  two  identical  calculations  for  the 
cosine  and  sine  transforms;  the  calculation  time  of  the 
squares  of  the  F.T.  and  of  the  square  root'  of  their  sum  is 
short  compared  with  that  of  the  F.T.  themselves). 


(4) 


i  jf 


Effect  of  an  Error  Upon  the  Measurement  of  the  Phase 
Difference 


The  reconstruction  of  the  spectrum  based  on  the  inter- 
ferogram  requires  the  knowledge  of  the  phase  difference  cor¬ 
responding  to  each  of  its  points.  This  is  best  accomplished 
by  simultt leously  recording  (with  a  two-pen  recording  device 
or  a  two-track  magnetic  recorder)  the  interf erogram  to  be 
studied  and  a  reference  signal  given  by  a  monochromatic  radia¬ 
tion  of  known  wave  length.  The  precision  with  which  6  is 
known  is  then  limited  only  by  the  noise  in  the  reference 
signal.  In  practice  this  last  consideration  does  not  matter 
in  problems  with  low  resolution,  the  only  ones  we  have  at¬ 
tempted  experimentally  up  to  now;  so  we  have  noi*  undertaken 
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this  study.  However,  laving  to  produce  and  record  this 
reference  signal  presents  a  certain  difficulty.  Another  much 
simpler  but  much  less  precise  procedure  consists  in  recording 
simultaneously  the  interf erogram  and  reference  time  pips  pro¬ 
duced  by  the  rotation  of  the  drive  screw**1.  Still  more 
simply,  one  can  drive  the  screw  and  the  recording  paper  by 
synchronous-  motors.  Wc  arc  going  to  study  the  deformations 
of  the  apparatirs-f unct ion  introduced  in  these  last  two  methods, 
by  mechanical  errors.  Our  conclusion  will  bo  that  they  can  be 
used  only' to  treat  problems  of  low  resolution,  for  example  in 
the  far  infrared. 

The  movable  mirror's  advance  is  never  rigorously  linear 
because  a  screw  is  never  perfect.  It  always  has  a  periodic 
error  of  period  0,  the  period  of  the  rotation  of  the  screw 
0  mf/V,  p  being  the  pitch  of  the  screw.  We  shall  confine 
ourselves  in  the  following  to  the  study  pf  the  "ghosts"  intro¬ 
duced  into  the  calculated  spectrum  by  an  error  Of  this  type t 
because  it  is  by  far  the  most  annoying. 

A.  CALCULATION  OF  THE  RELATIVE  HEIGHT  OF  GHOSTS.  In  a 
general  way,  d(t)  can  be  considered  as  the  sum  Of  two  terms 
(Fig.  31),  one  a  linear  function  of  time,  the  other  an  arbi¬ 
trary  function  of  time  with  mean  value  zero*  J  *  Vt  -  f(t). 

If  f(t)  is  a  periodic  function  of  period  0 ,  we  can  write  it 
in  the  form  of  a  Fourier  series 

f{t)  *  cos  2rnj  +  bn  sin  2»njj  ■  cos  (2#nj  - 

n-i  n»i 


(V,An  intermediary  precision  method  (sufficient  for  mid 
or  far  infrared)  uses  a  reference  signal  provided  by  the  moire 
fringes  obtained  with  two  gratings*  otic  of  which  is  attached 
to  the  moving  component. 


4  7 
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Fig.  31.  Variation  of  6 
as  a  Function  of  t. 


When  the  error  is  purely  sinusoidal 
of  period  S,  only  tha  first  term 
exists  and 

(111.17)  fi  =  Vt  +  cos 

Let  us  take  the  time  origin  at 
d  =  0;  the  expression  (111.17)  be¬ 
comes 

6  -  Vt  +  Cj  sin  (2*t^). 


Under  these  conditions  the  interf erogram  obtained  with  a  mono¬ 
chromatic  source  of  luminance  B,  of  frequency  (r^,  can  be  ex¬ 
pressed  as 

(III.  18)  I(t)  «  B  cos  2t0  (Vt  +  e  sin  (2»t/0)] 

■  B  cos  (Q#t  +  p  sin  wt), 

where  n#  «  2»<r0V  (mean  frequency)  and  p  *  2»ej/X0  (modulation 
index). 

One  is  thus  led  to  the  classical  problem  of  studying  the 
spectrum  of  a  frequency-modulated  sine  wave  ( [ 4 3 J  p.  407). 

This  problem  is  familiar  to  radio  engineers.  It  also  occurs 
in  the  study  of  ghosts  of  a  grating  which  has  been  drawn  with 
a  ruling  machine.  From  this  point  of  view  it  has  been  the 
subject  of  several  papers,  the  oldest  of -which  is  due  to 
Rowland  ([25]  p.  97,  [44],  [45]).  Wo  are  going  to  take  it 
up  in  our  particular  case  and  wo  do  not  limit  ourselves  to 
the  case  where  0  is  small  in  comparison  with  the  wave  length 
studied,  because  the  interferometer  screws  are  generally  of  a 
quality  very  inferior  to  that  of  the  screws  used  in  the  machines 
which  rule  gratings.  ./"**" 


4B 


NflVWi 


Equation  (III. 18)  can  be  put  in  the  fora 
(III. 19)  I(t)  *  B  cos  0#t  cos  (0  sin  wt) 

-  B  sin  not  sin  (.8  sin  wt). 

Using  the  classical  expansions  of  cos  (/J  sin  cot)  and 
sin  (0  sin  wt)  and  noticing  that  J_n(z)  -  (-l)n  Jn(s),  one 
gets 


(III. 20)  I(t)  * 


T  1 

B  J  (0)  cos  (0#  +  nw)  t. 

n“-* 


The  spectrum  of  the  signal  studied  is  composed  of  the 
principal  line  of  wave  number  and  of  an  infinity  of  ghosts 
located  on  both  sides  of  the  principal  ray,  of  intensity  Jn(0) 
and  of  wave  number  ±  n  0<r,  n  being  an  integer  and  J0  dlS 
Interval  between  any  two  of  these  ghosts*  . 

-  h*  * 

f 

being  the  pitch  of  the  screw. 

Since  the  Bessel  functions  become  negligible  when  the 
index  is  distinctly  greater  than  the  argument,  the  number  of 
ghosts  that  one  sees  is,  in  practice,  limited.  We  can  get  an 
idea  of  their  height  relative  to  the  principal  line  by  the 
following  comparisons. 

Suppose  that  the  screw  is  used  in  a  machine  ruling  grat¬ 
ings  and  that  the  gratings  so  obtained  are  used  in  a  Littrow 
mounting  at  grazing  incidence  (i.e.,  under  the  worst  conditions 
the  ghosts’  intensity  increasing  as  the  squar-e.  of  the  sine  of 
the  incident  angle).  One  can  show  that  the  relative  height 
of  the  ghosts  obtained  by  the  F.T.  of  the  interf erogram  is 
equal  to  the  ratio  of  the  amplitude  of  the  luminous  vibration 
difi<  <cted  in  the  ghost  to  that  diffracted  in  the  principal 
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line  whence  equal  to  the  square  root  of  their  relative  inten- 
It  is'  thus  much  greater. 

NUMERICAL  EXAMPLE.  Suppose  that  with  an  excellent 
machine  one  can  get  gratings  such  that  the  intensity  of  the 
first  ghost  far  an  incidence  of  60*  and  a  given  wave  length 
is  10‘>  times  the  principal  line's  wave  length.  Used  to 
record  an  interf erogram,  the  same'  machine  would  have  yielded 
for. the  same  wave  length  a  ghost  whose  relative  Intensity  is 
3.6  x  1(TV  -- 

Thet.o  facts  give  an  idea  of  how  much  more  difficult  it 
is  to  build  a  machine  sufficiently  precise  to  record  an 
interf erogram  and  to  take  a  F.T.  of  it  without  recourse  to  a 
reference  line  than  it  is  to  build  a  grating  ruling  machine. 

B.  EXPERIMENTAL  STUDY.  We  didn't  check  these  results 
directly  but  we  studied  our  interferometer's  screw  and  we 
deduced  the  relative  height  of  the  ghost  that  we  would  have 
gotten  in  taking  the  F.T.  of  an  interferogram  without  using 
the  reference  line. 

As  a  source  we  used  the  red  cadmium  line  and  we  recorded 
the  14,100  fringes  corresponding  to  a  variation  of  phase  dif- 
ference  of  1  cm.  We  measured  the  distances  between  every 
hundred  fringes  on  the  recording  with  an  estimated  error  of 
one-quarter  fringe;  the  displacement  of  the  moving  mirror 
was  adjusted  to  about  0.08  p .  It  was  then  possible  to  chart 
the  screw  errors  (Fig.  32).  A  periodic  error  was  clearly 
disclosed  with  period  equal  to  that  of  the  screw-  whose  thread 
./as  ;p  »  1  mm;  the  error  amplitude -wa3  about  c  *  0.S  p. 

The  ghosts  which  one  would  obtain  would  be  Off  »  x/2f>m 
S  cm-1  apart.  We  computed  the  relative  height  in  three  dif¬ 
ferent  cases. 

1* .  -  0.6438  p,  i.e.,  the  interferometer  is  to  be  used 

in  the  study  of  the  rad  cadmium  line.  The  modulation  index 
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P  is  4.88.  One  must  go  out  to  the  eleventh  ghost  to  find  an 
intensity  which  is  2/10,000  times  the.  theoretic  intensity  of 
the  principal  line  (Fig.  33a). 

2* .  A  ®1.5|*,  0*2.  The  line  studied  is  in  the  near 
infrared.  There  are  fewer  ghosts  (Fig.  33b). 

3* .  A#  «  15  jt,  p  *  0.2.  Only  three  ghosts  have  an  intensity 

greater  than  2/10,000  (Fig.  33c). 

I  • 

i 

Figure  33  shows  that  the>  spectral  interval  occupied  has 
width  100  cm"1,  60  cm"1,  or  jo  cm"1  according  as  the  line 
studied  is  in  the  visible,,  the  near  infrared,  or  the  mid~infrared 
part  of  the  spectrum.  If  the  limits  of  resolution  are  of  this 
order,  one  cannot,  in  practice,  recognize  the  ghosts. 

These  results  illustrate  very  well  how  necessary  it  is 
to  use  a  reference  line  when  it  is  desirable  to  exceed  a 
resolution  of  the  order  of  100,  in  the  visible  an J  near  infra* 
rad. 
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IV.  GENERAL  NOTfcS  ON  NOISE  IN  THE  FOURIER 
TRANSFORM  METHOD 

In  the  introduction  we  showed  by  a  simple  argument  the 
gain  in  the  signal/noise  ratio  that  could  be  expected  for 
equal  resolution  and  measurement  time  by  the  F.T.  method  in¬ 
stead  of  a  classical  method.  Recall  that  this  is  true  only 
in  the  infrared  region,  where  the  fluctuations  perturbing  the 
measurements  originate  in  the  detector  itself  and  not  in  the 
line  studied,  i.e.,  when  the  photon  noise  is  negligible  com¬ 
pared  with  detector  noise.  - 

The  noise  problem  of  the  F.T.  method  does  not  occur  in 
the  same  way  as  it  docs  in  the  classical  methods.  _In  fact, 
here,  the  result  of  the  experiment  is  the  interf erogram.  The 
interf erogram  noise  has  the  same  behavior  as  the  noise  in  the 
spectrum  obtained  by  a  classical  method,  but  the  noise  in  the 
“Spc^trr'um^bbtained  py  a  Fourier  transformation  is  the  F.T.  of 
the  noise  contained  in  the  interf erogram  and  has  a  special 
aspect. 

In  this  chapter  we  shall  study  the  behavior  of  the  noise 
in  the  spectra.  We  shall  show  that  the  s/n  ratio  in  the  spec¬ 
trum  varies  as  the  square  root  of  the  total  measurement  time  T 
and  how  one  can  predict  the  s/n  ratio  in  the  spectrum  using 
the  s/n  ratio  in  the  interf erogram. 

1.  CHOICE  OF  THE  BEST  CONDITIONS  FOR  RECORDING  THE 
INTERFEROGRAM ,  NOISE  IN  THE  INTERFEROGRAM 

The  procedures  for  recording  an  interf erogram  vary  accord¬ 
ing  to  whether  one  intends  to  make  the  F.T.  of  an  interf erogram 
with  a  harmonic  unelv.:c-r  or  by  a'  numerical  method.  In  the 
first  case  one  makes  a  continuous  recording  (for  example,  on 
magnetic  tape  [1,46]).  In  the  second  case  one  can  use  either 
a  pen  recorder  or  a  converting  system  which,  at  a. given  instant, 
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changes  a  voltage  into  a  number  written  on  punched  tape  *42], 
Whatever  method  is  used,  the  problems  relative  to  amplifica¬ 
tion  of  the  signal  furnished  by  the  detector  remain  the  same. 
Fellgett  was  the  first  to  insist  upon  the  necessity  of  intro¬ 
ducing.  only  perfectly  linear  elements  into  the  amplification- 
detection  system.  We  are  going  to  distinguish  several  possibi¬ 
lities  according  to  the  signal  frequency  to  be  recorded. 

1* .  If  the  detector  used  is  a  photomultiplier  and  if  the 
available  energy  is  low,  so  that  a  long  recording  time  is  re¬ 
quired  hence  a  low  velocity  V,  the  electrical  frequencies 
contained  in  the  signal  to  be  amplified  which  we  shall  now 
call  Fourier  frequencies  are  low  (of  the  order  of  one  cycle 
**  per  second).  One  can  then  use  a  direct  current  amplifier. 

3* .  Lead  sulfide  cells  have  a  non-uniform  noise  spectrum, 
.always  maximum  in  the  neighborhood  of  the  zero  frequency 
(scintillation  noise).  On -the  other  hand  their  response  to 
a  modulated  signal,  uniform  for  lower  frequencies,  decreases 
rapidly  for  higher  frequencies  (of  the  order  of  1,000  eps). 
Under  these  conditions  the  best  s/n  ratio  for  a  signal  occupy¬ 
ing  a  given  band  width  is  obtained  when  the  mean  signal 
frequency  is  of  the  order  of  100  to  1,000  cps. 

If  the  source  to  be  studied  is  Intense  enough  to  permit 
operating  at  such  speeds  that  the  frequencies  to  be  analyzed 
are  of  this  order,  an  alternating  current  amplifier  with  a 
sufficiently  large  band  width  can  be  used.  But  in  the  majority 
of  cases  the  usable  energy  is  low.  The  Fourier  frequencies 
will  be  very  low,  whence  the  necessity  of  chopping  the  lumi¬ 
nous  flux  with  a  modulator  (following  the  classical  method  in 
infrared  spectroscopy).  The  detection  of  the  electric  signal 
necessarily  will  have  to  be  a  synchronous  detection  116 J.  In 
fact,  we  shall  see  in  what  follows  that  with  the  F.T.  method 
*  one  can  detect  lines  which  give, in  the  interferogra%  sine 
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Fig.  34.  Schematic  Diagram  of  the  Recording 
of  an  Interf erogram.  Captions: 

Signal  lumineux  =  Light  signal 

Detecteur  *  Detector 

Ampli.  select if  =  Selective  amplifier 

Detection  synchrone  *  Synchronous  detector 

Filtre  passe-bas  -  Low-pass  filter 

Enregistr.  =  Recorder 

Signal  de  reference  *  Reference  signal 

Modulateur  *  Modulator 


waves  whose  amplitude  is  much  smaller  than  the  noise.  But  a 
classical  detection  method  (so-called  linear)  is  in  fact 
linear  only  if  the  s/n  ratio  is  much  greater  than  1  ( [47  J  p. 
359),  hence  the  absolute  necessity  of  synchronous  detection. 

The  detector-amplifier  system  to  use  in  this  case  is 
exactly  the  same  as  the  one  Which  serves  to  record  the  spec¬ 
trum  in  classical  methods  (Fig.  34).  The  electrical  signal 
before  synchronous  detection  is  amplified  in  a  selective 
amplifier  in  such  a  way  that  the  synchronizing  frequency  cor¬ 
responds  to  the  frequency  of  modulation  of  the  light  signal. 
The  band  pass  should  be  wide  enough  for  all  the  frequencies 
to  be  transmitted  without  phase  distortion  or  amplitude  at¬ 
tenuation.  After  synchronous  detection,  the  signal  passes 
through  a  low-pass  filter  before  being  recorded.  We  shall 
see  that  this  filter  plays  a  fundamental  r01e  when  the  F.T. 
is  made  numerically  using  discrete  values  taken  from  the 
interferogram.  When  it  is  made  by  a  harmonic  analyser  (which 
uses  the  whole  interferogram  and  not  just  discrete  values) 
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the  s/n  ratio  in  the  calculated  spectrum  does  not  depend  on 
the  time  constant  of  the  low- pass  filter.  The  rSle  of  the 
latter  is  simply  to  reduce  the  noi?e  which  must  not  saturate 
the  recording  device.  The  spectral  densities  A(i>)  and  0(v)> 
corresponding  to  the  frequency  v  of  a  signal  before  aiut  after 
filtering,  are  connected  by  the  classic  relation 

(IV. 1)  a(0  =  A(*)  |G(i> )  | *  , 


G(v)  being  the  complex  gain  of  the  filter. 

If  one  assumes  that  the  noise  spectrum  is  constant  in 
the  neighborhood  of  the  modulation  frequency  (white  noise) 
with  a  noise  density  ,  the  power  spectrum  of  the  noise  con¬ 
tained  in  the  interf  eronrahi  (or  in  the  spectrum  obtained  by  a- 
classical  method)  has  the  form 


(IV. 2)  »(iO  -  »#  |G(»Ol*. 

In  the  particular  case  of  a 
resistance -capacity  filter  with 
I  time  constant  RC  «  r  * 

(IV.  3)  ®<*0  -  ®/ll  +  (2tvt)*  J 

(Fig.  35) 

2.  APPEARANCE  OF  NOISE  IN  THE 
CALCULATED  SPECTRUM 

The  cosine  transform,  as 
described  in  Chapter  II,  is  a 
linear  transform.  The  recorded  interf erogram  Ir(t)  is  the  sum 
of  the  theoretical  interf erogram  I(t)  which  we  have  been  con¬ 
sidering  until  now  and  the  noise  .\(t)  which  can  be  regarded  as 
random  luminance  of  mean  value  jf-ro.  The  F.T.  of  I’(t)  is  the 
sum  of  the  spectrum  Bp(v),  F.T.  of  I(t),  and  noise  X(p)  which 
is  cosine  transform  of  x(t);  X(v)  is  a  quantity  of  the  same 
kind  as  »(»»).  It  is  a  density  of  the  random  luminance 


Fig.  35.  Noise  Spectrum  and 
its  Even  Part. 
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(luminance  per  cps),  with  mean  value  zero,  which  is  superim¬ 
posed  on  the  spectrum.  There  is  no  signal-noise (* *  inter¬ 
action  and  to  study  the  sum  signal  +  noise,  it  is  sufficient 
to  study  the  signal  separately  in  the  absence  of  noise,  which 
has  been  done  in  Chapter  II,  and  then  study  noise.  If  one 
represents  the  sum  signal  +  noise  £$,  i.e.,  the  result  of  the 
cosine  transform  of  I’(t)  as  a  function  of  the  signal  s  which 
would  l>e  the  result  of  the  cosine  transform  of  I(t),  one  gets 
Fig.  36  in  which  ay  represents  the  mean  square  error  of  the 

fluctuations  in  the  calculated  spectrum 
There  will  be  a  certain  definite  pro¬ 
bability  P,  related  to  the  shape  of  the 
statistical  distribution  of  noise,  that 
the  measured  signal  lies  between 
s  +  (Fig.  36,  curve  a)  and  s  - 
(Fig.  36,  curve  b). 

In  the  study  of  a  random  variable, 
two  problems  occur: 

1° .  the  determination  of  its  autocor¬ 
relation  function,  whence  the  calcula¬ 
tion  of  the  radius  of  correlation  from 
which  one  deduces  the  appearance  of  the 
curve  representing  the  random  variable; 

.  the  calculation  of  its  mean  power,  again  called  variance. 


of  the  Sum  Function 
Signal  +  as  a 

Function  of  the  Sig¬ 
nal  When  the  F.T.  is 
Made  by  Referencing 
the  Zero  Point. 


lB,We  shall  see  in  what  follows  that  when  the  spectrum 
is  obtained  as  the  square  root  of  the  sum  of  the  squares  of 
the  cosine  and  sine  transforms  of  the  interferogram  recorded 
from  -  L  to  +  L,  as  explained  in  paragraph  III.2(3)C,  the 
study  ot  the  signal/noisc  ratio  is  much  more  complicated. 


57 


NAVWEPS  REPORT  8099 


A  simple  argument  allows  cine  to  predict  the  appearance 
of  noise  in  a  calculated  spectrum  without  its  b^ing  necessary 
to  calculate  the  autocorrelation  function.  The  noise  x(t) 
recorded  on  the  interf erogram  between  instants  0  and  T,  which 

.  "I 

we  shall  call  oscillogram  noise  has  no  periodicity  whatsoever 
(Fig.  37a).  Its  radius  of  correlation  is  r.  However,  if  T  is 


Fig.  37.  (a)  Noise  in  the  Interf erogram; 

(b)  Signal  in  the  Interf erogram;  (c)  Noise 
+  Signal. 


eho„,.-r.  sufficiently  large  with  r.  pect  to  the  time  constant 
*  which  has  b<  on  used  to  ncs.u  *  (t) ,  then  practic  ally, one 
can  assume  (although  in  theory,  this  mav  not  be  entirely  e-act) 
that  a ( t )  is  a  periodic  design  of  period  T  and  can  make  a 
classic  Fourier  decomposition  of  it  ((47)  p.  507,  |4bJ).  In 
the  interval  (0,T),  x(t)  is  represented  by  the  Fourier  series 
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(IV. 4)  x (t )  *  ^  2T  cos(2wif 

n=-*° 


nt-v>n),  with  *n  =  $  , 


which  wc  can  put  in  the  fora 


\L  rA  B  | 

(IV.  S)  x(t)  =  Lj  cos  2*i«nt  +  -jft  sin  2irvr|t  |  , 

n=-°°  *-  1 

the  coefficients  An  and  being  given  by  the  expression 

f  T 

A  -  \  x(t)  cos  2wvt  dt  *•  C  cos  v> 
n  j.  n  a  n 


>n  -  $; 


x(t)  sin  2#»'nt  dt  *  sin  v»n* 


The  quantities  C  and  <p  or  A  and  B  are  pairs  of  random  co- 

n  -  n  n  n 

efficients  attached  to  each  oscillogram.  They  are  -such  that 


(IV. C)  An  =0,  Bn  =  0  and  A*  *  B*  - 

The  noise  X(r)  in  the  spectrum,  calculated  from  the  interfero- 
gram  I’(t),  is  given  by  the  expression 


(IV. 7) 


X(«0  B  Tcos l xCt )  ]  - 


It  is  the  cosim  transfoim  of  a  function  x(t)  which  can  be 
considered  a  sum  *>i  si  ml  waves  of  absolutely  .rbitrary  Irc- 
quenoios  n  and  phases  To  each  sine  wave  there  corre¬ 

sponds  in  the  spectrum,  the  F.T.  of  the  noise  oscillogram, 
a  line  which  is  the  convolution  of  the  Dirac  function,  of 
height  C^^T,  with  the  apparatus-function  corresponding  to 
the  phase  displacement  v*n-  Depending  on  the  phase  of  the 
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Fig.  38.  Spectrum  Obtained  From  the  Inter- 
ferogram  I’(t). 

corresponding  sine  wave  in  x(t),  the  apparatus-function  will 
have  one  or  the  other  of  the  shapes  represented  in  Fig.  30. 
X(«')  is  therefore  a  random  function  whose  radius  of  correla¬ 
tion  is  of  the  order  oi  1/T,  t f » . .  width  of  the  apparatus- 
function  (Fig.  3b)  ‘#). 

(#)this  when  one  has  made  an  apodication.  Without  apodi- 
zation  the  radius  of  correlation  of  the  noise  is  effectively 
divided  by  2. 
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3 .  SffUDV  OF  VARIANCE  AMD  S/W  RATIO  IN  THE  INTERFEROcrjtM 
AND  SPECTRUM 

We  shall  call  the  ratio  s/n  in  the  intorf crogram,  (s/n)j 
for  a.  line  of  frequency  i'  ,  the  ratio  of  the  amplitude  ft  of 

the  corresponding  sine  wave  in  the . interf crogram  to  the  mean 

/ 

square  error  a*  of  the  fluctuations  of  the  interf erogram. 


(IV. 8)  ax  =  /  x*(t),  (s/n)T=^/ax. 

ft  depends  only  on  the  luminance  of  the  source,  upon  the  width 
of  the  beam  and  upon  the  detector  used;  it  is  independent  of 
the  recording  speed. 

To  calculate  the  s/n  ratio  in  the  spectrum,  w<  are  going 
to  distinguish  two  cases,  according  to  whether  the  spectrum 
is  obtained  simply  by  making  the  cosine  transform  of  a  sym¬ 
metric  interf crogram  recorded  for  6  varying  from  0  to  L  during 
the  time  T  or  by  making  the  complete  Fourier  transform  of  an 
interf crogram  recorded  from  -  L  to  +  L  during  the  same  time  T. 

i  ■ 

(i)  Case  Where  the!  Touricr  Transform  is  Made  by 

Adjusting  the  Zero  Point 

Wo  have  seen  in  Chapter  2  that  the  height  of  the  correspond 
ing  line  in  the  spectrum  was  q/JT,  q  being  a  numerical  factor, 
the  mean  value  of  the  apo«.ii..ing  function  A(t).  It  equals  1 
when  there  is  no  apodisation,  0.5  with  the  weighting  function 
Aj  and  6.53  with  the  weighting  function  A^ .  We  shall  call 
the  signal/noisc  ratio  in  the  computed  spectrum  the  ratio 
(s/n)  ,  which  is  the  ratio  of  the  theoretical  height  of  the 
line,  v  il  enlaced  from  a  ■  inc.  wave  of  amplitude  ft  in  the  intcr- 
forogiau,  to  the  mean  square  error  a  of  the  fluctuations  in 
the  computed  spectrum. 


(IV. 9)  aX  i  =■  J  \J(iO  and  (s/n)s  «  qfiT/v* ( g } . 
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The  signal  thus  defined  in  the  spectrum  increases  proportion¬ 
ately  with  the  duration  T  of  the  measurement;  we  are  going  to 
show  that  <r  }  increases  proportionately  as  the  square  root 
of  T,  hence  that  the  s/n  ratio  varies  also  as  ^/T. 


A.  CALCULATION  OF  VARIANCE  IN  THE  SPECTRUM.  The  fluctua 
t-ions  due  to  noise  in  the  neighborhood  of  the  frequency  p^  in 
the  reconstructed  spectrum  have  the  form 

(IV. 10)  XM  =  2  C  x ( t )  A(t)  cos  t  dt. 

1  o  * 


From  this  one  deduces  that 


x  ( t )  x  ( t 1 ) 


A(t)A(t’)  cos  2irVjt  cos  2n»jt’ 


dt  dt* 


and,  noticing  that  in  this  expression  only  x(t)  and  x(t')  are 
rarijdom  variables 

(ivUi)  * 

■  CTfT _ 

4  y  \  x(t)x(t')  A(t ) A(t ’ )  cos  2wp  t  cos  2rp  t’  dt  dt'. 

^0  "  0  II 

The  function  x(t!)x(t 1 )  is  the  autocorrelation  function  of  the 
noise  in  the  interferogram.  Its  study  rests  entirely  on  a 
theorem  of  Khinchin:  the  cosine  transform  of  the  autocorrela¬ 
tion  function  of  a  stationary  random  variable  of  order  two  is 
its  power  spectrum.  It  follows  that 

P(0)  =  7?t7  x(t ’ )  *  y  ®  (v)  cos  2t*'0'd»»,  with  0  *  t  -  t*  and 

(IV. 12)  9  (p)  -  \  Pie)  cos  2wp$  d0, 

P  " .« 
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®p(v)  being  the  even  part  of  the  power  spectrum  c(i>)  a  ti.:- 
noise  contained  in  the  interferograin  (Tig.  35). 

From  Eqs.  (IV. 12)  one  gets  immediately  the  expressions 
giving  the  density  of  noise  and  the  variance  of  noise  in  the 
interferogram: 


(IV. 13)  ®p(0) 


p(0 )  dO  and  ff*  *=  p(0) 


®  (p)  dp. 

„«•  p 

Replacing  p(0 )  in  expression  (IV. 11),  we  get 
/■*  "+*°  /■*  In  T 

a*  *  4  V  *  (p)  d^  \  \  A(t )A(t ' )  cos  2ki>  t  cos  2*p  t  * 

- 1  * )  J.oB  P  w#  J0  *  * 

[cos  2jr»*t  cos  2*Pt*  +  sin  2*vt  sin  2ivt '  ]  dt  dt  ’ . 

Af ter _ expanding  and  eoi looting,  the  variance  can  be  written. 

,  i  *{  O 

(iv.  14)  a*(1)  “  i  j  +  +  .  \ 

+  f,(pJ+p)  +  fc*(»,1+p)l  dp, 

f(p)  and  v(i’)  being  the  apparatus  functions  in  cosines  and 
sines  del  ined  in  the  prt  ceding  chapter.  When  the  interferograin 
has  been  recorded  through  a  filter,  the  spectrum  ®(p)  of  the 
noise  that  it  contains  is  not  constant  as  a  function  of  the 
frequency.  But  one  can  consider  ®(r)  as  constant  and  equal 
to  ®(p  )  in  a  frequency  interval  which  is  large  with  respect 
to  the  si.-.e  of  the  function  fl.(v,),  and  «*(p).  Moreover, 

■  •  -Ho  p  -He  j  •  +• 

(IV.  15)  )  f*(p>  dp  *  J  *#(p)  dp  -  y*  A*(t)  dt, 

whence  the  new  expression  of  the  variance 
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Pig.  39.  Square  of  the  Apparatus-Function 
Exploring  the  Noise  Spectrum. 


(IV. 16)  «*X(I)  -p(%) 

Replacing  f*(>'J-i')  by  (see  (II. 7)),  we  get 

<w-17>  •!<>,  *  v-d  * 

-  S  »<•’»>  *’’’ 


a  result  which  may  be  stated  thus  (Fig.  39):  the  mean  power 
of  the  noise  in  the  spectrum  is,  within  a  factor  of  g  f*,  the 
chopped  energy  in  the  power  spectrum  of  the  noise  contained  in 
the  interf eroaram  times  the  square  of  the  normalized  apparatus- 
function  centered  about  i*  .  The  factor  l/2  comes  from  the  fact 
that  the  F.T.  made  under  the  conditions  considered  above  is  a 
synchronous  detection  and  one  finds  the  usual  gain  factor  l/2 
each  time  a  classic  detection  is  replaced  by  a  synchronous  de¬ 
tection  (49 J.  Indeed,  we  have  Seen  that  noise  in  the  inter- 
ferogram  can  be  put  in  the  form  of  a  sum  of  even  terms  and  odd 
terms.  When  the  cosine  transform  is  made  one  detects  only  the 
even  terms.  ^  -H# 

According  to  (IV.  15)  the  term  \  fs(v )  d *  can  be  j>ut  in 
the  for*  “** 
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f  *  C*' )  di> 


2  Q  T  , 


Q  being  the  mean  value  of  the  square  of  the  apodizing  function: 


Q 


_i_  r +T 
2T  j_T 


A* (t)  dt  . 


This  is  a  numerical  factor  which  is  equal  to  1  when  there  is 
no  apodizing,  0.33  when  the  inter ferogt'alm  is  weighted  by  a 
triangle  function  A4  and  0.406  with  the  weighting  function 
Aa ;  (IV. 16)  then  becomes 


(IV. 18)  <r*  ■*  2  Q  ®p(^1)  T  -  Q  . 

We  shall  see  in  the  following  chapter  that  the  filter 
is  always  chosen  in  such  a  way  that  the  spectrum  of  the  signal 
under  consideration  is  in  a  region  where  one  can  as  a  first 
approximation  identify  ®(rj)  with  ®(0)  (noise  density),  in 
order  not  to  have  excessive  phase  displacement  of  the  signal; 
under  these  conditions  a ^  can,  whatever  the  frequency  con¬ 
sidered,  be  put  in  the  form 

(IV. 19  <rx  -  /Q  /S?oy  /¥  . 


The  noise  density  ®(0)  is  a  quantity  which  depends  only  on 
the  detector  used.  Win  n  ...in:  changes  the  recording  speed  it 
is  necessary,  of  course,  t>..  modify  the  time  constant  of  the 
low-pass  filter  but,  when  doing  this,  one  modifies  only  the 
width  of  the  band  pass  of  the  filter  and  not  ®(0)  (Fig.  40). 
The  spectrum  (IV. 9)  then  becomes  , 


(IV. 20) 


(s/n)s 


■  .-affl--  *  a  _JL_  yr  . 

n'q*(o)t  v/q  y®(o) 


The  numerical  factor  q/x/Q,  which  is  the  result  of  dividing 

the  moan  value  of  the  weighting  function  by  the  square  root 
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1* 


Fig.  40.  Noise  Spectrum  for  Two  Values  of 
the  Time  Constant. 


of  the  mean  value  of  the  square  of  this  weighting  function,  is 
equal  to  1  when  there  is  no  apodizing,  0.87  with  the  triangle 
weighting,  and  0.84  with  weighting  by  the  function  A#. 

B.  DISCUSSION  OF  THE  RESULTS.  From  Eq.  (IV. 20)  giving 
the  expression  for  the  s/n  ratio  in  the  spectrum,  we  cam  draw 
three  important  conclusions: 

a.  When,  by  the  Fourier  transform  method,  one  wishes  to 
treat  a  given  problem  with  a  given  resolution,  the  s/n  ratio 
in  the  reproduced  spectrum  varies  as  the  square  root  of  dura¬ 
tion  of  measurement.  Thus  one  meets  again  an  absolutely  general 
result,  common  to  every  method.  By  making  the  comparison  be¬ 
tween  the  coefficients  of  proportionality  in  the  Fourier  trans¬ 
form  method  with  those  of  classical  methods  that  one  can  prove 
the  superiority  of  one  of  the  processes  over  the  others. 

Let  us  suppose  that  one  wishes  to  record  with  a  classical 
instrument  the  monchromatie  line  oi  negligible  width  which 
will  give  a  sine  wave  of  amplitude  p  in  the  interferogram.  We 
assume  that  the  instrument  has  the  same  extent  as  the  inter¬ 
ferometer  used  as  a  spectrometer,  that  the  transmission  factors 
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I'ig.  41.  Signal  Recorded  by  Fig.  42.  Spectrum 
a  Grating  Spectrometer.  of  the  Electrical 

Frequencies  Con- 

...  .  ...  talned  in  the  Re- 

_  corded  Signal. 

of  the  two  apparatuses  are  the  same* 10 >,  that  the  two  measure¬ 
ments  arc  made  with  the  same  detector  and  that  the  resolution 
is  the  same  in  both  cajcs.  Mndor  these  conditions,  the  signal 
recorded  as  a  function  of  time  reproduces  the  apparatus -function 
The  maximum  will  have  />*  at-  its  height.  Let  us  suppose  that  we  • 
have  taken  a  time  I  to  record  a  spectral  interval  equal  to  the 
limit  of  resolution,  that  is  to  say  equal  to  the  width  of  the 
apparatus-function.  In  Fig.  41  wo  show  the  apparatus-function 
oi  a  grating  spectrometer.  The  frequency  -spectrum  of  the 
electric  signal  thus  record* d  is  the  cosine  transform  of  it; 
it  extends-  from  -  l/T  to  4  l/T  (Fig.  42).  The  ideal  filter 
to  use  lot'  recording  under  these  conditions  is  a  low- pass  rec¬ 
tangular  filter  having  a  modulus  of  gain  equal  to  1  for  v  taken 
between  - .l/T  and  +  l/T  and  the  power  spectrum  •{»»)  of  the 


(,#,Tn  order  that  tt>  u an s pittance  be-  equal  to  one,  with 
the  Michelson  into rlr  rome ter  on.  ua  .u  use  the  exit  beam  as 
described  herein  and  th*.  i *.f  1>-oted  beam  which  is  sent  back  in 
the  direction  of  the  source.  It  suificcs  then  to  receive  the 
two  beams*  on  a  cell  and  to  take  the  difference  of  the  two 
electrical  signals  (16]. 
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it* 


fluctuations  in*" the  spectrum  is  thus 
determined  (Fig.  43).  In  fact  one  can¬ 
not  make  a  filter  whose  gain  curve  is 
that  of  Fig.  44  (curve  a).  One  can  use 
a  real  filter  having  the  same  area  (Fig. 
44,  curve  b)  and  the  power  of  the  noise 
in  the  spectrum  stays  the  same  even 
though  its  structure  may  be  different.  In  both  cases  the  vari¬ 
ance  in  the  spectrum  equals 


Fig.  43.  Ideal  Fil¬ 
ter  for  Use  With  a 
Grating  Spectrometer. 


(IV. 21) 


1 


®  ( v )  dp 
P 


C  +i/T 

J  ,  dv 

-l/T  p 


The  s/n  ratio  in  the  recording  then  takes  the  form 
s/n  *  0(1/  v/  »( 0) )  V/T .  This  expression  is  exactly  equal  to 
that  giving  the  s/n  ratio  in  the  computed  spectrum  using  the 
Fourier  transform  method  without  apodizing  (IV. 19).  Thus  if, 
it  is  a  question  of  studying  a  single  spectral  element,  the 
classical  methods  and  the  Fourier  transformation  method  give 
the  same  result. 


Fig.  44.  (a*)  Theoretical  Gain  Curve;  (b)  Gain 
Curve  of  a  Real  Filter  Having  the  Same  Area  as 
the  Theoretical  Filter  r  *  1/4  T;  (c)  Gain 
Curve  of  the  Filter  Used. 
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Notice  that,  with  the  real  filter  mentioned  above,  the 
high  frequencies  in  a  spectrum  undergo  inadmissible  attenua¬ 
tion  and  phase  displacement.  In  practice,  in  order  that  the 
attenuation  and  phase  displacement  of  the  high  frequencies  do 
not  exceed  a  certain  amount  which  has  been  fixed  in  advance, 
one  uses  a  filter  with  a  band  pass  which  is  much  larger  (Fig. 

44,  curve  c).  Then  the  power  of  the  noise  in  the  spectrum  is 
higher,  but  one  can  perform  a  smoothing  operation  on  the 
spectrum  and  for  all  intents  and  purposes  find  the  same  re¬ 
sults  as  before. 

When  the  classical  spectrometer  used  is  not  a  grating 
spectrometer,  the  spectrum  of  the  frequencies  is  no  longer 
necessarily  zero  for  frequencies  above  l/T  but,  the  record¬ 
ing  time  of  a  spectral  element _always  being  T,  the  maximal 
frequency  of  the  electrical  signal  that  one  may  transmit  with¬ 
out  deformation  is  always  of  the  order  of  l/T  and  the  preceding 
considerations  x  main  valid.  - 

Let  us  now  suppose  that  we  wish  to  study  a  spectrum  contain¬ 
ing  M  spectral  elements.  If  one  wishes  to  keep  the  same  s/n 
ratio  by  the  classical  method,  the  duration  of  the  measurement 
will  become  MT.  But  for  a  total  duration  of  measurement  equal 
to  MT,  the  s/n  ratio  in  ih*  spectrum  obtained  by  a  Fourier 
transformation,  according  to  (IV. 19),  s/n  *  /i  (l/ y/ »(0)  /MT . 

Hence  we  sec  the  gain  of  in  the  s/n  ratio  that  the  Fourier 
transform  method  permits  one  to  obtain  with  respect  to  classi¬ 
cal  methods  in  the  study  of  a  problem  with  a  given  resolution 
and  a  given  duration  of  time. 

b .  The  s/n  ratio  in  the  spectrum  is  independent  of  the 
time  constant  of  the  filter  through  which  one  has  recorded 
the  interf erogram.  This  is  a  result  to  be  compared  with  the 
fact  that  in  the  classical  methods  the  precision  of  the  meas¬ 
urement,  but  not  the  s/n  ratio,  is  independent  of  the  time 
constant  which  has  been  used  to  record  the  spectrum. 
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c.  For  a  measurement  time  equal  to  T,  the  s/n  ratio  in 
the  spectrum  varies  with  the  type  of  apparatus-function  chosen. 
It  is  maximum  when  no  apodizing  has  been  done.  It  is  lower  by 
a  factor  0.87  when  the  weighting  function  in  a  triangle  func¬ 
tion  and  by  0.84  when  it  is  the  function  .  Weighting  the 
interferogram  improves  the  shape  of  the  apparatus-function 
but  brings  about  some  losses  in  resolution  and  s/n  ratio. 

This  weakening  in  the  s/n  ratio  is  moreover,  slightly  lower 
than  what  one  gets  in  optics  by  using  absorbing  apodizers  in 
the  presence  of  photon  noise  and  a  fortiori  in  the  presence 
of  detector  noise.  The  simple  argument  which  consists  in  say¬ 
ing  that  in  the  case  of  numerical  apodizing  used  on  an  inter¬ 
ferogram,  one  does  not  diminish  the  s/n  ratio  because  the 
signal  and  the  noise  are  divided  by  the  same  factor  or  because 
the  number  of  photons  received  is  not  diminished,  is  there¬ 
fore  not  correct. 

c.  CALCULATION  OF  <R.  RESULT  OF  DIVIDING  THE  S/N  RATIO 
IN  THE  SPECTRUM  BY  THE  S/N  RATIO  IN  THE  INTERFEROGRAM.  <R,  the 
ratio  of  (s/n)g  to  (s/njj,  can  be  put  into  different  equiva¬ 
lent  forms  which  we  shall  discuss.  In  general 


(IV. 32) 


Theoretically,  the  interferogram  could  have  been  recorded 
without  a  time  constant,  without  the  s/n  ratio  in  the  calcu¬ 
lated  spectrum  varying.  Under  these  conditions,  the  mean 
square  error  <r^  of  the' fluctuations  in  the  interferogram  will 
b  infinite  and  <H  will  also  be  infinite.  -In  practice  the 
interferogram  is  always  recorded  through  filters,  and  now  we 
arc  going  to  show  that  the  quotient  <rx/<r^  is  nearly  propor¬ 
tional  to  the  inverse  of  the  square  root  of  the  time  constant 
of  the  filter  used. 


70 


NAVWEPS  RE! PORT  BQqq 


If  the  low-pass  filter  is  an  RC  filter,  the  variance  in 
the  spectrum  in  the  neighborhood  of  the  frequency  *»  ■  0  has 
the  form 

a 


•i,.,  ■  2  Q  T  y°>- 


(IV. 22b) 

The  variance  of  the  interf erogram  is  given  by 

1 


%  -  V0)  J_. 


l+(2jrvr)1 


.  V0> 


whence ( 41 } 


CIV. 23) 


X(0) 


4  Q  T  r  a* 


( 1 1  ^here-xs- a  very  simple  way  to  arrive  at  this  result 
directly  without  using  relations  (IV. 14)  and  (IV. 17)  which  bring 
in  the  convolution  of  the  noise  spectrum  with  the  square  of  the 
apparatus-function.  Although  the  calculation  is  simpler  using 
this  procedure,  the  first  way  of  looking  at  th^  problem  is  still 
very  useful  in  order  to  understand  it  physically  and  to  treat 
the  cas^  of  noise  in  the  numerical  F.T.,  as  will  be  done  at  the 
beginning  of  the  following  chapter. 

The  functions  y(t)  and  X(i>),  being  reciprocal  F.T.,  are 
related  in  the  following  way 

n  *  n  «° 

(IV. 24)  \  y*(t)  dt  -  \  X*(m)  di*  ... 

J0 


In  our  problem,  the  noise  interf erogram  has  been  weighted  by 
the  function  A(t)  in  such  a  way  that  y(t)  ■  x(t)  A(t),  x(t) 
being  the  noise  in  the*  interf  erogram.  The  fluctuations  X(iO 
in  the  spectrum  obtained  from  the  recorded  interf erogram with 
the  filter  of  gain  G(v)  are  related  to  those  that  one  would 
have  in  the  spectrum  obtained  from  an  interf erogram  recorded 
without  the  filter  Y(»*)  by  relation 

(cont.  on  following  page) 
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and  so  we  have 

=  *  -2- 

The  variance  in  the  spectrum  in  the  neighborhood  of  frequency 

y  has  the  form 

! 

- 2  Q  T  W  =  *5,..  V',»/yo>, 

-  > 

which  we  can  also  write  as 
-  _  1 


X(1)  1+(2wi»it)*  X(,) 


(l,i,cont) 

f  |x(v)|a  dr  -  f  |Y(r)|a  |G(r)|a  dr. 

Jo 

* 

Hence  (IV. 24)  may  be  written 

(IV. 25)  \  x*(t)  A*  (t)  dt  -  \  |X(r)|a  |G(r)|a  dr. 

o  Jo 

Let  v*  be  the  variance  in  the  interf erogram  and  cra  the  variance 
in  the  spectrum  in  the  neighborhood  of  the  zero  frequency. 
Noticing  that  Aa(t)  and  | G ( *0 ! *  are  quantities  which  vary  very 
slowly  compared  with  xa(t)  arid  Ya(r)  the  relation  (IV. 25)  be¬ 
comes 

CO  /*  «Q  p  *9 

(IV. 26)  <ra  j  A* (t )dt  -  <ra  J  |C(r)|a  dr  *  <rJ(0J  j  |G(r)|adr, 

since  the  li'JCvi  eon.  idc.r*  d  has  used  unit  gain  at  the  zero  fre¬ 
quency. 

Calling  Q  the  i.u  an  value  of  the  square  of  the  dpodizing 
function  and  replacing  |G(i*)l*  by  its  value  3/li+(2»rr)a I ,  one 
concludes  Q  ,T  ea  ■  aX(0;^4T’  w^enc®  *  *  $  *  *  **#  *n  0%m 

presslm  equivalent  to (XV. 23). 
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The  expression 
(IV. 23b) 


ff»  .  4.  ft  I  T_  „• 

X(l)  l+(2ri/|T)*  x 


allows  o>  e  to  calculate  the  variance  in  the  spectrum  in  the 
neighborhood  of  any  frequency  whatsoever  as  a  function  of  the 
variance  in  the  interf erogram. 

The  ratio  <U  varies  as  the  square  root  of  the  ratio  of  the 
duration  of  measurement  to  the  time  constant  of  the  filter  JT 
It  suffices  to  increase  T  in  order  to  make  it  as  large  as  one 
wants.  Let  us  take  the  case  where  the  F.T.  is  made  without 
apodizing:  q/  V'Q  then  equals  1  and  =  0 .  r>  \/ T,/ t\  ~ if 

T  =  4  x  104t  (which  results  in  an  interf crogi'am  recording  last* 
ing  three  hours  if  the  time  constant  used  is  r  =  1  sec), 

A  =  100.  We  deduce  from  this, that  in  order  to  have  the  s/n  =  2 
in  the  calculated  spectrum  it  is  enough  that  the  s/n  ratio  in 
the  interf erogram  by  l/sO.  The  corresponding  sine  wave  will 
be  totally  invisible  in  the  noise.  Hence  one  n^st  be  careful 
not  to  stop  recording  the  interf erograms  when  one  no  longer 
sees  fringes  and  when  one  has  the  impression  of  recording  only 
noise . 


D.  EXPERIMENTAL  VERiriC; ^IQHG .  Our  experimental  veri- 
irications  have  been  eoncertnU  with  three  precise  points.  By 
means  of  a  synchronous  detection  amplifier  and  low-pass  filter, 
described  previously,  we  have  recorded  the  fluctuation  of  a 
lead  sulfide  cell.  We  calculated  the  correlation  function  of 
the  fluctuations  x(t)  and  their  spectrum  by  making  the  cosine 
transform  of  the  autocorr-  lation  function.  Then  we  composed 
a  synthetic  interf erogram  by  adding  a  sine  wave  of  known 
period  and  amplitude  to  the  noise  under  consideration  and  we 
made  the  F.T.  of  this  interf erogram.  Thus  wc  have  been  able 
to  verify  the  agreement  between ‘‘the  predicted  s/n  and  the 
s/n  measured  in  the  spectrum,  as  well  as  the  variation  A  with 
the  s/n  ratio  when  one  passes  from  the  interferogram  to  the 
spectrum. 
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a.  Correlation  Function  and  Noise  Spectrum  in  the 
Interferogram.  We  wish  to  verify  that  no  other  sources ~of 
noise  have  been  introduced  other  than  those  considered  and 
that  we  would  be  able  to  calculate  from  the  expression  that 
we  had  when  the  interferogram  was  recorded  through  an  RC  filter 
We  chose  r  <*  2  sec  (Fig.  37a).  The  recorded  fluctuations  have 
a  mean  square  error  =  3  7.5  u,  u  being  the  value  of  an  arbi¬ 
trary  unit  which  servos  to  measure  the  deviations  on  the  re- 
•- „-i  ding(hcrc  millimeters).  The  correlation  function  was 
computed  by  the  IBM- 704  calculator  * 11 ) . 

If  the  recording  time  of  x(t)  were  infinite,  its  auto¬ 
correlation  function  p(0)  would  have  the  form  exp  (-t/r)  (Fig. 
40,  curve  a).  Since  the  recording  time  is  limited,  p(0)  is 
slightly  different  f-om  exp  (-t/r)  and  the  error-  increases  as 


(1,,It  is  possible  to  calculate  the  autocorrelation 
function  of  a  function  f(x)  which  is  known  to  have  a  bounded 
frequency  spectrum  of N width  Ar  using  equidistant  discrete 
Values  of  this  function.  It  is  sufficient  to  choose  the  dis¬ 
tance  l  between  two  points  in  such  a  way  that  the  repetition 
ircquoncy  l/'Zi-  of  the  spectrum  of  the  frequencies  of  the  new 
i unction  thus  defined,  is  equal  to  or  greater’  than  2  Av  If  the 
frequency  spectrum  of  T(x)  is  not  bounded,  the  approximation 
with  which  one  calculates  its  autocorrelation  function  by  this 
method  is  better  th«  smaller  h  is.  The  complete  discussion  of 
this  question  is  similar  to  the  one  that  will  be  made  in  the 
f  ollowing  chapter  to  drt.  vmiu<-  the  distance  7  between  two 
points  pin  ;d  from  the  inter i crogvaii*  in  order  to  have  the  maxi¬ 
mum  c  'n  ratio  when  one  makes  a  numerical  F.T.  Figure  45  shows 
a  function  !(.«.)  and  its  autocorrelation  function  CT(x)  computed 
numerically  using  discrete  values  of  f(x).  There  is  perfect 
agreement  between  the  calculated  and  predicted  Cf(x). 
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Fig.  45.  Illustrating  the  Calculation  of 
the  Autocorrelation  Function  Cf(x)  of  a 
Function  f(x). 


T/t  decreases.  Cui'vcs  b  and  c  of  Fig.  46  show  p(0)  for  T/r  * 

90  and  T/r  *  180.  Tfn.  fluctuations  are  considerably  reduced 
in  the  second  case,  as  predicted.  The  cosine  transform  of  p(0) 
gives  the  power  spectrum  of  the  noise  x(t)  (Fig.  47,  curve  t). 
It  departs  little  from  the  theoretical  curve  (Fig.  47,  curve  a) 

b.  5/N  Rat io  in  the  Calculated  Spectrum.  Wo  made  up  a 
theoretical  inci ri v  r ogram  I'(t)  (Fig.  ,37c)  formed  from  the 
superposition  of  the  portion  of  the  noise  x(t)  such  that  T/r 
«*  18CU-this  was  rxktis ickuvd  _fy*riiin*— and  of  a  sine  wave 

I(t)  -  37  COS  2*a.05t  (Pig.  3Tb);. 
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Under  these  conditions  the  s/n  ratio  in  the  interfero- 
gram  equals' .1  and  it  is  practically  impossible  to  distinguish 
which  of  the  two  graphs  37a  or  37c  is  pure  noise. 

The  spectrum  B^Cp)  of  I'(t)  has  been  numerically  calculated 
by  choosing  h  =1.2  (Fig.  38).  We  shall  see  in  the  following 
chapter  that  the  s/n  ratio  in  the  calculated  spectrum  is  then 
98/ltiO  of  the  maximum  s/n  ratio  that  would  have  been  obtained 
by  using  every  value  of  I’(t).  The  formula  (IV. 20)  permits 
the  prediction  of  (s/n)  .  According  to  Fig.  47, 

5 

®p(0.05)  »  4,950  u*sec,  where  (s/nj^  *  37  x  0.84  x 

x  - -  >/w5  «=  5.91. 

v/ 2X4,950 

The  s/n  ratio  measured -in-the  spectrum  equals  6,4.  It  has 
been  obtained  by  taking  the  ratio  of  the  measured  height  of 
the  line  at  v  -  0.1  to  the  mean  square  error  of  the  fluctua¬ 
tions  in  the  spectrum.  This  mean  square  error  is  the  mean 
value  of  the  square  of  the  fluctuations  in  the  neighborhood 
of  v  -  0.1.  The  noise  sample  has  been  chosen  sufficiently 
long  so  that  X(i>)  «  0.  Under  these  conditions  the  calculated 
equals  1,280  u  sec.  Notice  that  the  formula  (IV. 19)  allows 
one  to  predict  . 

«  •/  0 .406  x  2  x  4,950  x  365  ■  1,200  u  sec 
The  theoretical  signal  has  the  value 

s  *  q  0  T  -  0 .  £,  3  x  5  ?  :<  ji',0  -  7,059.6  u  sec 

It  is  probable  that  the  measured  s/n  ratio  lies  between  the 
two  limits  (s  +  o^)  *  6.R  and  (s  -  Q^/°y^  “  4.9.  ‘  We  found 

it  equal  to  6.4. 
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c.  Variation  in  S/W  Ratio  in  Passing  Froa  the  Inter- 
ferogram  to  the  Spectrum.  The  s/n  ratio  in  the  interf erogram 
being  equal  to  1,  the  measured  value  of  *  is  equal  to.  £.4. 

Its  theoretical  value  can  be  calculated  using  formula  (IV. 22) 
which  gives 

«  -  0. S3  .X  2^00  X  360  "  5.905.  . 

(2)  Case  Where  the  Fourier  Transform  is  Made  Without 

Locating  the  Zero 

We  showed  (§III.2(3)C)  that  it  was  possible  to  reconstruct 
the  spectrum  without  knowing  the  zero  point  of  the  interfero- 
gram.  It  is  enough  to  record  the  interf erogram  for  6  varying 
between  -  L  and  +  L  and  to  take  the  square  root  of  the  sum  of 
the  squares  of  the  sine  transform. and  the  cosine  transform,  in 
other  words  to  make  the  complete  F.T.  We  propose  to  study  the 
signal/noise  ratio  in  the  computed  spectrum  by  this  method  and 
to  compare  the  results  with  those  obtained  in  the  case  where 
only  the  cosine  transform  is  used.  When  one  takes  squares,  one 
introduces  nonlinear  elements  and  one  can  predict  that  it  will 
not  be  sufficient  to  consider  separately  the  F.T.  of  the  Inter* 
ferogram  without  noise  and  the  F.T.  of  the  noise;  there  will  be 
signal- noise  interaction. 

We  shall  suppose  that  in  the  two. cases  the  measurement 
has  lasted  the  same  total  time  T.  In  the  absence  of  noise  the 
two  methods  are  strictly  equivalent.  But  the  study  of  the 
signal-noise  ensemble  is  mueh  more  complex  when  one  takes  a 
complete  F.T. 

In~the  first  place,  the  noise  in  the  calculated  spectrum 
has  a  Caussian  distribution.  Its  mean  value  is  zero.  Wc  have 
called  "signal"  the  height  s  of  a  line  obtained  in  the'*  absence 
of  noise  taken  from  a  particular  interf erogram.  We  designate 
by  fi  the  variance  ej  in  the  spectrum  in  the  presence  of  noise 
(a  quantity  that  one  can  measure  in  a  spectrum  by  making  the 
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Fig.  48.  Appearance  of  a  Spectrum  Calculated 
by  Taking  the  Cosine  Transform  of  an  Intcr- 
ferogram  in  Which  the  Zero  Point  Has  Been 
Located. 

mean  square  of  the  fluctuations  in  a  region  where  there  is  no 
line.)  (Fig.  48).  There  is  a  definite  probability  P  that  the 
number^  measuring  the  height  of  the  line  lies  between  the 
s  +  \-rp  and  c  -  -/P,  for  any  s.  In  other  words,  the  operation 
is  perfectly  linear  and  the  mean  square  error  of  the  fluctua¬ 
tions  in  the  spectrum  is  independent  ol  the  signal.  Figure  51 
shows  the  variation  A/  y,rp  as  a  function  of  s/x/p.  Curve  a 
represents  the  mean  value  of  the  signal  Ji^p  *=  s  •/*>•  The 
curves  a '  ,$  /  \1>  =  (s/  \  P)  +  1 ,  and  a"  /  </p  *-  (s/  v/p)  -  1, 
mark  the  two  limits  between  whidh  there  is  the  probability  P 
of  finding  the  value  of  the  measured  signal  divided  by  sfp .  We 
have  called  the  signal/noise  ratio  s/n  the  quotient  s /  \fp . 

If  the  complete  Fourier  transform  is  taken,  the  results 
are  different;  in  order  to  compare  the  two  methods  more  easily, 
we  are  going  to  express  the  quantities  that  we  shall  call 
signal  and  noise  in  this  ease  as  a  function  of  s  and  p.  Being 
given  the  signal-noise  interaction,  one  cannot,  when  making 
the  complete  Fourier  transform,  deal  with  the  noise  alone  and 
then  with  the  sigr.al-noise  ensemble.  It  is  necessary  to  treat 
tire  signal  and  the  noise  together  and  from  them  to  deduce  what 
happens  in  those  regions  of  the  spectrum  where  the  signal  is 
zero.  Here  are  the  principal  results  that  we  put  in  evidence: 
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Fig.  49.  Spectrum  Calculated  by  Taking  the  Com¬ 
plete  F.T.  of  the  Same  Interf erogram  as  in  Fig.  48, 
but  in  Which  the  Zero  Point  Has  Not  Been  Located. 

1* .  In  regions  of  the  spectrum  where  there  is  no  line  (Fig. 
49,  section  ab) ,  the  noise  is  a  random  function  J(iO  that  one 
can  consider  as  the  superposition  of  a  constant  nonzero  com¬ 
ponent  J(i>)  *  i  and  of  a  fluctuating  part  z’(iO  of  mean  value 
zero  z  M75  *  0  ‘and  of  mean  square  value  z  (v) ; 

2* .  The  number  measuring  the  height  of  a  line  in  the  repro- 
duced  spectrum  no  longer  will  be  simply  a  number  taken  between 
s  +  a  and  s  -  a,  a  being  a  mean  square  error  to  be  determined. 
Each  value  &  can  always  be  put  in  the  Torn  ...  0 

(IV. 27)  i  ■  m  +  z  , 

m  being  the  mean  value  of  $  and  z  being  a  random  function  of 
mean  value  zero.  But  the  important  fact  is  that  m  no  longer 
is  equal  to  s  and,  moreover,  is  a  function  of  p.  On  the  other 
hand,  "z*"  is  likewise  a  function  of  s;  the  mean  value  m  becomes 
an  approximately  linear  function  of  s  when  s/ y/p  1. 

Let  us  state  the  meaning  of  m.  If  one  calculates  an  in¬ 
finity  of  spectra  using  an  infinity  of  interf erograms  recorded 
under  the  same  conditions  (constituting  an  infinity  of  experi¬ 
ments  belonging  to  a  single  category),  the  number  a  measuring 
the  signal  corresponding  to  a  certain  frequency  v  in  each 
spectrum  group  themselves  about  the  value  m,  mathematical  ex¬ 
pectation  of  J,  i.e.,  E 
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Wc  shall  call  "signal”  in  the  case  of  the  complete  F.T. , 
the  difference  m  -  i,  i.e.,  the  height  of  the  line  above  the 
mean  value  of  the  noise  (Fig.  49)  because  it  is  effectively 
this  quantity  which  will  be  accessible  to  the  experimenter. 
The  signai/noise  ratio  in  the  calculated  spectrum  is  measured 
by  the  quotient  (m  -  j)/ y/z'* . 

A.  STUDY  OF  THE  SIGNAL-NOISE  ENSEMBLE.  We  shall  study 
the  signal-noise  ensemble  at  eqch  stage  of  the  reproduction 
of  the  spectrum  by  assuming  that  the  sine  wave  in  the  inter- 
ferogram  corresponding  to  the  frequency  under  consideration 
has  a  displacement  if>6  with  respect  to  the  origin,  as  was  as¬ 
sumed  in  the  preceding  chapter.  The  final  result  clearly 
ought  not  to  depend  on  0 .  Wc  shall  then  study  successively 
the  results  of  the  cosine  ..transform  of  the  interf erogram,  of 
its  sine  transform,  and  of  the  square  root  of  sum  of  the 
squares  of  these  two  transforms. 

a.  Cosine  Transform.  The  height  of  the  line  resulting 
from  the  cosine  transform  of  the  interf erogram,  if  there  is 
no  noise,  will  be  (§11.2(3)0 

(IV. 28)  sc  *  s  cos 

The  corresponding  value  of  a  claculated  spectrum  s  can  be 

C 

put  in  the  form 

(IV. 29)  *  s  coo  l  a O'j). 

We  have  seen  before  that  the  noise  recorded  on  the  interf oro- 
gram  can- be  considered  as  a  supevpos  jti.-n  of  sine  waves  of 
random  amplitude  and  phase  and,  ye  can  write  that ..tho_noise 

component  at  the  frequency  v  h^s  the  form  ( 1 4 7 J  p.  31»1) 

Ci 

y(t>  .  cos(2ri»1t-v>) 

C  C 

“  2T  C0S  ^  COS  +  2T  sin  ^  sin 
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in  such  a  way  t^hat  the  result  of  the  cosine  transform  of  y(t) 
is  written  '  . 

X^ )  =  A  cos  tp. 

This  is  a  random  variable  whose  characteristics,  mean  value 
and  mean  square  error,  are  the  same  as  those  of  the  quantity 
X(v),  the  cosine  transform  of  x(t)  for  the  different  values  of 

Then  one  has 

(IV. 30)  A  cos  <p  ~  0  and  tr*  =  A*cos*  ^«A*/2  «*  p. 

A  (  1  ) 

Let  us  suppose  that  we  can  record  an  infinity  of  inter- 
ferograms  under1  the  same  conditions  and  that  we  plot  on  a 
graph  the  values  (ty  measured  from  each  spectrum  obtained  from 
each  interf erogram.  The  curve  joining  these  points  will  have 
the  general  shape  shown  in  Fig.  50a.  The  values  of  group 
“themselves  about  the  mean  value  ts  cos  the  mean  square  error 
of  fluctuations  being  equal  to  y/p. 

b.  Sine  Transform.  Obviously  the  results  are  analogous, 
and  in  the  same  way  one  can  define 

(IV. 31)  $s  •  s  sin  ?0  +  A  sin 

with 

A  Jin  'V  »  0  and  A*  sin*  <p  *  p  (Fig.  50b). 

i- — I 

<13)The  means  o.  ,X(»-  )  will  be  tal.«- n  from  an  infinity  of 
trials,  i.e.,  an  infinity  of  values  X(i>  ')  obtained  by  making 
the  cosine  transform  of  an  infinity  of  interferograms  recorded 
under  the  same  conditions.  One  has  the  relation  x(*  ^  *  xT*7 
■  0. 
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Fig.  50.  Study  of  the  Signal-Noise  Ensemble 
When  the  F.T.  is  Taken  Without  Locating  the  Zero  • 

Point  of  the  Interferogram. 

The  two  random  functions  are  «H  cot-  *p  and  ill  sin  <p  have  the-  same 
.uean  value,  same  mean  square  error,  name  autocorrelation  func¬ 
tion,,  but  their  coefficient  of  correlation  is  xero(l4>* 


* 14 * Their  correlation  function  is  alto  tcro  for  a  zero 
delay;  but  the  variables  41  sin  and  cos  v  arc  not  independ¬ 
ent  variables  and  their  correlation  function  is  different  from 
zero  tor  a  delay  which  is  not  zero,  which,  however ,  doesn't 
happen  in  the  present  case. 
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c .  Sum  of  the  Squares  of  the  Sine  and  Cosine  Transform. 

The  square  of  Sc  has  a  constant  component  and  a  fluctuating 
component  of  mean  value  zero.  '  The  two  quantities  which  are 
interesting  to  consider  are  the  constant  component  which  equals 

5*  and  the  variance  of  the  fluctuating  part,  i.e.,  (£*  -  i*)* 

(Fig.  50c): 

»  s*  cos*  +  «l*  cos*  (p  f  2  s  «  cos  *#  cos 

whence 

(IV. 32)  $*  -  s*  cos*  4>0  +  p. 

Thus  the  mean  value  of j$*  depends  not  only  on  s  and  on  but 
also  on  p,  hence  on  the  density  of  noise  in  the  inter ferograa. 
The  fluctuating  part  has  the  value 

-  £*  *  <K*  cos*  <p  -  p  +  2  s  *  coa  ##  cos 

.1 

Its1  variance  is  then— 

(IV.  33)  <a;-JJ)*  »  «4cos4v»  +  P*  +  4  s*cos**0  6*003*?  -  2p  4* cos*?. 


One  can  evaluate  the  term  «4  cos4  </>  using  the  following  clas¬ 
sical  theorem.  1 

If  x(t)  is  a  Laplacian  stationary  random  function  of  mean 

value  zero  with  autocorrelation  function  C(t),  i.e.t  such  that 

x*  *  C(0),  the  autocorrelation  function  of  its  square  x*(t) 
can  be  written  '51) 

x*(t)  x2(t-r)  »  2  C*(r)  +  C*(0). 

From  this  one  deduces  that 

x4(fc)  *  3  C*(0). 

i 
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In  our  case  since  <Ra  cos*  <p -*  p,  «4  cos4  9  *  3  p1,  (IV.  33) 
takes  the  form 

(IV.  34)  (<?*  -  £*)*  «  2  p*  +  4  p  s*cos  %o. 

The  fluctuations  will  be  the  greater  for  greater  p  and  also 
for  higher  signal.  The  relations  (IV. 32)  and  (IV. 34)  show 
the  signal-noise  interactions  very  well. 

In  the  same  way  one  calculates  the  constant  component  of 
$*  and  the  variance  of  its  fluctuating  component: 

(IV.  35)  =  s*  sin*  if. 0  +  p  and 

"  1 S  '  1 

=  2p*  4  4ps#sin%#  (Fig.  50d). 

2  1 1 2  V  s 

In  forming  the  sum  s  ~^c  +*s»  the  constant  components  are 
added  together  while  the  fluctuating  components  add  their 
variances1.  |  Hence  one  has  -  { 

(IV. 36)  £*  ■  s*  +  2  p  and 

«$*  -  J*)*  -  4 p*  -f  4ps*  (Fig.  50e) . 

d.  Calculation  of  eft.  In  the  las stage  of  the  calcula¬ 
tion,  the  square  root  of  $*  remains  to  be  taken;  we  propose 
to  determine  the  mean  value  of  $,  say  m  and  /?  the  mean 
square  error  of  the  fluctuating  component  of  if  as  functions 
oi  s  and  of  p. 

In  order  to  determine  those  quantities,  we  are  going  to 
make  use  01  the —two  relations  (IV.  36)  in  which  we  aret  going 
to  replaced*  and  A*  by-  the  ir  values  as-f  unctions  of  m  and  z 

(IV.  37)  ^*  =  m*  +  2ms  +  a*  -  m*  +  z*  (Fig.  50f) 

and  according  to  (IV. 36), 
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(IV.  38)  m2  +  2*  *  s2  +  2  p , 

(IV. 39)  (J2  -  &2)2  =  [)?  -  (s2  +  2  p))2 

-  V  -  2  v*(s*+2p)  +  (sa+2  p)2  =  V4  -  (s2+2  p>2. 

After  replacing  ^  by  its  value  as  a  function  Of  m  and  z  we  get 

Ft  T  _  4  ,  ,  2  a  ,  T 

We  can  substitute  for  m4  and  m2  their  values  as  functions  of 
s '  eno  p  taken  from  (IV.  58)  and  notice  that  "zT  *  3  (z  )*.  The 
second  equation  (IV. 36)  can  then  be  written 


(IV. 40)  (z2)2  -  z2(A  p  +  2  s*)  +  2  p2  4  2  p  3*  -  0, 

from  which  one  deduces 

_ 

(IV. 41)  z*  ■=  2-p  +  s2  -  1  vV  -  2  p  s*  +  2  p* 

t 

and  i 


(IV. 42) 


ra 


*  (s4  +  2  P  s2  +  2  p*),//4. 


The  results  arc  gathered  together  in  Fig.  51  (curves  b,  br,  b’’) 
in  which  w c  represented  ;v/ \  P  as  a  function  of  s / vP  >  i.e.,  of 
the  Lignal/noisi  ratio  in  the  spectrum  which  would  have  been 
obtained  simply  by  making  the  cosine  transform  of  the  same- 
interim  rogram  in  which  one  would  have  located  the  zero  point. 
Curve  b  indicates  the  variation  of  m/ \/p>  as  the  function  of 
a/  *  jp  ..ana  the  curv.  b'  and  b"  the  limits  between  which  there 


i.>  a  «i*'i  inite  pin- Lai: ility 


.  iVi  J.  » j 


the  measured  value 


2/  v  i 


W.  sou  immediately  that,  in  the  presence  of  noise,  m  is 
not  a  linear  function  of  s.  To  be  precise,  we  arc  going  to 
begin  by  examining  the  two  extreme  cases,  one  where  the  signal 
is  zero  (s  *•  0)  and  one  where  s /yp  is  large'. 


K 
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1*.  s  *  0.  The  relations  (IV. 41)  and  (IV. 42)  are  written 

z*  *  2  p  -  p  -/2  *  P  (2  -  /I)  and  m  -  2*/4  y/p. 

This  means  that  in  a  place  in  the  spectrum  where  there  is  no 
line  (Fig.  49,  portion  ab),  one  sees  noipe  made  up  of  a  con¬ 
stant  component  which  equals  2l/4v'p  and  oi  a  tluetuating 
*'<’\mponent  which  has  mean  square  error  (2  -  v/2)l//*vf>  ~  0.8  v’P> 
i.e.,  4/5  of  the  mean  square  error  of  the  fluctuations  in  the 
'spectrum  calculated  by  taking  the  cosine  transform  of  the  inter 
feroyram  recorded  from  0  to  I*. 
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2°.  s/ J~o  is  Large .  If  p/s a  can  be  considered  as  first-order 
inf initestimal,  the  relations  (IV.  41)  and  (IV. 42)  become,  by 
neglecting  terms  cf  the  second  order ,  • 

zJ  *  2  p  +  s*  -  s2(l  +  (2  p/s2)]l/a  =  P\  m  =  s  +  (p/2s), 

......  *  j  '  * 

Hence  for  lines  which  would  appear"with  a  large  s/n  ratio  in 
the  spectrum  calculated  by  the  cosine  transform  of  an  inter- 
ferogram  in  which  zero  phase  difference  is  located,  the  complete 
F.T.  gives  an  approximately  linear  result  and  the  mean  s.quare 

error  of  the  fluctuations  which  affect  $  is  then  the  same  as 

I. 

in  the  first  method. 

3° .  Values  of  s/  v ^p  in  the  Neighborhood  of  1.  It  remains  to 
examine  the  intermediary  cases,  for  values  s /v/p"  in  the  neighbor 
hood  of  1.  The  signal  m  -  j  is  weak  and  does  not  vary  propor¬ 
tionately  to  s  at  all. 

To  s/n  ratios  of  2,  3,  6  in  the  spectrum  that  one  would 
obtain  by  taking  thf  cosine  transform  (lines  of  frequencies 
v  ,  v  ,  v  in  Fig.  48)  correspond  t?ib  signal/noise  _ratios  of 
1,  2,  5  in  the  spectrum  that  one  would  obtain  by  taking  the 
complete  F.T.  (frequencies  vti  v  ,  v  in  Fig.  49). 

B.  EXPERIMENTAL  VERIFICATIONS.  These  have  been  carried 
on  only  in  the  case  where  s  =  0.  They  have  been  made  using 
the  recording  of  noise  x(t)  of  which  one  part  is  represented 
by  Fig.  37a  and  which  has  beert  used  for  the  experimental  veri¬ 
fications  of  the  preceding  chapter  (§111.1(3));  and  <8S, 
cosine  and  sine  transforms  of  x(t),  are  represented  by  Figs.  52a 
and  52b.  The  calculated  band  width  equals  0.1  cps,  i.e., 

20  times  the  correlation  radius  of  r^iso'.  The  values  of  the 
crests  lie  between  ±  8  u  see  (u  being  an  arbitrary  unit)  and 
the  mean  square  errors  measured  have  respectively  the  values 
v'pT  -  4.25  u  sec  and  *  4.20  u  sec.  The  functions  and 
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Pig.  52.  Study  of  the  Complete  F.T.  of  an 
Xnterferogram  Containing  Only  Noise. 

$*  (Fig.  52c  and  52d)  have  respectively  the  mean  valued*  “ 

18.06  u*  see*,  and  4^  “  17.79  u*  sec*,  very  similar  values. 

If  $c  and  $s  had  been  calculated  on  a  frequency  band  of  infi¬ 
nite  width,  the  variances  £*  and  would  be  exactly  equal. 

The  mean  value  of  the  sum  f*  *  <3*  +  ^*  *  35.85  ua  rec*. 

From  the  conclusions  of  the  preceding  paragraph  the  pre¬ 
dicted  value  of  the  variance  of  the  fluctuations  of  i*  -J* 
would  be  $*  -  written  t&*  -jJ*)*  *  2  p*  +  2  pj  •  1,296  u*  sec*, 
The  same  quantity  measured  in  Fig.  52c  equals  1,240  u*  sec*. 
The  predicted  value  of  the  mean  of  C?  is  5  ■  m  *  (/» v^)1^*  * 

4.9  u  sec  and  that  of  the  mean  square  error  of  the  fluctua¬ 
tions  of  J  -  m  is 

/ (J~-  m)*  *  yp\ 2  -  •/?)  ■  3.2,u  sec. 

The  same  quantities  measured  in  Fig.  S2d  have  the  values 
a  -  5.2  u  sec  and  </(£  -  •)"  ■  2.9  u  sec. 
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C.  CONCLUSIONS.  We  are  now  in  a  position  to  compare 
the  two  methods  more  completely  than  in  the  preceding  chapter 
The  second  method  is  never  strictly  linear  (in  the  same  way 
that  a  classic  detection  is  never,  in  the  presence  of  noise, 
rigorously  linear).  It  is  catastrophic  for  the  low  values  of_ 
signal/noise  ratio  but,  for  high  values  of  the  signal/noise 
ratio,  the  two  methods  are  practically  equivalent. 

One  can  show  that,  if  the  Fourier  transformation  is  made 
with  a  harmonic  analyser  (a  method  which  does  not  necessitate 
referencing  iiero)  the  study  of  the  linearity  and  of  (s/n)s 
ratio  leads  to  the  same  results  as  when  the  complete  F.T.  is 
made.  On  the  other  hand  let  us  recall  the  two  inconveniences 
of  this  last  method:'  The  displacement  of  the  moving  part  is 
'doubled  and  the  calculation  time  is  more  than  quadrupled. 


\ 
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■  V.  NOISE  IN  NUMERICAL  FOURIER  TRANSFORMATIONS. 

DIFFERENT  METHODS  FOR  OBTAINING 
THE  MAXIMUM  S/N  RATIO 

■I 

i 

When  the.F.T.  is  calculated  numerically  from  a  finite 
number  of  values  of  the  interf erogram,  the  s/n  ratio  in  the 
calculated  spectrum  is  always  lower  than  the  maximum  s/n  ratio 
that  can  be  obtained  by  using  a  continuous  succession  of  values 
of  I(t).  We  shall,  establish  the  e-.pression  giving  the  mean 
square  error  of  the  fluctuations  in  the  numerically  calculated 
spectrum;  we  shall  treat  the  important  question  of  the  choice 
of  the  time  constant  to  use  during  the  recording  and  we  shall 
discuss  diifci  i  nt  procedures  that  can  be  used  to  obtain  a  maxi¬ 
mum  s/n  rulio  from  a  given  interf erogram. 

1.  NOISE  IN  A  NUMERICAL  F.T.. 

Let  us  recall  the  fundamental  result  of  the  preceding 
chapter:  The  mean  power  of  the  noise  in  the  calculated  spec¬ 

trum  is  the  chopped  energy  in  the  spectrum  of  the  noise  contained 
in  the  interf  crop-*  ran  times  the  square  of  the  apparatus-function: 
now,  when  one  makes  a  numerical  F.T.,  we  have -seen  in  Chapter  I 
that  tin:  apparatus-function  is  composed  of  several  peaks  the 
distance  between  them  b«_  ing  1/fc  (Fig.  13),  /.*  being  the  interval 
of  the.  grating  lunction  us-  d  to  pick  points  from  the  interfero- 
gram;  the  expression  (IV.lv) 

<■»  +» 

aX  J  !*<%->  d* 

then  K  comes,  in  the  ease  of  the  numerical  Fourier  transforma¬ 
tion, 

(V.l)  ax*1  ^  \  3  (i>) 

One  thus  sees  that  several  peaks  of  the  app  ratuo-f unction 
explore  the  noise  spectrum,  even  in  the  ease  where  the  spectrum 
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Fig.  53.  Different  Peaks  of  the  Apparatus* 
Function  Exploring  the  Noise  Spectrum. 


to  ho  studied  is  limited  in  such  a  way  that  it  may  not  be 
explored  except  by  one  peak.  If  the  noise  spectrum  »(»/) 
contained  in  the  interf eregram  were  constant,  the  mean  power 
of  the  noise  in  the  calculated  spectrum  would  increase  in¬ 
definitely  with  the  number  of  peaks  in  the  apparatus-function, 
hence  the  absolute  necessity  of  using,  during  the  recording, 
a  filter  which  limits  the  noise  spectrum  in  the  interf erogram. 
Wo  shall  consider,  in  what  follows,  the  case  of  the  RC  filter 
(it  is  difficult  to  make  more  efficient  low-pass  filters  for 
the  very  low  frequencies  used  in  our  experiment);  the  noise 
spectrum  contained  in  the  interferogram  then  has  the  fora 

*  (p)  *  ®  (0)/ll  +  (2»vt)*J,  \  with  r  -  RC 

t  *  • 

and  in  the  neighborhood  of  a  frequency  in  the  calculated 
spectrum,  the  mean  poweF  of  the  noise  is  the  sum  of  the  chopped 
energies  in  the  spectrum  ®p(»0  multiplied  by  the  series  of 
poai  s  corresponding  to  the  square  of  the  apparatus-function 
.  ntcred  about  i> Q  ±  (m/fc)  (Fig.  53).  One  concedes  that  one 
will  have  used  nearly  all  the  information  contained  in  the 
recorded  interferogram  if  one  chooses  the  interval  ft  (measured 
in  units  of  time)  close  to  the  time  constant.  We  are  going 
to  state  this  idea  precisely  in  studying  the  variance  as  a 
function  of  interval  ft. 
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(1)  Mean  Power  of  Noise  in  a  Numerical  F.T. 

By  replacing  *p(0  by  its  value  in  Eq.  (V.l)#  we  get 


(V. 2) 


yo> 


m=-+* 


TT.S 


-Ho 


i+j  2^  [ 


f,(i»)  d»». 


Again  we  find  in  the  expression  for  <x*  the'  product 


"x,.,  =  *p(0)  '  .  d" 


which  gives  the  variance  in  the  calculated  spectrum  as  an 
integral,  in  the  neighborhood  of  the  0  frequency  which  we 
have  already  seen  in  the  preceding  chapter  to  be  (TV. 22b) 


(V,3)  ffx( o )  *  2  Q  T  ®p(0). 

In  the  case  of  a  numerical  F.T.  one  can  put  la  the  fora 

(V.4)  V*’T>’  . 


by  putting 
(V.  5) 


It  is  necessary  to  evaluate  this  series.  We  can  do  this  in 
a  particularly  easy  fashion  because  it  happens  that  the  F.T. 
of  the  function  l/Jl  +  (2irpr)iJ  is  easy  to  Calculate.  The 
calculation  is  based  upon  the  following  property:  let  P(r) 
and  Q(t)  be  two  functions  which  are  the  F.T.  of  each  other. 
Then  there  exists  between  them  the.  relation  . 

n*)  &  -  q(o). 
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»; 


<«*■"  t  -•*  * 

Fig.  54. --(a)  Grating  Function 
***/*’  Grating  Function  6^. 


the  Dirac  function  centered  about 
We  put 


Notice  that  the  series  (V.5) 
is  none  other  than  the  inte¬ 
gral 


,  (V. 6)  *  « 


i 

1+(2iti>)* 


Ri/k 


being  the  convolution  of 
a  Dirac  distribution  of  periodic 
support  of  interval  l/fc  with 
the  frequency  r  »  (Fig.  54a). 


P(»0  * 


1  +  (2»*r) 


7  V* 


00  -  »^(iO 


In  order  to  calculate  * 


1 


P(t>)dv  it  then  suffices 


to -find  the  F.T.  Q(t)  of  P  ( »* )  and  to  take  the  particular  value 

Q(0). 

The  F.T.  of  P(v)  is  the  convolution  of  the  F.T's  of  •£(*) 
and  of  The  F.T.  of  ®A(»0  is  known;  it  is 


cos 


1+(2»»»t) 


U-Le-t/r 

2r 


The  F.T.  of  R^fiO  has  the  expression 

TcoslRi m  *  **  cos  2ee#t, 

R.’  being  a  Dirac  distribution  of  periodic  support  of  interval 
h  (rig.  54b).  From  this  one  deduces  immediately  the  value  of 


*  "  2t  *  ^r 


ijJ  cos  2»i'#t, 
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that  one  can  put  in  the  form 

*  ■  *  [Re  (2  2  * 

u  m*0 


s  tip  .ink 


-0] 


or 

(V.  7) 


t-Afr  ' 


2r 


L 


Re  i 


2 


l-exp{-fc  [  (l/r)-2»i»» 


J5 '  >)]• 


After  expanding  and  collecting,  (V.7)  becomes 
(V.  8) 


a  _  JL  sinh(fc/T) _ 

2r  cosh (h/r)- cos 2tp^h  * 


from  which  the  expression  for  the  mean  power  of  the  noise  in 
the  case  of  a  numerical  F.T.  is 


(V.9) 


ax  "  ®p(0)  2  Q  T  • 


The  expression  (V.8)  of  *  (k,r)  permits  one  to  calculate  the 

% 

variations  of  the  power  of  the  noise  as  a  function  of  the 
interval  h ,  when  the  recording  has  been  made  with  a  certain 
time  constant  r.  In  the  particular  case  where  v  ■  0, 

##(fc,T)  then  becomes  a  function  only  of  the  quotient  k/r 
(Fig.  55,  curve  a): 


(V. 8b) 


*0(*A) 


h  sinhC^/T)  . 
2r  cosh(k/r)-l  * 


h/r  ( h  measured  in  time  units)  is  the  inverse  of  thej  number 
of  points  taken  from  the  interferogram  during  a  time;  equal 
to  r. 

The  curve  ♦  (J’/r)  has  a  horizontal  tangent  at  the  origin. 
For  v  -  0  and  fc/t  =  0,  *0(h/'r)  =  1  and  again  we  encounter  the 
familiar  expression  for  the  variance  corresponding  to  this 
case 

-  *  yo) « *  • 


95 


MAVWEPS  REPORT  8098 


Fig.  55.  Study  of  the  Variance  in  the 
Spectrum  as  a  Function  of  the  Ratio  of  the 
Interval  to  the  Time  Constant:  Curves  b, 
v%  •  0;  Dashed  Line,  ■  (3/2)  (3/2 *r). 

When  the  analyzing  frequency  is  an  arbitrary  frequency 
»*#  different  from  0,  one  can  define  a  function  *(*/t)  valid 
for  a  certain  value  of  the  product  i >0r.  It  can  be  written 

<V.  8c>  * v  fy  -  &  c,h(iM-coI?a...T.7vy  • 

This  gives  the  variations  of  the  (.power  of  the  noise  in  the 
spectrum  in  the  neighborhood  of  the  frequency  »9,  as  a  func¬ 
tion  of  the  quotient  fc/r  when  the  interferogram  has  been  re¬ 
corded  through  a  filter  with  time  constant  r.  Curve  b  of 
Fig.  55  has  been  constructed  for  *  l/4»,  i.e.,  the  fre¬ 
quency  i*#  considered  equals  1/4»t;  it  is  such  that 
•  (•»  )  -  0.8  •  (0)  (Fig.  56),  The  functions  and  •  T  are 

pep  e  f#t 

H 
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extremely  close.  They  have  the 

same  asymptote  ♦’  *  h/ 2r.  For  values 

of  h/t~ of  the  order  of  a  few  units  they 

are  practically  the  same. 

An  examination  of  these  curves 

furnishes  several  important  results. 

No  matter  what  the  frequency  of  analysis 

»»  is: 

0 

1” .  The  smaller  fc/r  is,  hence  the  more  points  one  uses  from 
the  interf erogram,  the  less  noise  there  is  in  the  spectrum: 

2 * .  For  values  of  h/r  above  a  few  units,  the  power  of  tha 
noise  is  proportional  to  the  interval  k; 

3* .  One  gains  practically  nothing  in  choosing  h/r  <  1  (father 
than  h/r  =  1),  i.e.,  in  choosing  more  than  one  point  per  time 
constant. 

— ^  These  results  can  be  stated  exactly  when  one  studies  the 
variation  of  the  s/n  ratio,  in  the  calculated  spectrum  based 
on  a  particular.. interferogram,  as  a  function  of  h/r. 

(2)  Variation  of  the  S/N  Ratio  as  a  Function  of  the  Interval 
k_ 

Noise  is  minimum  when  one  uses  a  continuous  sequence  of 
values  of  I(t);  it  is  then  characterized  by 

(V.10)  (c9x)m  -  2  Q  T  »  (0)  ♦„  r(0). 

Let  us  call  (c/n)^  the  maximal  value  that  the  s/n  ratio  then 
takes  in  the  calculated  spectrum: 

<»/n>H  '  *AV«  *  S/',J  «  1  V0)  Vi0)' '  • 

If 

The  s/n  ratio,  obtained  for  a  value  It  of  the  grating  interval, 
has  a  form 


Fig.  56. 
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»/n  -  s/ey  - 


-  s/viQT®p(0)*|;>T(Vr). 

The  quotient  S  *  (s/n)/(s/n)M 
measures  the  loss  in  the  s/n  ratio 
due  to  the  fact  that  one  usee  only 
discrete  values  of  I(t): 

S  -v/*™  (0)/*  TC*/r). 

o  •  , 

It  is  represented  by  a  curve  which;  in  principal,  is  valid  for  a 
single  value  of  the  product  But  being  given  tne  weak  varia¬ 

tions  of  *  witn  v ,  the  curves  S  will  be  practically  Identical. 
Figure  5  7  corresponds  to  *<0  ■  0.  One  sees  that  by  taking  the 
F.T.  by  using,  instead  of  a  continuous  sequence  of  values  of 
X(t),  discrete  values  taken  at  time  intervals  equal  to  the 
time  constant,  a  slight  weakening  of  the  s/n  ratio  in  tM 
spectrum  is  brought  about^it  is  multiplied  by  0.95. 


(3)  Experimental  Verifications 

These  have  been  made  based  on  a  synthetic  interferogram 
I*(t)  having  a  constant  component  I(t)  *  40  u  (arbitrary  unit), 
whose  F.T.  will  be  a  line  centered  on  the  frequency  0,  and  some 
fluctuations  of  mean  square  error  •  37.5  u  registered  with  a 
time  constant  r  ■  2  sec  and  represented  by  Fig.  37a.  The  total 
recording  time  is  T  ■  360  sec. 

The  F.T.  of  this  interferogram  made  by  using  a  continuous 
sequence  of. values  I(t)  is  composed  of  a  line  centered  on  the 
zero  frequency  and  some  fluctuations  whose  Man  square  error 
is,  according  to  (IV. 23), 

ox  «v/4  x  0.406  x  360  x  }  x  37.5  x  37.5  »  1,281  u  fee. 

The  F.T.  of  I*(t)  has  been  calculated  numerically  with  three 
different  values  of  ft. 
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gram  for  Different  Values  of  h .  __ 

1°.  1‘  -  1.2  sec.  Number  of  points  taken  from  the  interfero- 

gram,  N  =  300;  h/r  -  0.0.-  According  to  Fig.  57,  S  then  equals 
0.98  whence  the  predicted  value  of  the  mean  square  error  of 
the  fluctuations  in  the  calculated  spectrum:  (predicted)  * 

1,281/0.98  =  1,307  u  sec.  The  mean  square  error  of  the  fluc¬ 
tuations  measured  in  the  *  alculated  spectrum  (Fig.  58a)  by 
taking  the  square  root  of  the  mean  of  the  squares  of  the  ordi¬ 
nates  chosen  from  this  spectrum  every  10”  cps  (thus  three 
points  per  correlation  radius)  is 

<rx  (measured)  ~  1,284  u  sec. 

2* .  =  10.8  sc-c,  N  -  5  3  and  h/r  -  0.4,  whence  S  =  0.55:  _ 

ax  (predicted)  ~  1, 281/0.55  -  2,329  u  sec.  The  calculated 
spectrum  repeats  every  1/10.8  =  0.0924  cps  (Fig.  58b).  If  one 

i 

were  to  trace  the  spectrum  above  the  frequency  v  -  0.046  cps 
one  would  observe  a  symmetry  in  the  spectrum  with  respect  to 
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this  value.  To  measure  the  mean  square  error  of  the  fluctua¬ 
tions  we  have  thus  used  only  a  portion  of  the  noise  between 
v  *  0.004  cps  (starting  with  this  frequency,  the  line  centered 
on  0  has  a  negligible  height)  and  i>  *  0.046  cps: 

a x  (measured)  =  2,016  u  sec. 

3“ .  h  -  21.6  sec,  N  =  17  and  V./r  =  10.8.  whence  S  «  0.4. 

ffv  (predicted)  =  l,2Bl/0.40  =  3,200  tl  sec.  The  ealeu- 

A  h 

lated  spectrum  repeats  every  1/21.6  =  0.0462  cps  (Fig.  58c). 
The  curve  representing  it  is  therefore  symmetric  with  respect 
"ttr  the  ’frequency,  v  =•  0l.023iQ5  cps.  TJie  mean  square  error  of 
the  fluctuations  was  then  calculated  using  only  that  portion 
of  the  nol*se  taken  between  r  *  0.004  cps  and  v  m  0.023  cps. 

(measured)  *  2,392  u  sec. 

The  errors  can  be  explained  by  the  fact  that,  being 
given  the  frequency  of  the  repetition  of  the  spectra,  the 
mean  square  errors  have  been  calculated  from  noise  samples 
extending,  over  frequency  bands  much  too  narrow  to  hope  to  get 
a  correct  value  of  <r^.  If  the  noise  sample  were  sufficiently 
long,  the  mean  value  of  the  fluctuations  would  be  zero;  but, 
it  equals  -  24  u  sec  in  the  first  case,  +  133  u  sec  ill  the 
second  and  +  427  u  sec  in  the  third,  which  explains  why  tha 
error  is  larger  in  the  third  case. 

(4)  Problems  Posed  by  tha  Choice  of  f  and  by  Locating  the 
Zero  Point 

The  interf erogram  recorded  through  a  low-pass  filter  is 
the  convolution  of  the-  thioretical  intc-rf  erogram'  with  the  im¬ 
pulse  response  of  the  filter  which  is  necessarily  asymmetric. 
This  is  equivalent  to  saying-  that  the  filter  introduces  a  • 
phase  shift,  a  nonlinear  function  of  the  frequency  hence  of 
phase  distortion.  Let  us  choose  as  time  origin  that  instant 
where  the.  geometric  phase  difference  is  zero-.  If  the  phase 
shift  were  linear  one  would  have  _ 
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I(t)  =  \  B(i»)  cog  [2irt't  +  vK*')}  di>  » 

,;o 

with  *p(i>)  -*  ti'  ,  r  being  the  time  constant  of  the  filter,  whence 

r\  oo 

I(t)  =  \  B(i>)  cos  2jri>(t+r)  dr  »  I’ft+r), 

Jo 

i.c.,  the  interf crogram  would  be  recorded  with  a  certain  delay 
r,  but  without  any  deformation.  In  order  to  have  a  correct 
spectrum,  it  is  sufficient  to  choose  as  origin  for  the  analysis 
the  time  t  &  r. 

Since  with  ordinary  filters,  the  phase  shift  is  never 
linear,  it  is  no  longer  proper  to  speak  of  the  oero  point  of 
th>.  int»  rl crogram.  The  origin  can  bo  chosen  correctly  only 
for  single  frequency.  In  the  reproduced  spectrum,  the  image 
of  the  line  Corresponding  to  this  frequency  will  be  symmetric. 
The  images  of  tin.  lines  corresponding  to  all  the  other  fre-( 
quoncios  will  be  as  >  now.  trie,  the  asymmetry  for  a  given  fre¬ 
quency  being  greater  the  greater  the  time  constant  of  the 
filter  used.  We  arc  going  to  examine  how  it  is  necessary  to 
•  choose— the  time  constant  in  order  that,  from  one  extremity  to 
tla  other  of  the  spectrum , the  asymmetry  doesn^t  exceed  a  cer¬ 
tain  value  chosen  in  advance  and  how  one  determines  the  zero 
point  of  the  interf crogram.  Our  study  will  be  limited  to  the 
i  a. *  ol  Che  RC  filter,  the  most  frequently  used  with  very  low 
i:  vqu.  neies.  The  phase  displacement  *p(r)  introduced  by  the 
1  n  qUenry  i*  is  then 

.(.V .il)  *?{v)  ~  arc  tan  2*vt , 

wi.i-i.  ■.  .  I-,  s ponds  for  this  sum.  frequency  to  a  delay 

(V.  1 2 )  arc  tan  2i rvr. 

Let  u s  continue  to  v Loose  for  the  0  instant  the  one  where  the 
geometrical  phase  difference  is  zero.  One  can  locate  it,  for 
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Fig.  59.  Curve  a,  Interferogram  Re¬ 
corded  Through  an  RC  Filter;  Curve  b, 
the  Particular  Sine  Wave  Contained  in 
This  Interf erogram.  _ - 

example,  by  recording  at  the  same  time  the  interf erogram  to 
be  studied  and  the  one  originating  from  a  continuous  intense 
source,  the  latter  recorded  without  an  RC  filter. 

If  one  chooses  as  the  nero  point  of  the  interferogram 
the  one  which  corresponds  to  t  -  0,  the  delay  of  every  fre¬ 
quency  in  the  spectrum  will  be  given  by  the  formula  (V.12); 
only  the  apparatus-function  corresponding  to  the  0  frequency 
will  be  symmetric.  But  it  is  possible  to  make  a  better  choice 
of  the  origin  for  analysis  in  order  that  the  relative  phase 
shift  from  one  extremity  of  the  spectrum  to  the  ocher  may  be 

smaller.  * 

Let  8  be  the  instant  corresponding  to  the  point  that  we 
have  chosen  as  the  Hero  point  on  the  recorded  interferogram 
(Fig.  59,  curve  a).  With  this  new  origin  for  analyzing  the 
interferogram,  the  delay  for  a  frequency  v  present  in  the 
interferogram  (Fig.  59,  curve  b) ,  becomes  or  *  0 ,  which 

corresponds,  by  calling  y  the  period  of  the  phase-shifted 
signal, 

(V.13)  2  *  a/ if  *  arc  tan  2*vt  -  2 t 

a/ri  and  the  coefficient  c/l,  considered  in  Chapter  III,  which 
measures  the  error  made  in  choosing  the  zero  point  by  taking 
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as  unity  the  wave  length  of  the  signal,  are  equal.  In  what 
follows  we  replace  a/i\  by  c/x, 

(V.13b)  X  =  arc  tan  *tvT  - 

which  can  also  be  put  in  the  form 

(V.13c)  \  =  ~2ir  arc  tan  ^ irvT  -  'J~r  2ivt  . 

c/\  will  be  zero  for  v  =  0  and  for  a  frequency  which  will 

be  given  as  a  function  of  0  by  the  relation 

■  •  i 

(V.  14 1  vQ  *  2»0  atc  tan  2wot. - ^ 


Hence  in  the  spectrum  the  apparatus-functions  corresponding 
to  the  frequencies  0  and  will  be  symmetric.  Those  which 
correspond  to  the  other  frequencies  will  be  asymmetric.  ’ 

Since  the  expressions  (V.  13b)  and  (\£13c)  allow  us  to  calcu¬ 
late  the  coefficients  c/X  corresponding  to  them,  we  can  deter¬ 
mine  their  appearance  as  in  Fig.  30.  We  tried  to  represent 
the  variations^!  c/A  in  the  same  spectrum  as  a  function  of 
2*i>t  for  different  choices  of  the  zero  point,  i.e.,  for  dif¬ 
ferent  values  of  2vv  t. 

o 

In  Fig.  60  we  com  idci\-d  the  case  corresponding  to 

2nv  t  =  0,  2*i>  t  ~  1  arid  2m>  t  -  2.  Curve  d  represents  the 

0  0  0 

function  (l/2rr)  arc  tan  2ttvt.  The  lines  a,  b,  c,  correspond¬ 
ing  to  the  three  cases  considered  have  respectively  the  slopes 
0/2kt  =  0.159,  o/2vt  =  0.12  5  and  o/ 2jtt  *  0.099.  According  to 
(V.13c)  one  deduces  immediately  the  corresponding  values  of 
c/X  (curves  a',  b',  c’). 

An  i  vamination  of  the  figures  shows  that: 


1° .  For  0  -  0t  the  phase  shift  is  always  negative  (Fig.  60, 
curve  a1); 

2° .  For  the  other  values  of  0,  the  phase  shift  increases 
starting  with  0,  passing  through  a  maximum  (c/X)^  for  a 
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frequency  1/  ^  decreases  to  0  for  a  frequency  and  becomes 
negative.  For  a  frequency  ,  it  equals  -  (e/A',4.  Hence  if 
the  frequencies  contained  in  the  spectrum  to  be  analyzed  are 
all  taken  between  0  and  v* ,  c/\  for  these  frequencies  never 
exceeds  |(e/X)M!  and  all  the  corresponding  apparatus-functions 
will  have  a  coefficient  of  asymmetry  below  CKc/XJ^J  which  one 
can  measure  from  Fig.  23  or  2?. 

If  one  choost  c  v  =  1/2jtt  (Fig;  60,  curve  b’)  one  has 
P$  -  1.25/2  xt  and  C  =  0.06. 

For  v  *  l/xr  (Fig.  60,  curve  c’)»  C  *  0.09. 

From  Fig.  bu  wc  have  been  able  to  deduce  Fig.  61  in 
which  we  have  represented  as  a  function  of  the  maximum 
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Fig.  01.  Variations  of 
2irr  r  and  2m>  r  as  Func- 

M  0 

tions  of  the  Maximum  Co¬ 
efficient  c/A  Allowed  in 
the  Spectrum 


coefficient  c/A  that  one  allows 
in  the  reproduced  spectmm  the 
values  of  2kv  r  and  2*vut, 
being  the  limit  frequency  con¬ 
tained  in  the  spectrum  to  be 
studied. 

Curves  28  and  29  relating  C 
and  c/x,  you  may  recall,  were 
established  in  the  case  where  the 
weighting  function  is  A^ft).  It 
is  easy  to  establish  them  for  each 
type  of  apparatus- function.  Figure 
62  gives  directly  and 

as  functions  of  C  in  the  case  where 
one  uses  A#(t). 


Fig.  G2.  A  Plot  of  2w» '  r 

and  2n*  t  as  Functions  of 

o 

the  Maximum  Coefficient  of 
^Asymmetry  C..  Allowed  in  the 
Spectrum. 
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In  Conclusion,  when  one  does  not  want  the  relative  phase 
shift  from  one  extremity  to  the  other  of  the  spectrum  to  ex¬ 
ceed  a  certain  value,  two  problems  are  posed  simultaneously: 

The  choice  of  the  time  constant  and  that  of  the  zero  point  of 
the  interf erogram.  The  two  Known  quantities  are:  the  maxi¬ 
mum  frequency  contained  in  the  spectrum  to  be  studied  and 
the  maximum  coefficient  «/A  that  one  can  allow.  Figure  bl 
permits  one  to  determine  the  ' corresponding  value  of 
then  of  r  since  is  known;  the- same  figure  then  permits  one 
to  determine  from  the"  responding  value  of  2tv^r.  When. 

one  knows  v  ,  (V.14)  permits  calculation  of  the  corresponding 

o  _ 

value  of  0. 

This  is  the  rigorous  mothem  for  determining  the  time  ■- 
«  onstant  to  use  and  the  r.ero  point.  There  exists- an  approxi- 
:u  ife  method  to  resolve  the  latter  problem;  -it— consists  in 
f  ■  cording  at  the  same  time  as  the -intorf erogram  another  pro¬ 
duced  by  a  band  of  radiation  isolated  in  a  continuous  spectrum 
by  a  l alter,  whose  maximum  of  transmission  coincides  with  the 
spectral  region  to  be  studied,  and  recorded  with  the  same  time 
ro'irtant.  The  maximum  of  the  «  entral  fringe  of  the  correspond¬ 
in'!  inter] erogram  then  indicates  approximately  the  point  that 
one  ought  to  choose  as  the  origin  lor  analyzing  the  interfere** 
gram. 

* 

( )  dtudy  of  a  Particular  Case . 

Suppose  that  wc  study  a  spectrum  extending  from  0  to 
(Fig.  W) .  If  one  allows  C  =  0.0'^  (this  leads  to  an  asymmetry 
of  the  ord-.r  of  l/inn,  .  incc  C  0.0/;  for  a  p<  rlcctly  symmetri¬ 
cal  i  unction)  on.  d.  ne.  (Fig.  cl):  .  r  l/2jrrH  and 
v  -  0.5.1 ' 2rr  *  0.6-1 

llnd.-r  these  i -and  it  ions  the  noise  spectrum  is  such  that 
■r(j»N)  •-  i  »(0)  (Fig.  fo)  and  the  frequency  in  the  repro¬ 

duced  sp>.- cl  rum  is  reduced  fcy.  a  factor  i/\2*  If  one  mak  c  a  a 
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numerical  F.T. ,  by  choosing  the  minimum  interval  ft  for  which 
there  is  no  overlapping  of  the  optical  spectra,  hence 
ft  =  l/2  i'„,  the  quotient  ft/r  equals  »;  according  to  Pig.  57, 
f»  then  equals  0.7. 

Hence  when  one  studies  a  spectrum  extending  from  the  0 
frequency  to  a  maximum  frequency  »'M,  one  can,  in  picking  the 
minimum  number  n  of  points  in  the  interf erogram,  obtain  an 
a/n  ratio  which  equals  7/10  of  the  maximum  s/n  ratio  that  one 
could  obtain  with  the  given  interf erogram,  with  a  distortion 
of  the  apparatus-function  of  the  order  of  l/lOO,  providing  the 
constant  and  the  aero  point  of  the  interf erogram  are  properly 
chosen.  In  order  to  make  a  eompaiison  of  the  intensities  of 
the-  lines  in  the  calculated  cpoctrum,  it  would  evidently  be 
necessary  to  make  a  correction  taking  into  account  the  varia¬ 
tion  of  the  filter  gain  in  the  domain  of  the  frequencies 
under  consideration. 

Suppose  now  that  the  spectrum  of  the  Pouricr  frequencies 
is  limited  between  and  t»M  with,  for  example,  »m  •  4  *»^/5 
(lig.  63).  The  time  constant  to  be  used  remains  thr  same  as 
K  1  m\  if  wo  fix  the  same  tolerances  and,  in  order  to  obtain 
tin  Tame  s  n  ratio  from  an  inti  rfex  i  ^iam  recorded  during  the 
xidiiu  time  T,  the  intei’val  ft  ought  always  to  be  equal  to  l/2  Pjj. 
But  according  to  Chapter  II,  in  order  that  there  be  no  over¬ 
lapping  of  the  optical  spectra  it  would  have  been  sufficient 
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to  choose  an  interval  h  =  l/(2An)  =  5/2  The  number  H  of 

the  points  to  choose  in  order  to  obtain  the  signal/noise  ratio 
s/n  =0.7  (s/n)M  is  equal  to  five  times  the  minimum  number  n 
of  points  that  it  would  be  necessary  to  choose  in  order  that 
the  different  optical  spectra  calculated  do  not  inf ringe  upon 
each  other.  The  duration  of  the  calculation  hence  is  multi¬ 
plied  by  5.  We  shall  study  in  what  follows  different  processes 
permitting  the  reduction  of  this  calculating  time,  all  the  while 
keeping  this  same  s/n  ratio  in  the  calculated  spectrum. 

2.  STUDY  OF  DIFFERENT  METHODS  FOR  OBTAINING  A  MAXIMUM 

S/N  RATIO  WITH  A  MINIMUM  CALCULATION  TIME  FROM  AN  ‘  ' 
INTERFEROGRAM  RECORDED  DURING  A  TIME  T  ’ 

We  have  seen  in  Chapter  IV  that  thi  s/n  ratio  in  the 
spectrum  varies  as  the  square  root  of  the  duration  of  measure¬ 
ment.  Suppose  that  we  have  at  our  disposal  a  certain  time  to 
make  the  measurement  and  that  we  have  chosen  then  to  take  the 
numerical  F.T.  of  the  interf erogram  obtained.  We  have  seen 
in  the  preceding  paragraph  that  the  numerical  F.T.  of  an  inter- 
f erogram  recorded  through  a  low- pass  filter  always  gives  an 
s/n  ratio  below  that  which  one  could  obtain  by  using  the  en¬ 
tire  interf erogram.  One  way  of  illustrating  this  fact  con¬ 
sists  in  saying  that  several  peaks  of  the  square  of  the 
apparatus-function  explore  the  noise  spectrum  when  one  has 
chosen  the  minimum  number  n  of  the  points  oven  though  a  single 
peak  of  the  apparatus-function  explores  the  optical  spectrum. 

In  order  to  obtain  the  maximum  s/n  ratio  from  an  inter- 
f erogram  recorded  during  time  T,  at  the  same  time  using  the 
greatest  interval  r  possible  to  make  the  numerical  F.T.,  it 
is  necessary  that  the  noise  spectrum  be  ::tro  and  outside  the 
domain  oi  the  Fourier  frequencies.  Belov/  we  are  going  to 
consider  tins  and  to  d-. scribe  in  detail  four  different  methods 
for  obtaining  this  result. 
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a.  Use  of  a  Band-Pass  Filter.  The  low-pass  filter  which 

is  used  for  recording  the  interf erogram  is  replaced  by  a  band¬ 
pass  filter.  This  method,  contrary  to  the  three  following 
ones,  ha's  never  been  used  in  any  practical  way.'  j._ _ 

b.  Change  of  Frequency  by  Heterodyne  Before  Recording. 

This  method  which  has  been  proposed  by  Mertz  [54)  under  the 
name  of  ’’Heterodyne  Spectrometry’’  is  equivalent  to  the  pre¬ 
ceding,  but  in  fact  is  more  elegant  and  easier  to  put  into 
practice.  — - - 

c.  Static  Recording.  In  the  absence  of  noise,  all  the 
pieces  of  information  concerning  the  spectrum  are  contained 
in  a  number  of  points  n  in  the  interferogram.  One  easily 
concedes  that,  in  the  presence  of  noise,  one  obtains  the 
maximum  s/n  ratio  in  the  spectrum  if  one  uses  the  total  time 
T  fop  the -measurement  of  these  points.  Hence  the  method  con¬ 
sists  in  stopping  the  moving  mirror  of  the  interferometer  for 
each  of  the  n  interesting  values  of  the  phase  difference  in 
integrating  the  signal  during  a  time  T/n  for  each  one  of  these 
points  and  in  displacing  the  moving  mirror  very  quickly  with¬ 
out  recording  between  these  different  values  of  6.  The  mechan¬ 
ical  realization  of  this  type  of  displacement  being  somewhat 
difficult,  we  have  used  a  variant  of  this  method,  registering 
at  variable  speed:  rapid  displacement  of  the  mirror,  then 
slowing  up  in  the  neighborhood  of  the  useful  points. 

d.  Mathematical  Filtering.  In  this  last  method,  the  re¬ 
cording  is  done  with  a  simple  low-pass  filter  and  the  number 

of  points  N  to  choose  remains  determined  by  the  considerations 
developed  in  the  preceding  paragraph  (l:  is  of  the  order  of  r). 
But  before  proceeding  with  the  calculation  of  the  spectrum, 
one  makes  an  intermediate  calculation:  that  of  the  convolu¬ 
tion  of  the  recorded  interf erogram  I’ (6)  with  the  impulse 
response  of  an  ideal  filter  whose  band  pass  would  have  as  its 
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-width  that  of  the  spectral  interval  occupied  by rthe  Fourier 
frequencies.  The  result  is  a  new  interf erogram  Iw(d)  identi¬ 
cal  to  that  which  would  have  been  obtained  by  directly  filter¬ 
ing  the  signal  by  a  perfect  band-pass  filter.  It  can  then  be 
completely  represented  by  a  number  of  points  n  reduced  with 
respect  to  N  in  the  same  ratio  as  the  band  passes.  One  then 
makes  the  F.T.  of  I”(6 )  and  the  calculation  time  is  reduced 
very  nearly  in  the  same  ratio  n/N  (the  calculation  time  of 
the  convolution  is  very  short}  which  constitutes  the  clear 
advantage  of  the  method. 

(1)  Electric  Filtering 

The  filters  should  have  a  band  pass  centered  about  the 
mean  frequency  of  the  useful  spectrum.  The  phase-shift  con¬ 
siderations  introduced  by  this  filter  are  like  those  intro¬ 
duced  by  a  low-pass  filter.  In  this  case  one  will  find,  for 
a  suitable  choice  of  the  origin  for  the  analysis  of  the  inter- 
ferogram  recorded  between  0  and  L  and  a  simple  cosine  transform, 
a  single  frequency  for  which  the  image  of  the  line  will  be 
perfectly  symmetric.  We  have  not  studied  this  process  in 
detail  since  at  the  very  low  frequencies  that  we  use,  it  is 
difficult  to  make .band-pass  filters  and  the  same  result  can 
be  obtained  more  elegantly  by  the  process  of  changing  the 
frequency. 

—  -  Let  us  remark  that  the  process  of , electrical  filtering 

is  interesting  only  if  the  width  of  the  interval  Ap  of  the 

spectrum  of  the  Fom  ior  frequencies  and  the  frequency  »  ^ 

satisfy  the  relation  v  _  >  2Aj*.  In  fact,  wc  have  shown 

max 

(Fig.  18)  that  a.  long  as  i’.ua<  stays  below  2 An,  the  interval 
k  chosen  iO  that  th  iv  is  no  overlapping  of  the  calculated 
spectra  has  the  .same  width  V.  *•  1/2M  if  the  Fourier  frequencies 
extend  from  0  to  ■ 
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(2)  Method  of  Changincr  the  Frequency 

This  method  has  been  proposed  by  Mertz.  We  do  not  dis¬ 
cuss  it  here,  limiting  ourselves  with  recalling  very  briefly 
its  principle  and  giving  an  example  of  its  use  in  Chapter  VI. 

In  lieu  of  recording  I(t)  as  before,  one  records  simul¬ 
taneously  I'(t)  and  I"(t)  obtained— as- the  product  -of  I(t)  with11 
two  sinusoidal  signals  of  frequency  respectively  in  "phase  and 
90°  out  of  phase  with  I(t).  The  sum  of  the  cosine  transform  _ 
of  I’(t)  and  of  the  sine  transform  of  IM(t)  gives  the  spectrum 
wc  are  looking  for.  A  delicate  problem  which  presents  itself 
is  the  adjustment  of  the  phases  of  the  signals. 

If  the  phases  are  correctly  adjusted,  the  method  of 
changing  the  frequency  gives,  for  an  appropriate  choice  of  * 
in  the  interior  of  the  domain  ,  a  s/n  ratio  equal  to 

that  obtained  with  a  band-pass  filter  and  with  a  calculation 
time  divided  by  2.  But  it  is  necessary  to  make  the  simul¬ 
taneous  recording  of  the  two  interferograms  which  presents  a 
certain  complication. 

One  can  be  satisfied  with  recording,  a  single  interfero- 
gram  if  one  chooses  v  outside  the  domain  v v #  or  if  the 
signal  of  frequency  is  out  of  phase  by  r/4  with  respect  to 
the  Fourier  frequencies  contained  in  the  interferogram  before 
changing  the  frequency;  but  in  the  two  cases  the  calculation 
time  is  at  least  multiplied  by  2  and  the  s/n  ratio  in  the 
calculated  spectrum  is  reduced  by  the  factor y/2/2. 

If  the  signals  are  out  of  phase  by  an  angle  <p  *  2wvQ€, 
in  the  case  where  is  outside  v  ,  v  ,  the  s/n  ratio  in  the 
calculated  spectrum  by  taking  a  simple  cosine  transform  of 
the  interferogram  recorded  for  o  varying  between  -  L  and  +  L 
is  weakened  by  a  factor  cos  (pt  the  calculation  time  being 
double  that  which  one  would  have  if  the  signals  were  in  phase. 
If  v  is  inside  v in  order  to  obtain  a  correct  spectrum, 
the  calculation  time  is  multiplied  by  8  with  respect  to  the 
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case  where  the  signals  are  in  phase  and  one  can  make  the  same 
remarks  about  the  lack  of  linearity  in  the  procedure  as  were 
made  in  paragraph  IV.3(2)A.  In  the  case  where  <p  is  exactly 
equal  to  r/4,  the  calculation  time  will  simply  be  multiplied 
by  4,  but  the  s/n  ratio  will  be  reduced  by  the  factor y/2/2. 

( 3 )  Static  Recording 

Let  us  recall  that  the  interf erogram  I'(t)  recorded 
without  a  filter  when  the  phase  difference  5  varies  linearly 
as  a  function  of  time  is  the  superposition  of  the  interfero- 
gram  I(t)  that  one  would  have  in  the  absence  of  noise  and  the 
fluctuations  x(t)  for  which  one  may  assume  a  constant  spectrum 
{noise  density  ®(0)J. 

Let  n  be  the  minimum  number  of  points  to  choose  from  the 
interf erogram.  The  total  measurement  time  T  can  be  used  either 
to  record  I’(t)  as  has  been  considered  earlier,  or  to  determine 
uniquely  the  indispensable  n  values  needed  to  make  a  numerical 
F.T.,  each  being  measured  during  .the  time  T/n. 

In  the  second  case  the  operation  is  performed  in  the 
following  way:  the  interferometer  is  adjusted  to  the  zero  phase 
difference.  If  there  were  no  noise  the  signal  I'(0)  would  be 
constant  and  would  keep  the  value  1(0).  But  in  the  presence  of 
noise,  the  signal  I'(0),  when  the  phase  difference  keeps  the 
constant  value  6  —  0,  fluctuates  with  time.  The  mean  value 
which  we^shall  call  I"(0)  taken  for  a  time  T*  *  T/n  is 

r*<o)«^r^T  i'(o>  dt  . 

Let  Is1  be  the  phase  difference  corresponding  to  the  second 
point  to  be  measured.  Suppose  that  the  passage  from  6  =  0  to 
6  -  h1  can  be  made  instantaneously.  At  the  instant  T’,  then 
6  equals  ft*  and  during  a  time  T*  we  take  the  mean  of  the  rlgnal; 
the  number  found  is  I" (ft').  In  the  same  way  one  determines 
I"(2  *»),  ...»  I*(n  *»)♦  etc.  so  that  one  can  write  I*(0), 
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Fig.  64.  (a)  Square-Pulse  Function  D(t); 

(b)  Autocorrelation  Function  of  D(t). 


I” (T 1 ) »  . ..,  I” l (n-l)T* ] .  Hence  one  obtains  at  the  end  of 
the  measurement,  n  equidistant  discrete  values  of  a  new 
interf erogram  I"(t)  which  is  the  superposition  of  _I(t)  and 
the  fluctuations  x'(t)  for  which  we  are  going  to  determine 
the  spectrum  with  respect  to  that  of  x(t).  Vinokur  [55] 
studied  a  particular  case,  of  this  problem:  he  showed  that 
under  the  double  condition  that  the  signal  y(t)  be  periodic 
and  that  the  period  T  is  known  precisely,  the  process  called 
"summation"  which  consists  in  adding  up  the  successive  elonga¬ 
tions  v(t4kT)  leads  to  an  improvement  in  the  signal/noise 
ratio.  We  can  give-  a  quick  proof,  valid  in  the  general  case, 
showing  that  the  spectrum  ©'(iO  is  of  the  form  ®’(v)  «  ®(0) |G(p) I* , 
G(»0  being  the  filter  gain  produced  by  the  integration. 

The  noise  x'(t)  has  the  form 


(V. 15  )  x’(t) 


x(t)  D(t)  dt 


9 


D(t)  being  a  square  pulse  function  which  equals  l/T'  for  t 
taken  between  0  and  T’  and  is  zero  outside  of  this  interval 


(Fig.  64a). 

The  mean  square  error  of  x’(t)  can  be  put  in  the  form 


■  —i.  - ......  ,  •  » /■*  «e 

(V.  16)  p'(0)  =  x,8(t)  *  \  „  \  xTcTxTtT  D(t)D(t')  dt  dt 1 

*V  J0 


The  term  x(t}x(t 7 )  is  the  autocorrelation  function  C  (#i)  of 


x(t).  The  term  £  D(t)D(t')  dtT,  that  one  can  also  write 
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Fig.  65.  Noise  Spectrup.  After  Inte¬ 
grating  the  Signal.  ''  '  • ':~~ 


1/T’,  the 
D(t)s 


D(t-p>  dt  by  putting- 1’  «  t-pis,  within  a  factor  of 
autocorrelation  function  of  the  square-ptilaN "function 


T’ 

D(t)  D(t-p)  dt  •  T'C^Cn). 
(V.16)  can  then  be  written 


(V.17)  p*(0>»T»£  cx0»)  CD0»)  <*» 

From  Parseval's  theorem, 

(V.18)  p«( 0)  -  T*  J  ^os^x001  TcosICD(,l)1  d*v 

The  cosine  transform  of  the  autocorrelation  function  of  x(t) 
is,  by  definition,  the  power  spectrum  of  x(t),  i.e.,  »(0) 

(Fig.  6S,  curve  a].  The  autocorrelation  function  Cg(p)  is  a 
triangle  function  of  base  2  T',  of  height  l/T,a  (Fig.  64b). 

Its  cosine  transform  is  the  function  (l/T')  (sin*  »rT,/(»»»  T*)*l 
whence  * 

c(o)  -  V"  .(to  d, ,  ■ 

which  can  also  be  written 

CV.W)  $_r  .,<►)  -  f  ~  V0>  (^r^y  *'• 
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Fig.  66.  Noise  Spectrum  Explored  by 
the  Square  of  the"  Apparatus-Function.:  . 

The  integration  therefore  has  had  the  effect  of  filtering  the 
noise  by  a  filter  of  gain  G(*>)  such  that 

(V.20)  iGCOl*  -  T'  (**)  (pig.  65,  curve  b) . 

(tv  T  ’) 

The  noise  in  the  spectrum  will  be  the  F.T.  X(v)  of  x’(t)  taken 
by  choosing  k  ■  T*.  The  mean  square  error  (ov)*  will  then.be 

'  *  JL  • 

the  chopped  energy  in  the  noise  spectrum  times  the  square  "of v 
the  apparatus-functions  centered  about  v  *  m  (l/T’),  The  re¬ 
sult  of  calculating  (ax)*  ia  obvious »  -one— of  the  maxima  is 
centered  about  the  0  frequency,  the  others  about  a  m  (l/T1), 
the  values  of  the  frequency  for  which  *p(p)  is  taro.  From 
this  one  concludes  Immediately 

**<#)  * 2  Q  T  V0)  {Fig*  66>  * 

an  expression  which  is  identical  to  v IV. 22b). 


<l,>The  different  "lobes"  which  are  l/T  apart  in  the  filter 
introduced  by  Vinokur  in  the  case  of  the  summation  that  he  uses, 
come  from  the  fact  that  ho  dots  not  make  a  true  integration,  but 
a  sum  of  tho  discrete  value-  of  the  signal.  These  different  lobes 
are  close  to  the  different  peaks  of  the  apparatus-function  that 
one  obtains  in  the  case  where  the  F.T.  is  made  numerically  utlng 
discrete  values  of  the  interferogram. 

ns 


t 


t 


Fig.  67.  Grating 
Function  R,/x' 

‘  Centered  on  the 
Frequency 


Fig.  C8 .  Product  of  *p(iO 

by  R|/ji • 


The  signal/noise  ratio  of  the  calculated  spectrum  is 
therefore  the  same  as  when  the  inte.rf crogram  I,(t>1»  recorded 
in  a  continuous  fashion  during  time  T  and  when  thfe  F.T.  is 

made  using  all  the  value,  of  X'(t). 

It  remains  to  show  that  if  the  F.T.  is  made  for  a  freq¬ 
uency  u  different  from  0,  the  mean  square  sir or  of  the  flue- 
— tuitions  "preserve  s  the  same  value.  _ . . 


X«m  ) 


-u  yj 


sin -—»»  T* 
(wt>  T * )* 


dv, 


R  /<f,  being  a  Dirac  distribution  of  periodic  support  of  interval 
1/T’  centered  about  the  frequency  »»#  (Fig.  67).  If  put  . 


(»v  T’) 


,y»  R»/T' 


(Fig.  68)  , 


and  Q(t)  *  T^lPC^)]  h*ve 


(V. 21) 


Q(0) 


The  F.T.  of  POO  (Fig.  GO)  is  the  convolution  of  a  triangular 
function  of  height  1/T’  and  of  base  2  T*  with  the  function 
R»  •*  T*  RJI,  cos  2*><#t,  R^i >  being  a  Dirac  distribution  of 
periodic  support  of  interval  T*.  One  sees  im**«Hat*ly,  ac¬ 
cording  to  Fig.  63,  that  Q(0)  -  1  whenea 


216 


Fig.  69.  Determination  of  Q(0). 


ffXi  V  “  ^Xco)  * 

This  method  of  recording  is  difficult  enough  to  use 
even  given  the  requisite  precision  to  vary  the  phase  differ¬ 
ence.  Nevertheless  one  can  consider  it  for  the  far  infrared 
where  the  tolerances  are  always  larger.  In  the  near  infrared 
we  have  used  a  variant  of  this  method,  variable  speed  recording . 

The  spectral  domain  occupied  bo  *  369  cm’1  is  such  that  the 
interval  ft  chosen  in  order  that  there  be  no  overlapping  of  the 
Fouricf  frequencies  is  ~  "V 

1/(2  Ae)  *  13.519  X  10"4  cm’1  -  21  A#  , 

A  being  the  wave  length  of  the  reference  line,  here  the  red 

0 

cadmium  line . 

The  phase  differences  are  referenced  with  the  aid  of  pre¬ 
selection  counters  and  adding  counters  whose  indicators  increase 
1  unit  each  time  that  the  phase  difference  increases  by  A#.  Thp 
carriage  is  displaced  alternately  at  a  slow  speed  during 
which  6  varies  by  3  A.  and  at  a  speed  V.  16  times  larger  during 
a  variation  A  6  *■  18  A.,  the  changes  in  speed  being  triggered  by 
the  preset  counters.  The  operation  is  done  in  the  following 
way  (Fig .  i'J).  The  scro  point  is  selected  by  a  white  fringe 
while  the  interferometer  is  moving  at  a  speed  V^;  when 
6  »*  22.5  A#  a  preset  counter  set  on  18  triggers  the  speed  V^; 
when  the  counter  has. returned  to  0,  it  triggers  speed  and 
at  tha  same  time  a  preset  counter  set  on  3,  etc.  The  ordinates 
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are  picked  from  the  interf erogram  for  8  *  0,  21  42  X#. 

. ..,  n.  21  A#,  as  indicated  in  Fig.  87.  The  method  has  worked 
correctly,  but  wc  have  abandoned  it  because  the  signal/noise 
ratio  in  the  calculated  spectrum  is  insufficient.  The  reason 
is  that,  under  the  described  conditions,  the  ratio  of  the  use¬ 
ful  measurement  time  (useful  time  that  one  can  write  r  l/2  A#, 
r  being  the  time  constant  used)  to  the  total  time  T  is  too  low. 
There  are  two  reasons: 

,1*T  Being  given  the  quality  of  the  mechanism  used,  the  con¬ 
stant  rate  of  speed  is  not  established  instantaneously;  three 
fringes  must  file  by  at  and  2/6  of  the  course  is  run  at  a 
lower  speed; 

y 

2* .  The  preset  counters  used  cannot  count  more  than  25  im¬ 
pulses  a  second,  and  the  ratio  of  the  two  Speeds  eih'be  only 
16.  The  elapsed  time  between  the  measured  points  hence  is  too 
large  a  fraction  of  the  total  measurement  time. 

Numerical  Filtering. 

Since  this  question  is  the  subject  of  a  more  detailed 
publication  [56j,  here  we  shall  give  only  a  succinct  resume. 

If  one  takes  the  convolution  of  the  interf erogram  I’ (6), 
recorded  as  above  through  a  simple  low-pass  filter,  with  the 
impulse  response  3(fi)  of  an  ideal  filter  G(a)  wnose  band  pass 
coincides  with  the  optical  spectrum  to  be  studied,  one  obtains 
an  interferogram  I"(6)  whose  F.T.  coincides  with  that  of  I' (6) 
in  the  domain  of  the  filter  and  is  zero  outside.  In  order  to 
have  the  maximum  s/n  ratio  in  the  calculated  spectrum  it  is 
then  sufficient  to  pick  oft  from  I"(6)  the  minimum  number  of 
points  compatible  with  the  width  of  the  optical  spectrum  and 
its  position  with  respect  to  the  origin  (see  paragraph  II. 2(3)B) . 

Suppose,  for  example,  that  the  relative  positions  of  the 
noise  spectrum  and  of  the  optical  spectrum  are  represented  by 
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Pig.  70.  The  -impulse  response  3(6) 
(Fig.  71a)  corresponding  to  G(</)  (Fig. 
71b)  has  the  form 


Fig.  70.  Noise  Spec¬ 
trum  and  Optical 

Spectrum. 


3(6) 
with  a 


TeoslG('» 


sintberfi 

t£06 


cos  2»«r  6, 


m 


(»,-Hra)/2. 


If  the  functions 

3(6)  and  I' (6)  had  a  band-limited  spec¬ 
trum  it  would  be  possible  to  calculate  their  convolutions  ex¬ 
actly  using  discrete  values  of  these  fujtcdohsr  Since  their 


spectra  are  unbounded,  from  discrete  values  of  these  functions 
one  can  calculate  only  a  value  $(6)  near  I"(6),  the  approxima¬ 
tion  improving  as  the  interval  '  chosen  for  taking  the  convolu¬ 
tion  is  the  smaller.  One  shows  that  the  s/n  ratio  obtained  by 
taking  the  F.T.  of  9(6)  with  the  maximum  interval  ft  (see  para¬ 
graph,  II. 2(3l  B)  is  exactly  the  same  as  with  the  spectrum  obtained 
by  taking  the  r.T.  of  I’(6)  with  the  interval  k-1.. -Since ~thcr 
convolution  is  calculated  numerically  much  faster  than  is  a  F.T. , 
the  calculation  time  gain  for  the  same  final  s/n  ratio  is  of  the 
on  r  h  '/h ,  i.e.,  of  the  order  of  Vie  shall  se*  an 

application  of  this  method  in  Chapter  VI. 


Fig.  71.  (a)  Impulse  Response  3(&)  of  the 

Filter;  (b)  Theoretical  Cain  of  the  Numeri¬ 
cal  Filter. 
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3.  CONCLUSION 

We  have  just  seen  several  methods  enabling  one  to  obtain 
for  a  minimum  calculation  time  the  maximum  s/n  ratio  compati¬ 
ble  with  the  measurement  time  T:  use  of  a  band-pass  filter, 
method  of  changing  frequency,  numerical  filtering.  The  con¬ 
clusion  of  this  study  is  the  following:  Each  time  that  the 
relative  width  of  the  spectral  domain  to  be  studied  is 

small,  it  is  advantageous  to  record  the  interferogram  by  the 
method  of  changing  the  fro quencyj-'-'t-he  -frequency  being 
chosen  close  to  the  middle  of  the  occupied  interval  Av; 
nevertheless  this  method  presents  a  certain  instrumental 
complexity  and  it  is  sometimes  possible  to  do  without  it  and 
record  the  interferogram  through  a  simple  low-pass  filter. 
Under  these  conditions  it  is  advantageous  to  make  a  numerical 
filtering.  This  last  method  can  be  Tiseful,  even  if  the  rela¬ 
tive  width  of  the  spectral  domain  to  be  studied  is  large, 
-provided  that  one  is  interested  only  in  one  portion  of  this 
domain. 
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VI.  EXPERIMENTS  AND  RESULTS 

With  the  aid  of  the  Fourier  transform  method  we  treated 
two  problems  in  the  near  infrared:  the  study  of  the  emission 
of  the  night  sky  near  1.6  (i*  the  region  where  the  atmosphere 
is  transparent,  and  the  study  of  the  re-Bmission  light  of 
germanium  about  1.7  jt .  These  two  sources  had  already  been 
studied  with  grating  spectrometers  [57,58] .  But  their  low 
energy  level  permitted  obtaining  in  a  reasonable  measurement 
time  only  low  resolutions  of  respective  orders  of  15CT  and 
1,500,  and  with  a  signal/noise  ratio  which  was  not  sufficient 
to  makp  evident  the  feeble  components  whose  existence  was  in 
question  [59]  or  to  determine  the  ldw  of  displacement  of  a 
line  as  a  function  of  the  magnetic  field.  The  end  pursued  is 
thus  an  increase  in  resolution  and  an  amelioration  of  the 
signal/noisc  ratio.  For  the  two  studies  we  used  the  same 
Michelson  interferometer  that  we  arc  going  to  describe  sum* 
marlly. 

1.  DESCRIPTION  OF  THE  INTERFEROMETER- DETECTOR  ENSEMBLE 

The  interferometer  used  is  completely  classical  (Fig.  72) 

OPTICS.  Everything  optical  is  glass;  hence  the  instru¬ 
ment  can  be  used  only  in  the  near  infrared  up  to  about  2.7  f< , 
that  region  where  glass  becomes  absorbent.  The  two  mirrors- 
are  7  in.  in  diameter  and  arc  aluminized;  the  beam  splitter  and 
compensating  plate  measure  7  cm  x  8  cm;  they  have  the  same 
thickness,  about  1  fi .  They  are  held  by  the  same  support,  the 
face  of  the  team  splitter  which  is  half -aluminized  being  a- 
gainst  the  compensating  plate.  There  is  no  ruling,  either  on 
the  beam  splitter  or  on  the  compensating  plate.  We  have  a- 
dopted  this  arrangement  in  order  to  obtain  easily  a  good 
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Fig.  72. 


parallelism  between  the  two  plates,  which  is  necessary  in 
order  that  the  apparatus  be  perfectly  symmetric  and  that  the 
recorded  interferograms  have  a  zero  point.  But  the  solution 
is  not  perfect  (because  of  the  difference  of  the  dispersions 
and  of  the  phase  shift  of  the  reflection  on  the  semi-aluminized 
glass  side  and  the  air  side).  In  order  that  it  be  perfect,  it 
would  be  necessary  that  the  two  plates  be  glued  together  with 
a  substance  having  the  same  index  as  glass  which  would  hold 
..the  plate  together  although  it  is  vary  thin. 
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Fig.  73.  (a)  Diagram  of  the  Detector  Mounting; 

(b)  Optical  Diagram:  ,  Mirror;  L’ ,  Exit  Leris; 

L  ^^ield^Lens;  C,  Condenser;  R,  Detector;. 


DETECTOR.  The  detector  is  a  Kodak  cell  of  lead  sulfide.  ~~ 
The  noise  power  of  these  cells  increases  nearly  proportionally 
to  their  surface;  hence  it  is  advantageous  to  choose  a  cell 
having  the  smallest  surface  possible.  If  one  calls  S  the  area 
of  the  cross  section  of  a  beam,  R  the  aperture  of  the  exit  beam, 
2  the  area  of  the  cell  and  a>  the  solid  angle  of  the  beam  that  is 
received,  in  order  to  use  all  the  light,  it  is  necessary  that 
SR  =  so.  Several  arrangements  have  been  proposed  to  allow  the 
largest  possible  u>.  Williamson  reco.nmended  the  use  <j>f  a  "cone 
channel  condenser"  { GO  J ,  P.  Conncs  that  of  a  microscopic  con¬ 
denser  on  the  rear  face  of  which  is  fastened  the  cell  14].  We 
adopted  this  last  solution  (Fig.  73a).  The  numerical  aperture 

of  th.-*  ^ondenswr  u.  <!  is  rvsin  0  ,  =  1.4.  The  solid  angle 

max 

finally  accepted  by  the  cell  would  correspond  in  air  to  the 
value  of  w  *  1.45  x  2?rsd.  A  field  lens  makes  the  image  of  the 
mirrors  on  the  front  face  of  the  condenser  and  the  condenser 
makes  the  image  of  the  exit  hole  on  the  lens  (Fig.  73b)'.  The 
cell  is  cooled  to  the  tei«#perature  of  dry  ice. 


125 


NAVWEPS  REPORT  8099 


Fig.  74.  Cchenatic  Diagram:  A,  Miehelson  Inter-  \ 

ferometer;  A1,  Movable  Mirror  With  Cear  Drive;  B,  - 

,  Entrance  and  Exit  Holes;  C,  Lead  Sulfide  Detector;  ^  y 
D,  Photon  Multiplier;  E,  Cadmium  Lamp : - — ^  /  ■>' 

-f 

In  the  interferometer  the  reference  beim  follows  exactly  ^ 
the  same  path  as  the  light  beams  to  be  studied  (Fig.  74);  “it 
is  received  by  a  photomultiplier,  because  we  have  always  used 
as  a  reference  line,  a  line  taken  in  the  visible  part  of  the 
spectrum. 

MECHANICAL  PART.  The  moving  mirror  is  locked  on  a  car¬ 
riage  which  travels  on  parallel  slides.  The  quality  of  the 
slides  is  such  that  the  parallelism  is  preserved  to  about 
l/lO  of  a  fringe  of  a  red  cadmium  line,  for  an  8mm  transla¬ 
tion  of  the  moving  mirror  (hence  a  variation  of  the  phase 
difference  of  16mm)  with  a  5.  cm  usable  diameter  of  the 
mirrors. 

The  adjustments  for  parallelism  act  on  the  supports  for 
the  mirror.  The  gross  adjustments  are  made  following  a  clas¬ 
sical  procedure  of  working  directly  on  the  back  of_the  moving 
mirror  with  the  aid  of  three  leveling  screws  placed  120*  apart 
and  which  one  draws  out  once  the  screw  adjustment  is  made. 
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'  a 

* - 1 


(b) 

-ig.  / 5 .  (a)  Tine  Adjust¬ 

ment  Median  ism  of  the  Fixed 
.  Mirror,  (b)  Cross-Section  aa*. 


The  fine  adjustments  are  located  01.  the  fixed  mirror.  The 
mirror  is  mounted-on  a  flexible  part  which  can  be  oriented 
by  bending  and  which  has  the  form  shown  in  fig.  75a.  The  j 
cylindrical  part  B  has  an  8mm  diameter  and  is  10mm  long.  j 
Part  A  has  two  stool  aims  with  perpendicular  springs  and  at  ' 
45*  to  the  vertical  on  the  extremity  of  which  one  can  operate 
a  screw  with  a  fine  pitch  of  0.5mm  (Fig.  ,’jb).  The  torsion  ol 
B  which  results  therefrom  two  completely  independent 

adjustments.  One  complete  turn  of  the  screw  corresponds  to 
a  misad justment  of  3  half  fringe.  It  is  necessary  that  the 
ensemble  be  very  rigid  since  the  optical  beams  arc  relatively 
long  (about  15  cm).  Also  the  two  consoles  bearing  the  fixed 
mirror  and  the  sc pai  ating-c 'ompensat.r.g  ensemble  arc  locked  in 
a  heavy  cement  Support.  The  support  lor  thi.  two  slides  is 
locked  in  the  same  block  of  cement. 
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DRIVE  MECHANISM.  The  carriage  is  driven  by  the  aid  of  a 
micrometric  screw  with  0.5mm  pitch  of  any  quality.  It  is  use¬ 
less  to  use  a  screw  of  superior  quality  since  we  have  seen 
that  even  with  a  very  good  screw  the  phase  differences  have 
to  be  marked  with  the  aid  of  a  reference  line,  the  screw  is 
displaced  by  the  action  of  a  synchronous  motor  turning  at 
750  rpm  through  the  use  of  a  gear  box  and  _of  a  tangent  screw 

i 

giving  a  reduction  ratio  of  100.  The  gear  box  permits  one  to 
use  reduction  varying  by  powers  of  2  between  1  and  4,096,  in 
such  a  way  that  the  displacement  speed  of  the  moving  mirror 
varies  between  6.2*5  *  10  and  1.5  -X'  10  mm/ sec.  The -appara¬ 
tus  is  extremely  sensitive  to  vibrations  which  affect  the 
length  of  the  interfering  beams;  there  are  two  principal 
sources  of  vibrations:  those  coming  from  the  earth  and  those 
which  are  duo  to  th-  motor  and  all  of  the  components ‘ of  reduc¬ 
tion  and  transmission.  In  order  to  supipress  the  first  kind 
as  much  as  possible  the  block  of  cement  is  mounted  upon  rubber 
shock  absorbers, and  in  order  to  attenuate  the  second  kind,  the 
drive  axis  and  all  the  transmission  components  are  linked  to¬ 
gether  by  flexible  joints.  Of  course,  it  is  necessary  to  avoid 
temperature  variations.  Hcncc  the  interf erorfieter  is  placed  in 
the  interior  of  a  box  lined  with  insulating  material  and  covered 
with  aluminum  foil.  This  thermo-isolation  proved  to  be  suffi¬ 
cient  even  when  one  is  working  in  a  room  without  temperature 
control.  In  order  to  diminish  the  absorbtions  in  the  glass 
the  windows  of  the  enclosure  are  formed  by  entrance’  and  exit 
collimating  lenses. 

2 .  STUDY  or  Tin,  EMISSION  SPECTRUM  OF  THE  NIGHT  SKY 

The  intorpr  tat  ion  of  the  r.  cults  of  this  study  lias  al¬ 
ready  been  the  subject  of  two  publications  [61,62).  .Here  wc 
are  going  to  give  a  brief  description  of  the  conditions  under 
which  the  experiments  have  been  made. 
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Fig.  76.  Filter  Transmission  Curves;  the 
Ordinates  are  Transmission  Factors  in  Per¬ 
cent. 


The  setup  is  the  one  shown  in  Figs.  7 2b  and  74.  A  simple 
reflection  mirror  shown  in  Fig.  72b  permits  one  to  look  to  the 
northwest  at  20°  above  the  horizon.  Hence  there  is  no  optical 
intermediary  between  the  sky  and  the  entrance  lens  of  the 
interferometer.  The  diameter  of  the  hole  has  been  chosen  in 
such  a  way  that  the  apparatus  works  under  the  best  luminosity- 
resdution  conditions.  Under  these  conditions  the  angle  of  the 
maximum  opening  of  the  beam  when  we  wish  to  obtain  an  effective 
resolution  of  2,000  measures  3  x  10~*  radians  and  the  solid 
angle  of  the  beam  used  is  SI  ~  2.7  x  10-3  steradians.  The 
modulator  placed  before  the  entrance  hole  Bs  chops  the  flux 
at  a  frequency  of  125  cps.  The  electric  signal  coming  from 
the  lead  sulfide  cell  is  amplified  in  a  classical  synchronous 
detection  apparatus,  filtered  by  an  RC  filter  and  recording 
potentiometer. 

Wc  have  studied  the  region  near  1.6  p  with  a  lead  sulfide 
cell  in  the  region  near  1  p  with  a  Lallemand  photomultiplier. 

( 1 V-  Study  of  the  Region  in  the  Neighborhood  of  1.6  n  — 

A  filter,  placed  directly  before  the  cell  and  whose  trans¬ 
mission  curve  is  represented  by  Fig.  76a,  eliminated  wave 
lengths  below  lp.  As  reference  line  we  used  the  red  cadmium 
line  isolated  in  the  spectrum  of  an  Osram  lamp  by  colored 
filters.  The  infrared,  which  comes  from  the  Osram  lamp  and 
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which  can  reach  the  lead  sulfide  cell  was  cut  by  an  ensemble 
of  filters  whose  transmission  curve  is  given  in  Fig.  76b  and 
by  a  2.5%  copper  sulfate  solution  of  2  cm  thickness. 

Figure  77a  gives  an  example  of  an  interf erogram  obtained 
under  these  conditions.  Immediately  one  notices  an  anomaly: 
the  central  fringe  does  not  correspond  to  a  maximum.  But  it 
can  be  shown  that  in  an  interference  phenomenon  the  maximum 
intensity,  if  it  is  a  question  of  a  system  of  interferences 
with  a  brilliant  central  fringe,  is  always  obtained  in  the 
central  fringe.  ((4  5]  p.  659).  The  explanation  comes  from 
the  fact  that  the  recorded  interf erogram  is  taken  as  the  sum 
of  two  interf erograrns  with  opposite  phases:  one  which  comes 
from  the  emission  of  the  chopper  blades  in  a  black  box  whose 
temperature  is  about  15°C  while  the  cell  is  approximately  at 
-  80cC,  the  temperature  of  dry  ice.  Since  the  chopper  radi- 
ates  approximately  a;,  a  black  body,  the  radiation  spectrum 
which  it  emits  is  large.  The  corresponding  "thterf erogram  is 
reduced  to  some  fringes  in  the  neighborhood  of  the  zero  phase 
difference.  We  recorded  it  one  night  when  the  sky  was  covered 
and  when  the  emission  oi  the  night  sky  was  therefore  null  (Fig. 
78a).  The  corresponding  spectrum  (Fig.  78b)  is  cut  aft  the  long 
wave  side  by  glass  which  becomes  absorbent.  In  order  to  atten¬ 
uate  this  signal  strongly  in  phase  opposition,  it  is  sufficient 
to  place  before  the  coll  an  interference  filter  which  cuts  off 
beyond  2  n  and  lets  pass  the  radiations  near  1.6  p.  With  such 
a  filter,  the  interf erogram  is  returned  to  nearly  normal  with 
a  much  larger  central  fringe  (Fig.  79a). 

The  referencing  of  the  (  i o  point  can  be  done 'in  two  ways 
which  give  equivalent  precision:  cither  by  taking  the-  mean  ol- 
the  midpoints  of  the  segments  joining  th<  axes  of  the  two  sym¬ 
metric  fringes  with  respect  to  the  central  fringe  like  ab,  c«J, 
of,  as  indicatou  in  Fig.  77a,  or  else  by  registering  at  zero 
phastr  difference,  the  interf  erogram  of  a  white  source  with  the 
same  time  constant  (Fig.  80). 
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Fig.  78.  (a)  Interf erogram  of  the  Radiation 

-  Emitted  by  the  Modulator;  (b)  Spectrum  of  th« 
Radiation  Emitted  Sy  the  Modulator. 


We  give  below  three  spectra  obtained  at  resolutions  900 
(Fig.  81  and  /9b)  and  2,000  (Fig. 

Spectrum  of  Fig.  81.  Effective  resolution,  900.  Time 
for  recording  the  interf erogram,  2  hrs. 

Figure  81  gives  a  comparison  between  two  spectra:  one 
at  resolution  900  obtained  through  a  Fourier  transformation 
and  one  at  resolution  150  obtained  with  a  grating  spectro¬ 
meter.  This  last  is  in  fact  the  mean  of  eleven  spectrums, 
each  having  been  recorded  during  a  half  hour*16'5.  The 

<l6)Wallace  Jon-- s  has  recently  obtained  with  a  crating 
spectrometer  a  spectrum  at  resolution  400  by  using  a  new  germa¬ 
nium  detector.  In  the  two  experiments -compared  in  Fig.  81,  the 
detector  was  a  Kodak  cell. 
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Fig*  79.  (a)  Interferogram  of  the  Night  Sky  (February  28,  1959),  Lead  Sulfide  Cell 

With  Interference  Filter;  (b)  Spectrum  of  the  Night  Sky. 
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Fig.  80.  Interferogram  of 
the  Night  Sky  and  Whit* 
Source  Recorded  With  the 
Same  Time  Constant. 


Fig.  81.  Spectrum  of  the  Night  Sky: 

R  »  150,  — —  R  900.  (Bruit  "  noise] 


...  i 
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negative  maxima  that  one  sees  on  the  low  wave  number  side  of 
each  ray  which  has  a  strong  intensity  show  that  there  has 
been  an  error  in  determining  the  zero  point  of  the  interfero- 
gram. 

Spectrum  of  Fig.  79b.  Using  the  interferogram  of  Fig.  79a 
recorded  at  twilight  the  following  have  been  calculated: 

Maximum  phase  difference  obtained,  1.96mm; 
calculated  limit  of  resolution  taking  into  account 
the  width  of  the  beam,  6.25  cm'1; 

measured  limit  of  T'efsolution ,  6.3  cm-1; 
spectral  domain  studied,  1,250  cm'1  between  .5*550 
and  6,800  cm'1;  number  of  spectral  elements,  210; 
recording  time,  40  min. 

We  were  limited  in  che  choice  of  the  recording  speed  by 
the  response  time  of  the  recording  pen  which  took  0.8  sec  to 
go  from  one  extremity  of  the  scale  to  the  other.  The  limit 
speed  of  the  pen  hence  is  30  mm/sec.  With  the  recording 
frequencies  chosen,  i.e.,  0.5  cps  in  the  middle  of  the  spec¬ 
trum,  and  the  maximum  amplitude  of  the  fringes,  i.e.,  150mm 
for  zero  phase  difference,  the  maximum  speed  of  the  pen  was 
only  150  mm/se<j:.  It  is  very  much  lower  everywhere  else  in  > 
the  interferogram  and,  under  these  conditions,  the  recording  , 
function  can  be  considered  as  perfectly  linear. 

The  time  constant  used  was  r  =  0.125  sec. 

Figure  82  shows  the  relative 
position  of  the  spectrum  of  the 
Fourier  frequencies  and  of  the  noise 
spectrum. 

Conditions  under  which  the  F.T. 
was  taken:  numerical  F.T.  with  an 
interval  h  =  0.36  sec  (the  distance 

separating  two  points  corresponds 

! 

to  a  variation  of  a  phase  difference 
of  a  half  wave  length  of  the  red 


Fig.  82.  Relative 
Position  of  the  Noise 
Spectrum  and  of  the 
Spectrum  of  Fourier 
Frequencies. 
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cadmium  line).  Under  these  conditions  k/r  m  3  and,  according 
to  Fig. '57,  the  signal/noise  ratio  is  weakened  by  a  factor  0.7. 
The  total  number  of  points  taken:  6,120. 

Number  of  calculated  points:  840,  l.e...  4  per  spec¬ 
tral  element. 

Calculation  time  on  the  IBM  704:  1  hr,  20  min. 

Mean  square  error  of  the  fluctuations  of  the  inter- 
ferogram  a  «  7.6  u.  This  error  has  been  calculated  by  taking 
the  square  root  of  the  mean  of  the  squares  of  the  ordinates 
picked  off  the  interferogram  by  choosing  4  points  per  noise 
correlation  radius.  The  portion  of  the  interferogram  chosen 
corresponds  to.fi  ■  ^max'  ®ne  *s  certa*n  that  for  this  phase 
difference,  the  signal  in  the  interferogram  is  zero.  Actually 
the  diameter  of  the  hole  has  been  chosen  such  that  the  product 
&R  be  maximum  (III. 2).  Under  these  conditions  at  the  maximum 
phase  difference,  the  modulation  fell  to  aero. 

The  predicted  mean  square  error  of  the  fluctuations  in 
the  spectrum  (see  IV. 23b). . 

»x  »  ^0.9  X  4  x  0.406  X  X  -  211  u  sjL. 

Mean  square  error  measured  in  the  spectrum:  212  u  sec. 

The  calculation  was  made  in  the  neighborhood  of  3,500  cm’*1 
because  in  this  region  the  glass  is  absorbent  and  one  is  certain 
that  the  signal  there  is  zero. 

Such  good  agreement  must  be  considered  fortuitous. 

Spectrum  of  Fig.  77b.  Using  the  interferogram  of  Fig.  77a 
recorded  without  a  filter,  it  was  calculated: 

maximum  phase  difference  obtained,  L  ■  3.84  mm; 
theoretical  limit  of  resolution,  fia  •  3  cm”1; 
measured  limit  of  resolution,  3  cm"»; 
number  of  spectral  elements  studied,  400} 
recording  time,  2  hrs ,  30  min; 
time  constant,  r  »  0.25  sec; 
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k  “  0.72  sec,  whence  k/r  »  3; 

number  of  points  chosen,  11,928; 

number  of  output  points,  1,S00;  _ 

calculation  time  on  the  IBM  650,  150  hrs(1,) 

If  this  calculation  had. been  made  on  the  IBM  704,  one 
could  have  been  able  to  make  a  numerical -filtering  to  isolate 
the  interesting  region  of  width  »  1,100  cm"1  between  5,550 
and  6,650  cm"1.  In  this  case  the  calculation  of  the  convolu¬ 
tion  would  have  -been  made  using  11,928  input  points,  but  the 
interval  k'  chosen  in  order  to  make  the  F.T.  would  have  been 
only  of  the  order  of  1,000.  ” 

Mean  square  error  of  the -fluctuations  of  the  interfero- 
graa:  *  a  •  13.5  u. 

Mean  square  error  of  the  fluctuations  of  the  spectrum 


9  x  4  x  0.406  x 


9.000 

4 


\*/« 

180.8 j  «  1,000  u  see. 


Measured  mean  square  error:  1,012  u  sec. 


(2)  Study  of  the  Region  Around  1  a 

The  detector  is  a  Lallemand  photomultiplier  chilled  with 
dry  ice.  Hence  for  the  maximum  noise  reduction  (here  it  is 
photon  noise)  it  is  necessary  to  isolate  the  region  to  be 
studied  by  a  filter.  For  this  result  we  have  used  a  combina¬ 
tion  of  a  colored  filter  which  cuts  all  the  visible  up  to 
8,000  X  and  of  an  interference  filter  whose  transmission  curve 
is  given  by  Fig.  83.  On  the  long  wave  length  side,  the  spec¬ 
trum  is  limited  by  the  region  of  sensitivity  of  the  cell.  The 
interferograms  are  of  the  type  represented  by  Fig.  84a.  We 
calculated  two  spectra  using  the  same  interferogram,  one  using 


,(lT ’This  calculation  was  made  at  the  University  of 
Toronto's  IBM-650  calculator  which  has  magnetic  drum  and 
fast  ^access  memories.  ' 
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Fig.  83.  Transmission  Curve  of  the  Inter¬ 
ference  Filter  Centered  on  1  ji;  as  Ordinate, 
Transmission  Factor  in  Percent. 

*4 : 


an  almost  triple  length  of  the  interf urogram:  1.22  mm.  Total 
recording  time:  2  hrs. 

The  resolutions  obtained  are  respectively  of  the  order  of 
350  and  1,000  (Fig.  84b).  A  nearly  identical  spectrum  had 
been  obtained  with  a  grating  spectrograph  for  an  exposure  of 
48  hr. 

( 3 )  Resume  of  the  Results  Obtained  in  the  Stueiy  of  the 

Spectrum  of  the  Emission  of  the  Night  Sky 

The  conclusions  that  can  be  drawn  from  the  study  of  the 
spectra  have  been  discussed  elsewhere  [62];  here  we  shall  give 
only  a  resume: 

all  the  observed  lines  can  be  attributed  to  the 
radical  OH; 

the  rotational  structure  of  the  bands  (3.1),  (4.2), 
(5.3)  and  (6.4)  is  completely  determined; 

one  has  been  able  to  set  up  a  table  of  observed 
frequencies  for  the  vibration-rotation  lines  of  the  radical 
OH,  a  task  which  is  difficult  to  do  based  on  experiments  per¬ 
formed  in  the  laboratory  because  of  the  difficulty  of  construc¬ 
ting  suitable  sources; 

the  attempt  to  determine  a  rotation  temperature  using 
the  intensities  o l  the  lines  shows  that  agreement  cannot  be 
obtained  between  the  calculated  temperatures  using  branches 
P,  Q  and  R  simultaneously.  A  possible  explanation  is  that 
the  transmission  probabilities  calculated  by  Benedict,  Plyler 
6  Humphreys  [63]  and  used  in  the  determination  of  temperature 
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based  on  the  Intensities  of  the 
lines  are  false.  Besides,  the 
authors  had  predicted  this  pos¬ 
sibility  because  they  neglected 
the  vibration-rotation  inter¬ 
action. 

Vf  3.  STUDY  OF  THE  RE- EMISSION 
RADIATION  OF  GERMANIUM 

We  have  also  studied  by  this 
method  the  intrinsic  recombina¬ 
tion  radiation  in  germanium.  A,! 
complete  study  of  the  results  j 
obtained  will  be  the  subject  of  a  thesis  by  0.  Parodi.  Here  l! 
we  shall  confine  ourselves  to  presenting  the  first  spectra  j 
obtained. 

A  first  study  was  made  by  C.  Benoit  a  la  Guillaume  fc 
0.  Parodi  who  had  studied  with  a  grating  spectrometer  (Fig.  85) 
the  recombination  radiation  emitted  by  a  germanium  crystal  in 
which  an  excess  of  minority  carriers  had  been  created  by  opti¬ 
cal  injection.  At  the  temperature  of  liquid  helium*  they  dis¬ 
covered  three  lines.  Lines  A  and  C  are  due  to  the  annihilation 
of  an  indirect  exciton  with  simultaneous  emission  of  a  longi¬ 
tudinal  acoustic  photon  (A)  or  optical  tranverse  photon  (C). 

The  origin  of  the  line  B  which  appears  only  below  2*K  was  not 
exactly  determined.  We  have  been  able  to  prove  that  it  was 
a  line  due  to  recombination  at  centers  associated  with  dis- 
locations.  Its  intensity  varies  with  the  point  where  the  opti¬ 
cal  injection  was  made  5n  the  crystal;  it  disappeared  when  the 
sample  was  placed  in  a  magnetic  field. 

The  theoretical  limit  of  resolution  of  the  grating  spec¬ 
trometer  do.  was  4  cm"1.  The  width  of  line  A  measured  in  the  ' 

0 

spectrum  in  Fig.  85  equals  14  cm  . 

One  could  hope  for  several  advantages  of  the  F.T.  method: 


*VW 


Fig.  85.  Spectrum  of  Re¬ 
combination  Radiation 
corded  by  a  Gratir 
Spectrometer,  R  -  1,400, 
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a.  To  explore  a  more  extended  spectral  domain  between 

5,600  and  6,100  cm”1  without  increasing  the  measurement  time. 
In  fact  the  question  of  tfhis  “ime  is  of  first  importance  in 
this  experiment  because  lighting  the  crystal  heats  it  and  the., 
consumption  of  helium  necessary  to  keep  the  temperature  below 
2*K  increases  rapidly.  /'  .< 

b.  To  increase  the  signal/noise  ratio  which  is  necessary 
for  two  reasons: 

The  recordings  with  a  grating  spectrometer  did-  not  allow 
sufficient  precision  in  deducing  the  law  of  variation  of  the 
displacement  of  the  line  A  as  a  function  of  the  magnetic  field 
in  which  the  crystal  is  placed; 

2° .  One  could  hope  to  disclose  the  eventual  presence  of  lines 
around  5,950  cm”1  whose  existence  had  been  suspected  by  Haynes 
|59|; 

c.  To  increase  the  resolution  in  order  eventually  to 

determine  the  structure  of  the  line  A.  I 

Even  though  this  problem  may -not  be  a  problem  with  a 
great  number  of  spectral  elements,  like  that  of  the  night  sky^ 
it  appeared  desirable  to  approach  it  with  this  method. 

We  tried  several  ways  of  recording  the  interferogr&s. 

(1)  Recording  at  a  Constant _ Speed 

Since  the  spectrum  studied  is  much  simpler  than  that  of 
the  night  sky,  the  interf erogram  is  very  different  (Fig.  86a). 
It  has  the  appearance  of  a  sine  wave  whose  amplitude  is  slowly 
modulated: 

maximum  phase  difference,  L  *  2.1mm; 
theoretical  limit  of  resolution,  “5.7  cm’1; 

measurement  time,  40  min. 

In  this  case  again,  tile  recording  speed  is  limited  by  the 
maximum  speed  of  the  recording  pen 
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Fig.  86.  (a)  Interf erogram  Recorded  at  Con¬ 

stant  Speed,  Rq  ■  1,000;  (b)  Spectrum  Obtained 

From  This  Interf erogram. 

The  spectrum  calculated  numerically  from  this  interfero- 
gram  is  shown  in  Fig.  86b: 

width  at  mid-height  of  the  line,  14  cm’1; 
signal/noise  ratio  for  line  A,  of  the  order  of  100. 

Again  one  finds  the  lines  A,  B,  C. 

Furthermore  the  presence  of  a  line  D  is  disclosed  around 
5,655  cm  1  which  can  also  be  interpreted  as  a  dislocation  line 
(it  appears  only  below  2°K  and  disappears  in  the  field)  and  of 
a  line  E  around  5,960  cm’1  whose  origin  is  not  explained  and 
which  apparently  is  double.  In  order  to  make  sure, it  is  nec¬ 
essary  to  increase  the  resolution,  hence  the  maximum  phase 
difference  obtained.  It  cannot  he  a  question  of  making  a  re¬ 
cording  at  a  constant  speed  because: 

1* .  The  measurement  time  would  be  much  too  larga  with  the 
recorder  we  used;,  — — •  —  T  — 


2* .  The  number  of  points  to  pick  from  the  interf erogrrem 
would  be  too  large. 

Two  solutions  can  be  considered:  recording  at  a  variable 
speed,  the  method  of  changing  the  frequency. 

(2)  Recording  at  Variable  Speed-—  _ 1__ 

We  have  already  explained  in  the  preceding  chapter  that 
this  system  had  been  tried,  had  correctly  functioned  (Fig.  87) 
but  that,  being  given  the  characteristics  of  our  mechanical 
setup,  the  useable  time  was  not  a  sufficient  fraction  of  the 
total  time  and  the  s/n  ratio  was  insufficient. 


(3)  Method  of  Changing  the  Frequency 

This  procedure  allows  the  recording  of  lower  electrical 
frequencies  (hence  an  exploration  speed  which  can  be  con¬ 
siderably  increased  and  a  measurement  time  very  diminished) 
and  at  the  same  time  the  bringing  into  better  coincidence 
the  noise  spectrum  with  the  spectrum  of  the  Fourier  frequenelea 
(hence  a  smaller  number  of  points  to  choose) .  — 

The  spectral  domain  occupied  is  Ac  ■  S00  cm**1  between 
5, €00  and  JB^IQO  cm’ 1 .  If  one  chooses  the  heterodyne  wave 
number  <r(  67160  cm" 1 ,  the  spectrum  to  be  studied  after  chang¬ 

ing  frequency  is  found  between  0  and  500  cm*1  and  can  be  ob¬ 
tained  by  a  simple  cosine  transform  of  the  new  interferogram 
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Fig.  68.  Spectral  Domain  Occupied  Before  and 
After  Changing  the  Frequency. 

(Fig.  88).  For  the  sarpe  speed  of  the  moving  carriage  the 
maximum  frequency  of  the  signal  to  be  studied  will  then  be 
divided  by  12.  i 

Recall  that,  being  given  the  irregularities  in  the  car¬ 
riage  drive  and  frequency  instabilities  which  _re suit  from  it, 
there  is  no  question  of  using  an  electrical  signal  furnished 
by  a  low- frequency  oscillator  in  order  to  make  a  change  of 
frequency.  It  must  be  furnished  by  the  interferometer  itself 
under  the  same  conditions  as  the  signal  to  be  studied.  We 
have  not  used  the  signal  coming  out  of  an  optical  line  of 
wave  number  6,100  cm”1  but,  noticing  that  the  wave  number  of 
the  green  mercury  line  is  18,312  cm”1  and  that  18,312  cm”V3  * 
6,104  cm”1,  one  finds  that  it  is  sufficient  to  use  tie  third 
harmonic  frequency  corresponding  to  the  green  mercury  line. 

SOME  DETAILS  OF  THE  EXPERIMENT.  The  diagram  of  the  set¬ 
up  is  shown  in  Fig.  89,  the  apparatus  itself  in  Fig.  90. 

Reference  Line.  As  reference  line  we  have  us^d  the 
green  mercury  line.  The  signal  furnished  by  the  photomulti¬ 
plier  is  sent  on  the  one  hand  to  a  s.  ■  of  multivibrators 
which  divide  its  frequency  by  3,  on  the  other  hand  to  the 
writing  relay  of  the  recorder.  Moreover  this  relay  can  write 
cither  a  square  pulse  for  each  reference  fringe  or  a  signal 
having  a  frequency  three  times  smaller. 
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pressure  mercury  arc  (lamp),  chopper,  photomulti¬ 
plier,  frequency  divider,  fringe  counter,  writing 
relay;  lead  sulphide  cell,  automatic  gain  control, 
amplifier,  Inverting  relay,  low  pass  filter,  re- 
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j~  Fig.  91.  Sauare^Pulse  Function. F(t) 

of  Period  l/»#. 


Fig.  92.  Spectrum  Obtained  by  the  Coslam 
Transform  of  the  Product ~X(t)  P(t). 

Change  of  Frequency.  We  have  not  multiplied  the  initial 
. intorf erogram  by  a  sinusoidal  signal  cos  Srvt,  but  by  a  square 
pulse  function  F(t)  (Fig.  93)  which  one  <ja n  put  in  the  form 

F(t)  *  cos  2»v#t  -1  cos  2jt  3i»#t  +  ^  cos  2»  Sk#t  +  ...  , 

which  is  much  easier-  to  generate  because  all  one  has  to  do  it 
invert  the  signal  at  equal  time  intervals  1 19  * . 

Hence  after  multiplication,  the  frequencies  to  be  analyzed 
occupy  the  intervals  0,  being  the  maximum  frequency  cor¬ 
responding  to  <tk  *  500);  3  3  k#;  5  ^  5  (Tig. 

' 1C 'The "inversion  of  the  signal  would  take  place  at  equal 
intervals  of  time  if  the  advance  of  the  carriage  were  perfect; 
practically  it  is  done  at  time  intervals  corresponding  to  the 
passing  of  1.9  fringes  oi  the  mercury  11ns. 
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Fig.  93.  (a)  Interf erogram  Obtained  by  the 

Method  of  Changing  Frequency;  (b)  Spectrum  De- 
duced~Trom  This  Interf erogram,  R#  «  6,000. 

[vitesse  lente  (rapide)  *  slow  (fast)  speed) 

It  is  necessary  to  take  into  account  the  parasite  spectra  in 
the  choice  of  the  interval  h.  Their  frequency  is  higher  than 
that  of  the  terms  in  vQ  +  which,  as  we  explained  in  the  pre¬ 
ceding  chapter,  are  filtered  by  a  low-pass  filter.  Moreover  we 
shall  make  a  numerical  filtering  which  completely  annuls  their 
effect. 

•i  '  .  *  * 

Suppression  of  the  Mean  Level.  The  signal  furnished  by 
the  lead  sulfide  cell  includes  the  modulated  part,  which  we 
call  the  interf erogram  and  with  which  we  have  been  concerned 
up  until  now,  and  a  constant  component  which  after  changing  the 
frequency,  gives  an  important  signal  at  the  frequency  v  . 

Hence  it  is  necessary  to  eliminate  this  constant  signal  by  add* 
ing  to  it  an  equal  signal  in  the  opposite  sense.  A  residue  of 
the  mean  signal  on  the  interf erogram  gives  rapid  oscillations 
which  are  superposed  on  the  slow  fringes  of  the  new  interf ero¬ 
gram  (Fig.  93a).  This  signal  at  the  frequency  is  moreover 
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Fig.  94.  Relative  Position  of  the  Noise  Spectrum 
and  of  the  Fourier  Frequency  Spectrum,  r  ■  1  sec. 

strongly  attenuated  by  the  low- pass  filter  which  is  supposed 
to  eliminate  completely  the  terms  of  frequency  v  + 

Locating  the  Zero.  In  order  that  the  marking  be  done 
with  sufficient  precision  in  the  neighborhood  of  zero,  the 
speed  is  16  times  smaller  than  the  normal  recording  speed;  -i 
at  the  passage  of  the  white  fringe  furnished  by  a  continuous  ' 
source,  one  releases  a  counter  which  counts  the  fringes  of 
the  mercury  line.  Its  indications  set  off  the  frequency  multi 
plication;  the  same  movement  divides  by  3  the  frequency  of  the 
writing  relay  of  the  reference  fringe.  A  last  operation  pro¬ 
vides  that  at  any  instant  what eve?  one  may  pass  to  a  speed  16 
times  larger. 

R*  . orded  Interf erogram  (Fig.  93a).  It  has  not  been  re¬ 
corded  under  the  best  conditions: 

maximum  phase  difference  attained,  L  ■  11.7mm; 
theoretical  limit  of  resolution,  6e#  •  1  cm*1; 
recording  time,  1  hr,  30  min;  :  — 

electrical  frequencies  of  the  signal  to  |>e  studied 
taken  between  0  and  0.1  cps  (Fig.  94); 
time  constant  r  *  1  sec.’ 

Spectrum  of  Fig.  93b.  Figure  93b  shows  the  spectrum  ob¬ 
tained  from  the  interf erogram  of  ^ig.  93a: 

value  of  the  interval  fc,  0.5  sec,  where  h/r  m  0.5 
and  the  s/n  ratio  in  the  spectrum  is  practically  the  same  as 
if  one  had  used  all  the  points  of  the  Interferogram; 
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Fig.  95.  Intrinsic  Recombination  Radiation  in 
Germanium,  Rfl  »  6,000. 


numbeiLJ^jixaiJxt-S  chosen,  12,672; 
number  of  calculated  points,  401.  " 

A  numerical  filtering  (width  of  the  filter  was  chosen 
equal  to  500  cm-1)  resulted  in  the  reduction  of  the  number  of 
points  used  for  making  the  F.T. 

The  calculation  was  made  only  in  the  regions  of  the  lines 
A,  C  and  E.  Calculation  time:  22  min. 

An  examination  of  the  spectrum  shows  two  interesting 
results:  ■_  5 

the  line  A  is  distinctly  asymmetric; 
the  line  E  is  divided  into  two  lines  at  5,952  and 
5,963  cm*1; 

detailed  examination  of  the  line  (Fig.  95)  shows 
that  its  width  is  of  the  order  of  2  cm”1;  thus  it  has  been 
studied  with  an  excessive  resolution.  Also  in  later  experi¬ 
ments  whose  results  are  not  given  here,  we  contented  ourselves 
with  a  resolution  of  the  order  of  4,000.  Since  the  noise  in 
the  spectrum  of  fig.  93b  appeared  to  us  to  be  due  in  a  great 
part  to  the  fluctuations  of  the  mean  level  of  the  interf erogram, 
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we  stabilized  the  mercury  arc  and  reduced  the  potential  of  the 
lead  sulfide  ce'l  as  is  indicated  in  Fig.  89. 

4 .  CONCLUSION 

The  use  of  the  Fourier  transformation  method  allowed  us 
to  go  from  a  resolution  of  150  to  a  resolution  of  2,000  in  the 
study  of  the  radiation  emitted  by  the  night  shy.  In  the  case 
of  the  recombination  radiation  in  germanium,  a  problem  with  a 
small  number  of  spectral  elements,  the  results  are  less  spec¬ 
tacular.  Hovi£M&¥-, r— we  have  been  able  to  demonstrate  the  asymmetry 
of  the  line  due  to  the  annhiliation  of  an  cxciton  with  emission 
of  «  longitudinal  acoustic  photon  and  the  presence  of  two  lines 
around  5,960  cm"1  whose  origin  is  not  accounted  for  yet. 
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VII.  COMPARISON  OF  DIFFERENT 
SPECTROSCOPIC  METHODS 

The  Fourier  transformation  cpeqtroscopic  method  that  we 
hav<  just  treated  requires  very  different  techniques  from 
those-  of  classical  methods.  However  there  are  no  fundamental 
dif l'c  r-r-nces  among  all  the  methods.  Whatever  method  is  used, 
the  reproduced  spectrum  is  always  the  cosine  transform  of  the 
even  part  of  the  autocorrelation  function  C (6)  of  the 

incident  luminous  vibration.  It  is  this  function  v(S)  that' 
up  till  now  we  have  called  "interf erogram" . 

Wo  arc  going  to  show  how,  in  classical  methods  (prism, 
grating,  or  Fabry-Pri'Ot  spectrometer)  the-  flux  that  one  mea¬ 
sures  at  the  exit  of  the  instrument  is  the  F.T.  of  ?(*),  Ob¬ 
tained,  for  each  frequency,  as  the  sum  of  a  particular  series 
of  n  equidistant  discrete  values  of  £(6),  n  being  the  number 
of  interfering' beams .  In  the  Fourier  transformation  method, 
n  -  2;  the  flux  that  one  measures  for  the  value  h  of  the  phase 
difference  is  the  sum  of  only  two  particular  values  of  T’(6): 

£*(0)  and  and  it  is  necessary  to  make  the  Fourier  trans¬ 

formation  a  posteriori. 

One  particularly  rapid  way  of  making  a  synthesis  of  dif¬ 
ferent  spectroscopic  methods  consists  in  applying  to  the  treat¬ 
ment  of  the  luminous  vibration  by  a  spectrometer,  some  of  the 
more  elementary  results  relative  to  the  determination  of  the 
spectrum  of  fluctuating  magnitudes.  Indeed  the  electromagnetic 
field  at  a  point  in  space  is  a  random  function  of  time. 

Follgett  (14,1:;],  Gobble  ( 2u  1  and  Hertz  ( t  4  ]  have  already 
brieilv  indicated  this  ..av  of  treating  sp  ctroseopie  instru¬ 
ments.  The.  the  ore  ti  'al  problems  relative  to  the  random  magni¬ 
tudes  And  to  their  spectra  have  been  treated  in  a  satisfactory 
way  for  the  first  time  by  Wiener  and  Levy  in  a  study  of  Brownian 
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movement.  Since  then,  these  studies  have  been  pushed  very 
far  from  a  purely  mathematical  point  of  view  and  the  results 
have  been  applied  to  the  theory  of  turbulence  and  to  the  treat- 

•  aent  of  background  noise  in  the  theory  of  communications. 

~  Before  developing  the  particular  case  of  optical  spectro¬ 

scopy  we  are  going  .to  recall  some  definitions  and  theorems 
concerning  the  spectral  analysis  of  a  random  variable. 

1.  ROLC  OF  ft  LINEAR.  HONOGElItCUS  FILTER  IN  SPECTRAL 
ANALYSIS .  • 

One  can  associate  with  a  variable  \(t)  ( deterministic  or 
random  f unction  of  time]  a  power  spectrum  only  if  the  mean 
power  E{|X(t)|*},  ilcfincd  as  the  arithmetic  mean  of  the  quan¬ 
tity  |x(t)  |*  measured  at  tl  >  ■same  instant  t  for  an  infinity  of 
trials  macroscopically  identical,  is  independent  of  time.  Then 

<  for  each  frequency  component  .!'  obtained  using  X(t)  lit  e  Fourlmir 

i.  ...  .. 


b(r)  e2****  dr 


'*  •  j  * 


it  will  be  possible  to  assign  a  definite  B(r).  >r.s 

The  mean  power  is  tho  sum  —  .  ? 

jCVII.l)  E  {|X(t)|*}  -  B(i»)  dr  . 

Wc  are  interested  only  in  random  functions  fulfilling  these 
conditions;  there  are  tho  functions  one  meets  in  optics.  They 
are  functions  with  stationary  covariance  and  their. study  rests 
for  a  larcc  part  on  the  theorem  of  Khin<;hin  according  to  which 

the  F»-T.  of.  th*.  -'crrolation  function  C  tr)  *  E{x(t)  **(t-r)}  Is 

*  .  \ 

prooiicjy  B(r): 


Vs  '■  •  Jr , 
*v  .  v~ 


'  (VII. 2) 


B(r)  *  j  ^  Cx(r)  *'2*,,’T  dr 


Vs>  . 


1 


HAVWKPS  REPORT  8099 


A  spectral  analysis  is  made  with  the  aid  of  a  filter.  We 

shall  call  the  linear  filter  «T  ((4a]  p.  347)  a  linear,  homo* 

*geneous  transformation  with  respect  to  time  which  makes  Y(t) 

>  •  .  ' 
correspond  to  the  variable  X(t)  in  such  a  way  that 

1  449 

(VII.3')  Y(t)  -  X(t)  *  P(t)  -  \  X(i)  P(t-S)  dd  , 

P(t)l  being  the  percussional  response  of  the  filter,  i its 
response  to  the  Dirac  impulse 

The  complex  quantity  G(r)  with  modulus ,g(p).  which  is 
called  filter  gain  is  the  F.T.  of  P(t)  .  • 

(VII. 4)  C(k)  -  J  P(t)  dt  . 

Two  classical  calculations  lead  to  putting  the  mean  power, 
after  filtering,  in  one  or  the  other  of  two  fundamental  forms* 

v  ■*'  y~ 

(VII. 5)  E  ||Y(t)|#}  -  ^  »<r)  |C(r)|*  dr  ,  *  7  ^ 


<l#,In  the  electromagnetic  theory  of  light  (65]  a  ligh  • 
pulse  ii.  a  train  of  very  short  waves,  which  getsslu1.  *er  when  xts 
distribution  into  given  frequencies  given  by  a  Fourier  expansion 
gets  larger  (Pig.  96), 


Fig.  96.  Light  Impulse  and  its  Fourier 
Decomposition. 


I 
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(VII. 6)  E  Cx(r)  Cp(r)  dr  , 

Cp(r)  being  the  autocorrelation  function  of  P(t). 

These  two  equations  are  equivalent;  moreover  one  can  pass 
from  one  to  the  other  by  Parsoval’s  relation,  after  noticing 
that  theF.T.  of  Cx(fT  is  the  spfectral  decomposition  B(m)  looked 
for  and  that  the  F.T.  of  Cp(r)  is  thr  square  of  the  modulus  of 

the  gain  of  the  filter  jf.  . 

Moreover  we  can  notice  that  when  the  function  P(t)  is 
real,  it  is  simply  the  real  part  RcCC^r)}  «^(r)  which' enters  - 
intc  Eq.  (VII. 6);  it  then,  becomes  4 ) 

(VII. 7)  E  ||Y(t)f*]  »  Re[cx(r)]  Cp(r)  dr  .  . 

The  real  part  of  a  correlation  function  being  even,  its  F.T. 

Is  also  an  even  f miction  and  the  relation  (VII. 51  then  takes 


■*  -  . 


'In  fact 


.I",  EjlYa)!*}  -  E|Y(t)Y*(t)}  -  e||  lY(t)Y*(t)  +  Y*(t)Y(t)lj 

..  *  E{  2  ii  X(«)P(t-e)dd>.*(«')P*(t-v‘)dO’  + 

'  +2fl  X*(O)P*(t-0)dOX(0’)P(t-«,)de,j. 

In  the  case  where  P(t)  is  a  real  function,  P(t-#)  ~  P*(t»P) 
and  one  can  group  terms:  ■  *> 

.  EjlY(t)i*V  - 

pVt-oyPU-fl’Jdpdo’  -  \  \  [cx(r)+C5^-T>]cp<T^dT  ■ 

‘  i.r  ReK<r)]  CPlT>  dT  * 


the  form_ _ .. 

(VIZ. 8)  B  *f|v(t)  |*^  -  \  *"  B(r>  |C(«)|#  dr  , 


Bp(p)  being  the  even  part  of  the  spectrum  B(a)  as  was  defined 
In  Chapter  II. 

Formula  (VII. 8)  shows  that  the  mean  power  is  the  chopped 
energy  in  the  spectrum  times  the  energy  transmission  curve  of 
the  filter  ICf  r) 1 * .  ^^IfJHP  is  a  rectangular  band-pass  filter, 
one  has  10(11)1*“  1  where  v  is  taken  between  *•  ^  and  v*  and 
G(r)  “  0  elsewhere.  The  ideal  filter  for  making  the  spectral 
analysis  will  have  zero  width  but.,  if  i»  and  v  *  are  vjiy  close 
together one  can  say  that  “the  measured  power  is  proportional  • 
to  the  spectral  density  corresponding  to  a  frequency  taken 

between  ¥.  and  .  Then, 

V-  » 

'  * '&1L  4  y  ■ 

(VII. 9).  E  J|V(t)|a|  -  B(r#)  t*  . 

One  explores  the  spectrum  by  varying  the  band-pass  of  the  filter 


2.  GENERALITIES  ABOUT  THE  SPECTROMETERS  CONSIDERED 

XlTfilters 

A  cla*  >ical  spectrometer  (prism,  grating,  Fabry-Ferot)  or 


^  two-boar  interferometer  used  as  a  spectrometer  is  composed 
of  an  optical. part  which  at  ev**ry  instant  receives  luminous 
radiations  and  transmits  a  certain  proportion  of  them  to  a 


detcctor^amplif icr-measuring  setup* which  is  provided  with  a 
time  constant  and  measures  the  moan  power  carried  by  the  trans¬ 
mitted  radiations.  Let  us  first  consider  the  case  of  the 
spooXtonieter  .used  with  infinitely  thin  .  ntranm.  and  .exit  dia¬ 
phragms,  the  useful  “practical  VaSe"  hcTn?/  •  V.sily  deduced  from 
this  (Fig.  97a).  Let  X(t)  be  the  value  o£  the  electric  field 
in  the  interior  of  the  entrance  diaphragm  (point  source)  and 


■  «/•_..  4...  vx'ji  r  ■*.*  •  "VV  .. 

'V  *;**..*■  t  4*  '*  *+•  ■*  •>  j"V •  »v> 


-kc— g>i 

tMn'traf 

,w  x«**  -»  ••  »■«*»-  ,  <“•**“*"«  1  -.««.»••» 

.  "  ^  1  -  u|4Ktrril  dr  rn^tirr  | 

F13.  97.  (a)  Schematic  Diagram  of  the  Op¬ 

tical  Part  of  a  Spectrometer;  (b)  Schematic 
Diagram  of  the  Operation  of  a  Spectrometer. 

Y(t)  the  value  of  the  olectrr.c  field  in  ti.e  interior  «f  the 
exit  diaphragm  where  the  detector  is  placed  (Tig.  J/b) .  The 
signal  furnished  by  tnc  latter  is  amplified  then  transmitted 
to  a  recorder,  the  ensemble  having  a  time  constant  so  that 
one  measures  approximately  the  quantity 

t+'c  V- 

\  |V(t)|*  dt ... 

4t  . 


We  shall  suppose  that  in  what  follows  that  the  time  X  is  Pig^: 
enough  so  that  we  need  not  distinguish  between  *  -  • 

r t+  *  i  — - - r  J 

j  |Y(t)r  dt  and  the  mean  power  |Y(t)|a  flwhich  will  be 

rigorous  only  if  *  -*«•)  r~  This  approximation completely  neglects 
noise,  which  we  do  not  consider  in  this  chapter.  The  Laplacian 
stationary  character  of  X(t)  being  given,  | Y<tr > | *  has  the  same  „ 
value  as  c{ I Y(t) I*} . 

The  problem  is  to  determine  with  wnat  approximation 
E{  I  Y<t  )-!*},  which  we  shall  henceforth  call  *,  the  flux  measured 


by  the  detector1 


Is  the  spectrum  of  X(t)  such  as  wo  have, 


4,1,In  the  theoretical  case  considered  here,  the  diaphragms 
arc  infinitely  th in  and  ♦  is  a  line,  tn  all  practical  casas, 
the  beam  will  have  width  U  and  one  will  measure  a  flux* 


defined  it  in  paragraph  VII. 1,  i.e,,  to  what  approximation 
the  instrument  takes  the  F.T.  of  the  autocorrelation  function 
of  X(t).  . 

(1)  Expression  for  the  Flux  ♦ 

The  spectrometer- breaks  up  the  incident  radiation  into 
different  components  (spatially  in  the  case  of  grating  or  prism 
spectrometer,  amplitudinally  in_  the  case  of  the  Fabry-Perot 
spectrometer)  which  it  transmits  with  optical  delays  in  arith. 
me tic  progression  and  different  attenuations  before  combining 
them  in  such  a  way  that  Y(t)_ean  be  put_in  the  form 

(VSX.10)  Y(t)  -  X(t)  ±  X(t-rt)  +  ;a  X(t-2T|)  +  „.  ♦ 

'  +  x[t  -  (n-l)  Tj)  , 

where,  kr^  is  the  delay  of  the  (k+1)  beam  and  f> ^  is  its  trans¬ 
mission  factor.  This  sum  is  equivalent  to  the  Integral*1*  * 


(VII 


/■»  /»  +• 

1. 11)  Y(t)  -  j  X(t-Th.(r)dr  «  J  X(t-r)P(r)R(r)dr  » 


R(r)  being  a  Dir^x.  distribution  of  periodic  support’ of  interval  ~ 
Tj  and  P(r)  a  distribution  function  which  takes  on  the  particu¬ 
lar  values  j>],  fn  for  t  *  0,  r^,  ....  (n-l)  r%  and  can 

eventually  become  zero  for  t  >  n  rt  if  the  number  of  transmitted 


beams  before  the  recombination  is  finite, 
is  commutative ,  (VII. 11)  is  equivalent  to 


>ince  a  convolution 


(VII. 12) 


Y(t) 


X(f)  j>(t-i>  df 


In  this  form,  the  analogy  with  (VII. 5)  is  complete; 
is  the  percuss ional  rOspons-i  of  the  spectrometer  considered 


In  the  case  of  the  prism,  the  number  of  interfering 
beams  is  infinite  and  the  interval  of  tha  function  R(r)  is  sero. 
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as  a  linear  filter.  It  is  £  real-valued  function.'  Therefore 
one  can  write  immediately,  by  comparison  with  (VII. 8)  and 
(VII. 7),  ■  1  . 

r 

♦  *  Bp(*0  lGCa») I*  d*  and  ♦  »  j  ^(r)  CpCr)  dr  . 

Noticing  that  v,  t,  r  and  r  are  related  to  the  wave 
number  a,  to  the  length  l  measured  in  the  wave  train  and  to 
the  geometric  delays  6  and  ft  by  the  relation  <r  =  » >/ 1 ,  l  *  tc, 
fi  *  rc  and  *  »  r^c,  ]the  two  preceding  equations  take  the  form 


(VII. 13) 
and 

* 


*-r 


B_(a)  |G(m)|*  dm  , 

W  »  *  *  *' 


*  - 1,1 

A..-  . 


c 


£*(«>  Cp(6)  dA  . 


"  "  These  two  equations  have  been  established  tor  a  particular 
value  of  delay  ft  between  two  consecutive  beams;  the  functions 
P-{/)  and  p(Z)  and  the  flux  are  functions  of  ft.  In  the_  follow¬ 
ing  wc  shall  use  the  two  relations  (VII.  13)  aiti  (VU.34)  (which 
arc  in  fact  equivalent].  It  is  in  considering  th»  a  i  ond  rela¬ 
tion  that  we  shall  see  how  the  instrument  makes  the  T.T.  of 

an  operation  whose  result'' is  interpreted  by  the  first . 

relation!. 

The  function  Cp(d)  is>  the  autocorrelation  l unction  of 
the  product  P(Z)  Rft(J).  Hence  it  Is  still  ,thc  product .of 
Dirac  distribution  of  periodic  support,  of  the  samT  interval 

ft ,  with  an  even  l  unction  A(4)~whic)t  i3  the  autocorrelation 
function  of  P(l)i 

(VII. IS) 


Cp(A)  -  *(#)  R^A) 
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W e  shall  see  in  whait  follows  for  each  type  of  spectronieter , 
the  particular  form  of  A(d)  that  wc  have  been  talking  about 
up  till  now  will  play  the  same  role  as  the  weighting  function 
introduced  in  the  Fourier  transformation  method. 

Replacing  Cn(d)  by  its  value  in  (VII. 14)  we  gat 
P 


(VII. 16)  ♦(*)  «  j  C^(«)  . 

which  is  equivalent  to 

(VII. 17)  ♦<»)  -  A,  +  2  Aa  *  ...  +  2  fnml  An  - 


k=n-  j 

s 


kc-(n-i) 


v’k  \+J  * 


It  is  a  finite  sum  of  discrete  values  of  the  autocorrelation 
function  of  the  incident  radiation,  multiplied  by  the  weight¬ 
ing  function  A(d)  appropriate  for  the  instrument.  It  is  easy 
to  show  that  this  is  a  Fourier  transl ormution.  In  fact,  noth 
ing  changes  in  the  value  of  the  sum  (VII. 17)  if  one  multi¬ 
plies  each  term  by  the  factor  cos  2rafck*  in  which  has  been 
chosen,,  equal  to  m  being  an  integer,  in  such  a  way  that 

the  coaims  ate  always  equal  to  1. ''The  new  sum 

k»n+» 

(VII.18)  ♦(»)  -  2  cos  JfffjU 

k— (n-») 

f  • 

is  equivalent  to 

(VII. 19)  ♦<*)  -  \  *%*{*>  *<4>  R,,(a>  cos  2 toh6 


This  is  exactly  the  expression  (11.12)  giving  thi  value-  uf  a 
spectral  density  calculated  by  taking  the  numerical  F.T.  ol 
the  inter! erogram.  By  putting  F(e)  ■  T pOS  |A(6)  **(&)!»  one 
can  write  (VII. 19)  in  the  fora 


£/'>*►**  -  ^ *•■'.■»*:  *’•  *’*  1 


<>A\  ■ 


,  *  .  L  »  . 

«  i  i 
JL:  JL<  JL 

o 

Fig.  98,  Apparatus-Function. 
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tr 

T 


.U.L 


(VII. 20) 


.  r  +* 

*“  2  J  (F(*fc-»)  +  FOr^-fa)]  dir  • 


./  A* 
.!•  /  *  ^ 


\t  '.  ;* 


:M*?  • ; 


-In  comparing  (VII. 13)  and  (VII. 20)  one  sees  that 

♦  i  , 

(VII.21)  |G(c)  j*  «  |  +  FC^-to))  .  'I 

The  spectral  density  measured  lor  the  wave  number  hence 
is  the  chopped  energy  in  the  spectrum  multiplied  by  tluv  function 
|c(«r)  |a  which  spcctroccopists  call  apparatus-functia*. 

*  .  .  •  *.  w 

'ji 

(2)  form  of  tlie  Apparatus  "Function  _ 

Vkc  have  emphasised  the  analogy  between  the  expressions 
(VII.  3  9^)  and  (II. li).  1’i.uilarly  the  equations  (VII.21)  and  (11.14) 
show  in  the  classical  speelr^mcteTs  as  in  the  ease  of  the 

numerical  F.T.  of  an  inter U r.ujram,  the  apparatus-function  is 
composed  of  two  scrips  of  peaks  separated  by  a  distance  l/* 
land' centered  about  ^  and  -  9^  ,  each  one  of  the  peaks  being 
tho'r.T.  of  the  weighting  function  A(4)  o t  the  interferogram. 
'However  th:r.  arc  tone  little  .ill Urcr.ce:.  in  the  easr  where  the 
nujaoiiiw.1  r.T.  i*  taken  as  was  indicated  in  Chapter  II.  In  the 
case  pi  a  classical  i pi vtrometc r  since-  efc,  the  wave  number  for 
which  the  F.T.  is  made,  is  an  integer  multiple  of  l/*,  there 
always  is  a  peak  centered  about  the  origin  (Fig.  Vtt)  and  the 


y  ■ 
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Orating  Functions  of  Different  Intervals. 

two  series  of  peaks  of  the  apparatus- function  are  nixed  to¬ 
gether  in  such  a  way  that  one  can  write 

< VII. 22)  B  <<rfc)  »  j  •  Bp(e)  FO^-e)  d»  . 

When  one  explores  the  spectrum,  that  is  when  one  looks 
lor  the  spectral  density  corresponding  to  a  number  jh,  in  the 
neighborhood  of  ofc ,  one  changes  the  delay  between  two  consecu¬ 
tive  brain-',  which  becomes  k  • .  Hence  one  changes  the/  interval 
of  the  grating  I  unction  U^(d)  which  serves  to  pick  off  points 
in  the  lunction  ^(i)  1 1’ ig .  99).  The  new  value  B(ofcI)  is 
tluivlorc  obtained  from  a  new  scries  of  values  of  {^(4).  The 
pew  apparatus-function  is  not  deduced  from  the  first  by  two 
translations  as  is  the  case  when  one  makes  the  numerical  F.T. 
of  an  interfejogram  using  a  single  series  of  equidistant  values 
picked  off  on  this  interferogram  as  explained  in  Chapter  II. 

The  distance  ¥  between  two  consecutive  peaks  is  then  equal  to 
1  * ,  L-ut  one  of  the  peaks  remains  centered  about  the  origin 

(Fig.  100).  Moreover,  in. the  strictest  sense,  A(6)  is.  . unc¬ 
tion  of  h  and  the  form  of  aach  maximum  /aries  with  a;  of  course 
the  variations  will  be  nagllgible  in  the  spectral  domain  ex¬ 
plored  . 


I 

I 
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Fig.  iOO.  Deformation  j}f  the  Appot  at  us- Function 
With  the  Value  fc  of  the  Interval;  One  of  the  max- 
ilia  Remains  Centered  on  the  £ero  Frequency. 


We  have  <  n  that  wh<  ivr:,c  t  aK  s  tv  i<ui:u  rieal  r.T.-  oi 
;  an  interfnrogram,  one  ciioo-c. precisely  the  interval  /•  in  such 
a  way  that  a  "ingle  peak  explores  the  sp<:  -  v rum.  Tn  the  case 


V# 


'  ■*’ 
"  #■ 


■*  a-':: 


,  it?: A 1* 


of  classical  spectrometers,  it  is  m  to  take  the  pro- 

"«*.i,  ...  ■  ,«* 

‘'caution  eventually  to  limit  by  an  optical  filter  it hc  spectrum..’ 

<3  .  •  •’  *'■  :’<V 

**tp  be  studied  in  such  a  way  that  th..  spectral  interval  occupied. 


,'C% 
•V 


I'.Vlff  l>e  at  most  equal  to  1  /h  (Tig.  18). 

In  general  one  calls  the  limit  of  i 


...  y . —  — ...  — ....  ....  - -  ...  ■■  elution  thf  width  «lt 

-height  of  one  ol  the  peaks  of  the  appjratu.' -function.  *  This 
•"y**!  ’  '•  width  will  be  the  narrower  the  wider  its  F.T.  j>.‘‘ A(» ) ,  is,  Iw  nee 

j-”  •  .  ••  •*  -j£ v*.  i  i  'v,’T*v  v 

••-€•5  one  will  Krow  a  greater  portion  of  its  autocorrelation  function* 

Xt*--  ,  “  ■  <'  '  '  1*7. 

3,  FORM:.  OF  APPAKATUC-T UNCTIONS  OF  CLALPU  AL  riMXTROHFTkR:! 

i  It  suffices’to  determine  p(  O  in  each  typy  «q  opt  etro-  . 

...•meter  in  order  to  discover  immediately  its  theoretical  appaj‘|P» 
i  Vtiifl  function.  •  '  VvM ' ... 

0  *  , 

i(l)  Crating  Spectrometer  .  . 

Consider  a  grating  w>th  n  el.it:  infinitely  thin  and  th-  ^  f 
...vibrat  ien.*  issuing  from  th  co  n  rllts  in  the- direction -01  ”• 

(rig.  101a);  p(  O  is  then  a  Dirac' cir, tv i Lotion  of  interval 
b  -  a(sin  tf  -  sin  <*’),  a  b  ing  the  ditancu  between  fy»o  lines; 
of  the  grating,  1  hired  by  a  aqupre  guise  function  of  width  v 
L  -  (ji-1)  ;  (Fir.  lOib).  If  the 'grating  receive  fan  impulse' 
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r  Fig.  101.  (a)  Schematic  Diagram  of  a 

Grating  Spectrometer;  (b)  Impulse  Re¬ 
sponse  of  a  Grating  Spectrometer;  (c) 
Autocorrelation  Function  of  the  Impuls# 
Response  of  a  Crating  Spectrometer. 


.....of  light  it  transmits  n  impulse  sin  the  direction  O'.  Cp(*) 

is  a  Dirac  distribution  of  interval  h,  modulated  by  a  triangu¬ 
lar  function  A (0)  which  is  aero  c.utside  the  interval  •  L»  *  L 
(Fig.  101c)..  Hence  the^expression  of  the- apparatus-function  ■ 


i. 


F<e)  -  |0(«)l*  -  *coely*>J 


we'  in 
£  L* 


sin  » 

t 

u‘ 

n  '  **) 

*  L 

7* 

■% 

A  .. 


•  .  .  f' 

(2)  rabrv-Fcrot  Spectrometer 

Let  c  be  the  distance  separating  the  two  plates,  r  the 
amplitude  reflection  factor  of  .the  surfaces  of  the  plates 
<#ig.  102a).  in  this.caae  (Fig.  102fc), 

p(l)  *  a1*/*  R^O)  «  «xp  £  Log -^)  R*<0  » 
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Fig.  102.  (a)  ichcmatir  Ujagram  oi  a  I’ahry- 

Pcrot  bp**  v  trine  ter;  (b)'  Impulse  Response  ol* 
a  Fabry-l’croC  bpcot  rout*.*  tc:r;  (c)  Autocorrcla-  j  *h  f;‘-  /' j;.r '  i  ,< 
tioji  function  of  the  Impulse  Response  of  a  J  . .  ;  .i 
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Fabry- I'c  rot  Upootromet*.^;  ( J)  Apparatus ■» 
Function  of  a  Fabry-Perot  bpc^VJHiiaiW#;. 
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This  is  the  superposition 
of  the  resonsance  functions 
which  repeat  at  multiples  of  ^ 
l/h  ,  This  is  an  Airy_  function 
whose  width  at  mid-height  is 
do  -  (l/*>(l-r*)/»r  (Fig.  102d). 

(3) 


Fig.  103. 

a  Prism. 


i  t 

Impulse  Response  of 


Prism  Spectrometer 
t£>rd  Rayleigh  [65]  studied  the  propagation  of  an  impulse 
ir  a  disp.rsive  medium  andjnade^  an  application  of  it  in  the 
case  of  a  prism.  He  showed  that  the  percussional  response  of 
a  prism  in  a  given  direction  is  a  train  of  sinusoidal  waves 
of  frequency  i»l  and  of  length  L  *  t  (dn/dX)  Xj  with  *  c/v 

if  *  is  the  length  of  the  base  of  the  prism  (Fig.  103)  and 


i  -»’fc 

’t 


I  sin  rL(ff 
«•>  -  |GC.)I*  -  {- • 


In  this  case  Cp(4)  -  A(6). 
zero  outside  of  the  domain 


A(5 )  is  a  triangular  function, 

\  L,  +  L.  The  function  ( < )  has 
an  infinitely  small  interval  because  there  if  an  infinity  of 
rays  which  iufn. e-fere.  The.  spectrum  which  one  obtains  with 


v  •  ^ .  th$ > prism  ?  unique. 


4  .  FOURii:K  Til  Ah  5F0RI  i  ATI  ON  METHOD 

In  a  two-wave  interferometer,  the  incident  vibration  is 
divided  into  two  vibrations  ol  equal  amplitude.  If  cue  instru¬ 
ment  receives  an  impulse ,  it  transmits  two  staggered  impulses 
in  the  time  r|  *  »/c  (Fig.  104a);  Cp(4)  has  ar  particularly 
simple  form.  vt  is  composed  of  three  verticals:  two  at 
A  *  a  ft  of  hr.  ig nt  1  and  another  of-  double  height  at  4  *  0 
(Fig.  104b)  in  i  in  h  a  way  that 

♦(*)  -  j  ^*x(4)  C  (4)  d4  -  2  *rx(0>  +  2  fxfO  . 
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Fig.  104.  (a)  Impulse  Response 

of  a  Micheison  Inter! erome ter; 

(b)  Autocorrelation  Function  of 
the  Impulse  Response  of  a  Michel** 
son  Interferometer. 


*• 


Here  4»{fc)  is  the  sum  of  two  particular  values  of  This 

is  not  enough  to  make  a  Fourier  transformation.  In  fac.  ) 
is  composed  of  a  constant  part  <Z?X(0>  of  a  variable  part  which 
gives  precisely  )  when  one  varies  h  linearly  as  a  function 

of  time.  .  ; 

To  reconstruct  the  spectrum  from  the  interf crogram,  one 

should  be  able  to  treat  the  interf crogram  exactly  as  the  clas- 
~slcal  spectrometers  do,,  by  forming  for  the  calculation  of  the 
spectral  density  corresponding  to  each  wave  number  the  sum  of 
a  particular  series  of  discrete  and  equidistant  values  ol  the 
interf crogram.  One  should  be  able,  of  course,  to  choose  A<6) 
at  will  in  such  a  way  that  one  could  obtain  exactly  the  same 
theoretical  apparatus-function  as  with  a  classic  spectrometer, 
with  the  two  series  of  peaks  mixed  together  and  a  pc^k  con¬ 
stantly  centered  on  the  zero  frequency.  But  practically,  it 
is  much  simpler  to  choose  a  single  series  of  values  ,ot  ^(6) 
and  numerically  to  take  the  Fourier  transformation  by  varying 
the  cosines  as  explained  in  Chapter  II. 

5.  CONCLUSION,  ^  "  , 

All  spectrometers  take  the  F.T.  of  the  even  part  *x(*> 
of  the  autocorrelation  function  of  the  luminous  vibration. 

The  resolution  is  limited  because  ?x<*>  ^  known  only  between 
0  and  L,  the  maximum  phase  difference  between  the  interfering 
beams.  When  the  disperser  is  a  prism,  the  F.T.  .of  .-*.*>  A!» 
mat  c  by  using  a  continuous  saqu-ncc  of  values  between  0  and  L, 
and  the  spectrum  obtained  is  unique.  In  the  other  classical 
spectrometers  the  F.T.  is  obtained  as  the  sum  of  discrete 
equidistant  values  of  ?x(*)i  thre  apparatus-function  is  therefore 
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composed  of  a  series  of  peaks,  each  one  of  them  being  the  F.T. 
of  the  weighting  function  A(5)  appropriate  for  the  instrument 
which  multiplies  In  the  Fourier  transformation  method 

A  j 

it  is  C*^(6)  which  is  measured.  One  can  either  make  the  F.T. 
by  using  a  continuous  sequence  of  values  of  **x(6)  as  for  the 
prism,  in  which  case  the  spectrum  obtained  is  unique;  or  nu¬ 
merically  calculate  the  F.T.  using  discrete  equidistant  values 
of  If  one  uses  a  different  .series  of  values  for  each 

wave  number,  since  one  can  choose  A(S)  at  will,  one  can  obtain 
exactly  the  same  apparatus- function  as' with  any  classical 
spectrometer.  When  the  Fourier  transformation-  is  made  based 
on  a  single  series  of  values  of  the  apparatus-f unction 

is  composed  of  two  series  of  peaks  which  are  displaced  in 
opposite  directions  during  the  exploration  and  the  limit. jcu 
resolution  is  fixed  from  one  extremity  of  the  spectrum  to  the 
other. 

t - 
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_  VIII,  CONCIUSIOH 

In  conclusion,  It  seems  to  us  that  the  &pectroraetric 
method  using  Fourier  transformations  Is,  In  t lie.. infrared,  the 
most  powerful  of  the  spectrometric  methods .  Indeed  it  unit  *■- 
the  fundamental  advantage  of  the  classical  interferometric 
vethods,  namely  the  large  useable  solid  angle,  and  that  of  tiir 
"multiplex*  methods ,  in  which  the  whole  measuring  time  is  us*.o 
tor  the  determination  of  each  spectral  element.  Its  superiority 
over  the  other  methods  Increases,  with  the  number  of  spectral 
elements  to  be  studied.  It  permits  us  to  treat  in  the  hear 
infrared  two  problems  at  least  one  of  which  would,  have  been 
impossible  to  resolve,  by  other  methods;  In  both  cases,  the 
gain  which  o^n  be  attained  is  expressed  as  a  considerable  in- 
croase.in  rmlolotlou,  , 

‘  But  its  use,  as  sood  as  the  desired  rejiolutions  pas 3 
^several  thousands,  poses  some  problems  which  are  not  yet  re* 
Solved*,  obtaining  a  large  ph#a»  oiffct-cnco  without  putting'  ' 
the  i.nterferometer  out  of  order  and  calculating  the  F  .?•  of 
£he  interferogram.  In  order,  to  perform  this  last  operation, 

■'■T  "r*  «t  tKe  present  tiSh  two  devices  are  used:  harmonic  ran*lytemB 

i  ^ i f‘*  .  '  .  ^  #  vy.  ;  *.  j 

which  do  not  allow  one  to  exceed  re  solutions  of  the  order  ol 


A  ■ 

r 

-** 


K  'k-’W 


.  -  f 


i  ' 


£ , 
/■>  ■ 


100and4\git*l  computers. 

\  calculation  t inu.  for  a  spectrum  by  a  numerical. 

Pout  lev  transformation  is  proportional  to  the  number  N  of 
,  input  points.  In  the~ab*«!»cc  of  noise.  It  will  be  c.'  the 
Girder  of  the  number  K  of  spectral  elements  *i>  rho  r.pcctrv^.ii 
t(i  be  studied.  The  .*  t  udv  of  eh-  f  a\  U»j  .v  <n  wtiiu  i*oi:  ■■  in' 

J.C'.-.s  tuav  in  or  'i.v  t 7T~<  l 


i  :i  p. 


the  calculate!  .  p.  'f .  ch  n-  oo 
the  maximum  Vinnal,  no  loo  ratio  in  the  spectrum,  if  thn  inter 
ferogram  has  b-.eri  r»-c<.  r.v,  i  through  a  low- pas *..  filter,  N  can 


St  'yory. much  larger  chair. ,14;  Several  ways  of  recording  the 
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inter*  -Togran  or  of  calculating  •  How  one  to  obtain  the  spectrum 
under  ♦  ho  best  conditions  and  ..t  the  same  time  to  use  a  number 
N  near  1‘  to  make  the  Fcurier  1  ransform:  static  recording, 
method  01  changing  tl.c  frequency,  numerical  filtering.  What¬ 
ever  rue  the  d  ,‘s  used  wntn  treating  a  problem  with  a  great 
number  of  spectral  elements  by  Fourier  transformation,  one 
meets  two  difficulties:  the  first  which  deals  with  the  need 
to  pick  oi f  a  great  number  of  points  from  the  interf  erograra 
can  be  mitigated  by  using  a  converter  system  which  records 
directly  on  punched  cards  the  necessary  values  of  the  inter- 
ferogram.  But  the  second  which  is-  concerned  with  the  long 
calculation  time  and  the  delay  whi^h  occurs  between  record¬ 
ing  the  inter ferogram  and  obtaining  the  spectrum  is  inevitable 
The  regular  use  of  this  method  for  problems  at  even  medium 
resolution,  hence  appears  difficult  as  long  as  ohe  does  not 
have  a  high-speed  calculator  available  for  Fourier  transforms 
which  allows  one  to  obtain  the  spectrum  very  rapidly  after 
recording  the  inter ferogram. 

In  concluding  I  must  express  my  appreciation  to  Professor 
Jacquinot  who  started  me  on  this  research. 
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0.  Parodi  for  the  part  they  have  taken  in  the  experimental 
work  and  for  numerous  and  fruitful  discussions,  P.  Giacomo  for 
the  help  he  gave  me  at  the  time  of  the  construction  of  the 
interferometer,  R.  Chabbal,  all  the  technicians  and  workers 
of  the  Ainrf  Cotton  Laboratories,  and_J.  Arsac  for  the  advice 
he  gave  me  in  the  study  of  questions  concerning  noise. 
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